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ABSTRACT 

Procedures  to  Predict  Vertical  Differential  Soil  Movement 

For  Expansive  Soils.  (December  1 997) 

Donald  David  Naiser,  Jr.,  B.S.,  Texas  A&M  University 

Chair  of  Advisory  Committee:  Dr.  Robert  L.  Lytton 

Damage  to  lightly  loaded  structures,  paving  and  service  piping  in  areas  of 
expansive  clay  soils  has  occurred  throughout  the  world.  The  cause  of  this  damage  has  been 
the  inability  to  accurately  model  expansive  soil  movement  so  that  foundations  are 
adequately  designed  to  withstand  the  movement.  The  amount  and  rate  of  differential  soil 
movement  for  expansive  soils  is  due  to  a  combination  of  soil  characteristics,  namely: 
suction  compression  index,  unsaturated  permeability,  and  diffusivity.  Currently, 
geotechnical  engineers  run  tests  to  measure  the  soil  properties  required  to  estimate 
differential  soil  movements.  However,  there  seems  to  be  apprehension  toward  attempting 
these  soil  movement  calculations  due  to  the  perceived  complexity  of  the  calculations  or  a 
simple  lack  of  understanding  of  the  theory.  The  procedures  delineating  the  step  by  step 
process  used  to  calculate  suction  profiles  and  volume  strains  of  expansive  soils  is 
presented.  These  procedures  include  the  methodology  to  predict  soil  heave  and  shrink 
underneath  shallow  foundations  which  generate  maximum  center  lift  and  maximum  edge 
lift  slab  distortion  modes.  The  main  contributions  of  this  research  are:  equations  and 
procedures  to  calculate  the  equilibrium  suction  profile  and  depth  to  constant  suction  for  a 
particular  soil  profile  and  location,  equations  to  calculate  the  horizontal  velocity  flow  of 
water  in  unsaturated  soils,  the  methodology  to  predict  differential  soil  movement  shortly 
after  a  slab  has  been  constructed  and  before  the  soil  under  the  slab  has  reached  an 
equilibrium  moisture  content,  and  the  procedures  to  apply  differential  soil  movement 
theory  to  soil  profiles  with  multiple  layers  and  moisture  effect  cases  to  be  used  for  shallow 
foundation  design. 
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(4,(Zj)  Percent  passing  2  micron  for  the  ith  soil  increment  (pg  21) 

vs  Percent  passing  number  200  sieve  for  layer  s  (Appendix  A-84) 

v(Zj)  Percent  passing  the  number  200  sieve  for  the  i"1  soil  increment  (pg  21) 

7i  Osmotic  suction  (pg  1 3) 

7uF(Zi)  Final  osmotic  suction  value  a  depth  zi  (pg  19) 

^(Zj)  Initial  osmotic  suction  value  a  depth  z{  (pg  19) 

0 .  Volumetric  water  content  at  depth  z^  corresponding  to  hdiy  (pg  32) 

0m  Volumetric  water  content  at  depth  zm  corresponding  to  h,,,  (pg  32) 

6r  Residual  volumetric  water  content  for  the  soil  at  depth  z^  (pg  33) 

0s  Saturated  volumetric  water  content  for  the  soil  at  depth  z^  (pg  33) 

0wl  Volumetric  water  content  at  depth  z^  corresponding  to  h^,,  (pg  32) 

p  Shape  factor  in  Gumbel's  probability  equation  (pg  28) 

aF(Zj)  Applied  octahedral  normal  stress  at  depth  Zj  (pg  19) 

ct,  Stress  level  below  which  there  is  no  suppression  of  suction  volume  strain  by 

overburden  pressure  (pg  1 9) 

\\)  Total  suction  (pg  13) 
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LIST  OF  DEFINITIONS 


Cation  exchange 
Cation  exchange  activity 
Cation  exchange  capacity 


Center  lift  distortion  mode 


Activity  ratio  The  ratio  of  the  plasticity  index  to  the  percentage  clay 

content  (Skempton  1955).  Percent  clay  is  defined  as  the 

percentage  less  than  2m  of  the  soil  material  passing  the  U.S. 

No.  200  sieve  (Pearring  1968). 

The  interchange  between  a  cation  in  solution  and 

another  cation  on  a  surface  active  material  (Holt  1969). 

The  ratio  of  the  cation  exchange  capacity  to  the  percentage 

clay  content  (Pearring  1968). 

The  total  amount  of  exchangeable  cations  that  a  soil 

is  capable  of  absorbing,  measured  in  milliequivalents  per  1 00 

grams  of  soil  (Mojekwu  1979). 

This  is  a  slab-on-ground  distortion  mode  resulting  from  the 

soil  at  the  perimeter  being  more  dry  (perimeter  soil  may  be 

shrinking)  than  the  soil  beneath  the  slab  (soil  beneath  the 

slab  is  either  at  equilibrium,  shrinking  at  a  rate  slower  than 

the  exterior  soil,  or  heaving).  Refer  to  Fig.  2  (Castleberry 

1974  and  Mitchell  1980). 

The  soil  particles  smaller  than  2\x  which  are  derived  from 

the  chemical  decomposition  of  rock.  It  consists,  chiefly,  of 

clay  minerals  but  small  amounts  of  quartz,  feldspar,  organic 

matter,  soluble  salts  and  amorphous  materials  are  also 

present  (Lambe  and  Whitman  1969). 

Depth  of  available  moisture    The  maximum  depth  of  moisture  available  for  use  by 

transpiring  vegetation,  which  is  stored  within  the  soil  zone 
up  to  the  depth  to  constant  suction  (Ga.y  1994). 

Depth  to  constant  suction        The  depth  in  a  soil  profile  to  which  there  is  no  longer  a 


Clay 
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Edge  lift  distortion  mode 


Edge  moisture  variation 
distance 


significant  seasonal  suction  change  (the  seasonal  change  in 
suction  is  less  than  0.2  pF  (Lytton  1997)). 
This  is  a  distortion  mode  resulting  form  the  soil  at 
the  perimeter  being  more  wet  (perimeter  soil  may  be 
heaving)  than  the  soil  beneath  the  slab  (soil  beneath  the  slab 
is  either  at  equilibrium,  heaving  at  a  rate  slower  than  the 
exterior  soil,  or  shrinking).  Refer  to  Fig.  3  (Mitchell  1980). 
For  soil  under  an  impervious  membrane,  the  edge 
moisture  variation  distance  is  the  distance  from  the  edge  of 
the  barrier  to  the  point  where  soil  suction  does  not  vary  with 
season  moisture  variation  (Jayatilaka  et  al.  1992). 

Equilibrium  suction  profile     A  suction  profile  generated  by  decreasing  the  equilibrium 

suction,  h,,,  a  centimeter  in  suction  for  every  centimeter 
increase  in  depth  (McQueen  and  Miller  1968). 
Any  impervious  barrier  that  is  attached  and  sealed  to 
the  edge  of  the  foundation  and  extends  some  distance,  dhbar 
horizontally  from  the  edge  of  the  slab.  A  sidewalk  attached 
and  sealed  to  the  edge  of  a  foundation  is  considered  a 
horizontal  barrier.  Refer  to  the  sidewalk  shown  in  Appendix 
C-132. 

The  water  content,  expressed  as  a  percentage  of  the  weight 
of  the  dry  soil,  at  the  boundary  between  the  liquid  and  plastic 
states.  The  water  content  at  this  boundary  is  defined  as  the 
water  content  at  which  a  groove,  1mm  wide,  cut  in  the  soil 
sample  with  a  standard  grooving  tool,  closes  for  a  length  of 
Vi  inch  when  a  brass  dish  containing  the  soil  is  dropped  25 
times,  at  a  rate  of  2  drops  per  second,  through  a  distance  of  1 
cm  on  a  standard  hard  rubber  base  (Mojekwu  1979). 

Mean  principal  stress  The  slope  of  the  volume-total  stress  curve  in  the  plane  of 


Horizontal  barrier 


Liquid  limit 


XX 


Plasticity  index 


Post  construction 


compression  index  zero  suction  (Lytton  1994). 

Plastic  limit  The  water  content,  expressed  as  a  percentage  of  the  dry 

weight  of  the  soil  at  which  the  soil  becomes  crumbly  and 
ceases  to  be  plastic.  The  water  content  is  defined  as  the 
water  content  at  which  the  soil  just  begins  to  break  apart  and 
crumble  when  rolled,  by  hand,  into  threads  one-eighth  of  an 
inch  in  diameter  (Mojekwu  1979). 
The  difference  between  the  greater  moisture  content 
of  the  liquid  limit  and  the  less  moisture  content  of  the  plastic 
limit.  It  is  defined  as  the  range  in  water  contents  through 
which  the  soil  remains  in  the  plastic  state  (Mojekwu  1979). 
This  refers  to  the  transient  case  for  calculating  volume 
strain  of  expansive  soils  immediately  after  construction, 
before  the  soil  under  the  slab  has  reached  an  equilibrium 
moisture  content. 

This  refers  to  the  case  for  calculating  volume  strain  of 
expansive  soils  after  the  soil  under  the  slab  has  reached  an 
equilibrium  moisture  content, 
the  negative  gage  pressure  relative  to  the  external  gas 
pressure  on  the  soil  water  to  which  a  pool  of  pure  water 
must  be  subjected  in  order  to  be  in  equilibrium  through  a 
semi-permeable  membrane  with  the  soil  water. 
The  absolute  value  of  the  range  between  the  maximum  and 
the  minimum  suction  values  being  considered  in  soil 
suction  profile,  (typically,  U0  =  h^fe)  -  Km(^)  or 

Suction  compression  index     The  slope  of  the  volume-total  suction  curve  in  the  plane  of 

zero  suction  (Lytton  1994). 
Suction  profile  When  the  state  of  suction  is  measured  at  intervals  of  depth 


Post  equilibrium 


Soil  suction 


Suction  amplitude 


XXI 


down  the  soil  profile,  the  resulting  relation  between  suction 
and  depth  is  termed  the  suction  profile  (Mitchell  1981). 
Thornthwaite  Moisture  Index  TMI  is  defined  as  a  climatic  index  that  measures  the  water 

balance  of  an  area  to  determine  if  there  is  a  deficiency  or 
excess  to  ascertain  if  the  climate  should  be  classified  as  arid 
or  humid  (Thornthwaite  1948). 

Any  impervious  barrier  that  extends  some  distance,  d^ 
vertically  downward  from  the  foundation  edge. 
The  total  amount  of  vertical  movement  between  two 
distinct  locations  separated  by  a  characteristic  distance. 


Vertical  barrier 

Vertical  differential 
Soil  Movement 


INTRODUCTION 

When  a  slab  barrier  is  placed  on  an  expansive  soil,  the  immediate  effect  is  that 
evaporation  and  precipitation  cannot  effect  the  soil  beneath  the  slab.  During  subsequent 
weather  events  occurring  to  the  surrounding  soil,  moisture  differentials  develop  between 
the  uncovered  soil  and  the  soil  beneath  the  slab  barrier.  These  differentials  are  largest 
when  the  climate  is  described  by  a  distinct  wet  period  followed  by  a  hot  dry  period.  Unless 
a  foundation  system  is  installed  to  prevent  or  reduce  moisture  movement,  moisture 
redistribution  occurs  within  the  soil  mass  to  reduce  these  differences  in  moisture  content 
(Holland  and  Lawrance  1980).  Damage  caused  by  this  type  of  moisture  differential  is 
considerable  and  warrants  investigation  in  developing  a  more  scientific  and  accurate 
methodology  to  address  these  problems. 

The  ability  to  predict  differential  soil  movement  is  critical  for  engineers  to 
formulate  a  picture  of  what  the  soil/structure  interface  will  be.  This  information  is 
necessary  to  adequately  design  foundations  or  other  ground  supported  structures.  In  slab- 
on-ground  design,  the  principal  interest  is  in  making  an  accurate  estimate  of  the  range  of 
vertical  movement  that  must  be  sustained  by  the  foundation.  The  procedures  contained 
herein  are  methods  for  defining  envelopes  of  maximum  heave  and  shrinkage  so  that 
adequate  slab-on-ground  designs  are  made.  For  highway  and  airport  pavements,  canals, 
and  pipelines,  the  wave  spectrum  of  differential  movements  versus  wave  lengths  are  the 
desirable  design  characteristics.  Structural  floors  suspended  above  expansive  clays  must 
be  provided  with  a  gap  that  exceeds  the  total  expected  heave.  Drilled  piers  (or  shafts)  must 
be  designed  to  resist  simultaneously  a  vertical  movement  profile  and  a  horizontal  pressure 
profile,  both  of  which  change  with  wetting  and  drying  conditions.  Retaining  structures, 
basement  walls,  rip  rap,  and  canal  linings  must  be  designed  to  withstand  lateral 
movements.  Finally,  all  foundations  must  be  designed  against  the  time-dependent  vertical 
and  horizontal  curvature  that  is  generated  by  down  hill  creep  (Lytton  1994). 


This  thesis  follows  the  style  of  the  Journal  of  Geotechnical  Engineering. 


Each  type  of  soil  movement  listed  above  and  their  effects  on  different  types  of 
ground  structures  is  of  sufficient  importance  and  complexity  to  warrant  research  of  its  own. 
Addressed  herein  are  the  effects  of  vertical  differential  soil  movement  on  commercial  and 
residential  types  of  slabs-on-ground.  Vertical  differential  soil  movement  in  these  cases  is 
defined  as  the  total  amount  of  vertical  movement  between  two  distinct  locations  separated 
by  a  characteristic  distance.  The  characteristic  distance  is  described  as  the  edge  moisture 
variation  distance  plus  the  distance  a  moisture  effect  source  is  away  from  the  slab  edge. 
The  model  used  to  predict  differential  movement  of  expansive  soils  is  based  on  constitutive 
soil  relationships  derived  primarily  from  works  developed  by  C.  W.  Thornthwaite  (1948), 
W.  R.  Gardner  (1958),  R.  L.  Lytton  (1973,  1977,  1994,  1995,  1997),  J.  R.  Pearring  (1968), 
J.  H.  Holt  (1969),  E.  C.  Mojekwu  (1979),  J.  L.  Nieber  (1981),  R.  G.  McKeen  (1981,  1992), 
P.  W.  Mitchell  (1980,  1984),  D.  A.  Gay  (1994). 

It  has  been  well  documented  that  the  magnitude  at  which  soil  shrinks  and  swells  is 
not  exclusively  based  on  the  soils  characteristics.  The  magnitude  of  differential  soil 
movement  in  expansive  soils  is  dominated  by  the  following  variables:  soil  characteristics, 
pre-construction  and  post  construction  vegetation,  slope  of  perimeter  surface  grade,  shape 
and  age  of  slab  and  depth  of  soil  moisture  active  zone  (Castleberry  1974,  Mathewson  et  al. 
1975,  Mathewson  et  al.  1980). 

The  methodology  used  to  predict  vertical  differential  soil  movement  involves  using 
soil  suction  profiles  under  the  foundation  and  at  the  edge  of  the  foundation.  These  suction 
profiles  are  developed  using  the  following  variables:  depth,  time,  local  surface  annual 
weather  and  vegetation  conditions,  and  soil  characteristics  such  as  the  suction  compression 
index,  unsaturated  permeability,  and  soil  diffusivity.  Soil  suction  profiles  combined  with 
overburden  stresses  are  used  in  volume  strain  equations  developed  by  Lytton  (1994)  to 
predict  the  volumetric  strain  for  an  incremental  layer  of  soil.  The  incremental  layer 
movements  are  summed  up  and  total  vertical  differential  soil  movement  is  predicted  for  a 
particular  column  of  soil.  These  procedures  are  detailed  in  the  chapter  on  "The  Procedures, 
in  a  Software  Design  Format".  These  procedures  address  foundation  perimeter  moisture 
effect  cases  that  are  either  very  common  with  many  foundation  problems  or  have  a 


significant  effect  on  vertical  differential  soil  movement.  Specifically,  the  foundation 
perimeter  moisture  effect  cases  that  are  addressed  include:  (1)  bare  soil  at  the  surface,  (2) 
grass  at  the  surface,  (3)  flower  bed  at  the  surface  with  a  defined  moisture  effect  zone  depth, 
(4)  trees  at  the  surface  with  a  defined  root  zone  depth,  and  (5)  moisture  effect  zones  that 
have  been  measured  in  the  form  of  suction  profiles.  Included  are  methods  for  predicting  the 
effects  that  vertical  or  horizontal  barriers  have  on  vertical  differential  soil  movement. 
Additionally,  methods  to  predict  differential  soil  movement  shortly  after  construction  and 
before  the  soil  has  reached  an  equilibrium  moisture  content  under  the  foundation  are  the 
highlights  of  this  research  and  are  presented. 

In  "The  Problem"  chapter,  the  theory  and  methodology  used  to  predict  differential 
soil  movement  is  defined  and  described  in  order  by  which  the  theory  was  developed. 
These  bits  and  pieces  of  theory  are  applied  in  step  by  step  processes  described  in  the 
procedures  sections  of  this  document  to  address  all  of  the  different  moisture  effect  cases, 
with  and  without  barriers  imposed. 


THE  PROBLEM 

The  Magnitude  and  Description  of  the  Problem 

Since  early  1950,  damage  to  lightly  loaded  structures,  paving  and  service  piping  in 
areas  of  expansive  clay  soils  has  occurred  throughout  the  world.  Estimates  show  that  one 
in  five  Americans  are  affected  by  this  type  of  damage  while  only  one  in  ten  are  affected  by 
major  floods.  Structural  damage  caused  by  expansive  clay  soils  in  the  United  States 
annually  exceeds  that  caused  by  earthquakes,  hurricanes,  and  floods  combined  (Jones  and 
Holtz  1973.)  Damage  to  residential  and  commercial  structures  caused  by  soil  movement  is 
estimated  at  660  million  dollars,  almost  one-third  of  the  total  monetary  damage  received  by 
all  structures  as  a  result  of  soil  movement.  On  a  state-by-state  basis,  California  and  Texas 
account  for  35%  of  all  expansive  soil  losses  in  the  United  States.  This  becomes  obvious 
from  a  map  of  the  location  of  much  of  the  expansive  soils  in  the  United  States,  see  Fig.  1. 
Much  of  this  soil  is  located  in  highly  populated  regions  of  both  California  and  Texas. 
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Fig.  1.  United  States  Map  of  Expansive  Soils  After  Wiggins  (1976) 


When  it  is  realized  that  most  residences  and  light  commercial  buildings  in  these  two  states, 
which  together  had  430,000  new  home  starts  in  1978,  are  constructed  with  slab-on-ground 
foundations,  the  need  to  develop  a  theoretical  design  approach  based  on  physical  and 
natural  law  for  slab-on-ground  design  in  expansive  soil  areas  becomes  evident. 

The  problem  with  slabs-on-ground  built  on  expansive  soils  becomes  manifested  as 
excessive  structural  distortions  in  response  to  the  subsoil  undergoing  volume  changes 
induced  by  a  change  in  subsoil  moisture.  If  foundations  are  not  of  a  stiffness  to  resist  or 
span  the  gaps  caused  by  differential  soil  movements,  and  if  the  superstructure  is  not  of  a 
suitable  flexibility,  the  foundation  distortions  resulting  from  the  soil  movements  can  be  of 
a  magnitude  to  cause  excessive  cracking  of  walls  and  distortion  of  internal  fittings.  This 
type  of  differential  shrink  and  heave  of  soils  give  rise  to  two  distinct  distortion  modes 
which  cause  damage  to  the  superstructure.    Fig  2  shows  typical  superstructure  damage  that 
would  occur  in  the  case  of  center  lift.  Center  lift  is  a  distortion  mode  resulting  from  the 
soil  at  the  perimeter  being  more  dry  (perimeter  soil  may  be  shrinking)  than  the  soil  beneath 
the  slab  (soil  beneath  the  slab  is  either  at  equilibrium,  shrinking  at  a  rate  slower  than  the 
exterior  soil,  or  heaving.)  Fig  3  shows  typical  superstructure  damage  that  would  occur  in 
the  case  of  edge  lift.  Edge  lift  is  a  distortion  mode  resulting  from  the  soil  at  the  perimeter 
being  more  wet  (perimeter  soil  may  be  heaving)  than  the  soil  beneath  the  slab  (soil  beneath 
the  slab  is  either  at  equilibrium,  heaving  at  a  rate  slower  than  the  exterior  soil,  or 
shrinking.)    The  engineering  profession  lacks  a  universal  method  of  analysis  within  a 
theoretical  framework  that  utilizes  fundamental  properties  of  soil  and  footing  behavior  to 
enable  a  footing  design  that  is  satisfactorily  accurate  without  resorting  to  empirical 
methods  (Mitchell  1980).  Since  the  late  1970's  strides  have  been  made  by  the  Post 
Tensioning  Institute  (PTI)  to  adopt  a  standard  method  for  designing  ground  structures  in 
expansive  soils.  The  Post  Tensioned  Institute  has  adopted  the  methodology  discussed 
herein  and  stipulates  that  post  tensioned  slab  designs  shall  provide  necessary  rigidity  to 
avoid  damage  to  the  buildings  or  houses  when  the  supporting  soil  expands  or  contracts 
(PTI  1996).  However,  a  theoretical  approach  to  predict  vertical  differential  soil  movement, 
as  presented  herein,  has  not  been  widely  adopted  by  today's  practicing  engineers. 
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Fig.  2.  (Mitchell  1980)  Center  Lift  Distortion  Mode  Fig.  3.  (Mitchell  1980)  Edge  Lift  Distortion  Mode 


Expansive  Soil,  on  a  Molecular  Level 

The  Building  Blocks  of  Expansive  Soil 

Why  are  "expansive  soils"  so  expansive?  The  answer  can  be  found  in  the  material 
on  a  molecular  level.  The  magnitude  of  the  surface  area  per  mass  and  the  molecular 
composition  of  expansive  soils  plays  a  major  role  in  the  soil's  ability  to  shrink  and  swell  to 
such  an  impressive  degree.  "Expansive  soils"  are  typically  the  fine  grained  soils  as 
classified  in  the  unified  classification  system.  Fifty  percent  or  more  of  a  soil,  by  weight,  in 
a  sample  that  passes  the  number  200  sieve  is  classified  as  a  fine  grained  soil.  Fine  grained 
soils  are  further  broken  down  into  categories  of  silts,  clays  and  organics  having  either  a 
high  or  low  plasticity.  The  ability  to  break  down  and  classify  a  soil  according  to  the  unified 
soil  classification  system  provides  an  engineer  with  a  wealth  of  information  about  the 
behavioral  characteristics  and  how  a  particular  soil  may  effect  design.  An  in  depth 
explanation  of  soil  characteristics  and  the  unified  soil  classification  system  can  be  found  in 
Lambe  and  Whitman  (1969). 

The  following  clay  structural  unit  explanation  is  idealized  and  simplified  and  are 
presented  to  provide  a  picture  of  the  basic  behavioral  characteristics  that  give  rise  to  clay's 
ability  to  undergo  large  changes  in  volume.  The  structural  units  that  make  up  a  clay 
particle  are  silicate  sheets.  Primarily,  clays  are  made  of  sandwiches  of:  silica  sheets, 


gibbsite  sheets,  and  brucite  sheets.  These  sheets  are  combined  in  two  distinct  groups: 
two-layer  sheet  minerals,  and  three  layer  sheet  minerals.  The  most  important  and  most 
common  two-layer  sheet  mineral  is  the  kaolinite  clay  particle,  see  Fig.  4.  Kaolinite  is  a 
two-layer  mineral  which  consists  of  a  silicate  sheet  bonded  to  a  gibbsite  sheet  by 
hydrogen  bonding  and  secondary  valence  forces.  Another  example  of  a  two  layer  sheet 
mineral  is  serpentine  which  is  formed  by  brucite  bonding  to  silica,  see  Fig.  5. 


WWWW   WW 


6-OH      -6 


4-A1      +12 


-10 


(b) 


Fig.  4.  (Lambe  and  Whitman  1969)  The  Structure  of  Kaolinite 
(a)   Atomic  Structure  (b)  Symbolic  Structure 


In  the  actual  formation  of  these  silicate  sheets  there  may  be  some  isomorphous 
substitution.    Isomorphous  substitution  is  when  one  kind  of  atom  is  substituted  for 


WW  W  WW  w 


6-OH 

-6 

6-Mg 

+  12 

4-0 
2-OH 

-10 

4-Si 

+  16 

6-0 

-12 

z 


(6) 


Fig.  5.  (Lambe  and  Whitman  1969)  The  Structure  of  Serpentine 
(a)   Atomic  Structure  (b)  Symbolic  Structure 


another  atom  that  is  very  similar  in  structure  and  size.  In  other  words  aluminum  may  be 
substituted  for  silicon  if  there  is  a  shortage  of  silicon  and  an  abundance  of  aluminum 
during  the  formation  stage  of  the  clay  mineral.  This  sort  of  substitution  may  give  the 
mineral  a  net  charge  deficiency  and  may  cause  slight  distortion  in  the  crystal  lattice  due 
to  the  substituted  atom's  size  difference.  Isomorphous  substitution  that  occurs  giving 
the  clay  mineral  a  large  negative  charge  will  give  the  mineral  the  capacity  to  attract 
polar  charged  water  molecules.  A  double  layer  of  molecules,  such  as  water  molecules, 
will  try  to  form  around  a  clay  mineral  to  neutralize  the  negatively  charged  clay  particle. 
This  great  affinity  for  water  due  to  isomorphous  substitutions  provides  an  explanation 
as  to  why  clay  minerals,  like  montmorillonite,  can  have  large  volume  changes  at  high 


volumetric  water  contents.  In  Fig.  6  the  clay  particles  are  shown  with  the  fully 
developed  double  layers  they  would  have  in  pure  water. 


(a) 


1  I 


(b) 


Fig.  6.  (Lambe  and  Whitman  1969)  Soil  Particles  With  Water  and  Ions 
(a)  Sodium  Montmorillonite  (b)  Sodium  Kaolinite 

In  kaolinite  there  is  a  very  small  amount  of  isomorphous  substitution.  However 
in  the  three  layer  sheet  clay  minerals  there  is  much  more  isomorphous  substitution. 
Three  layer  sheets  are  formed  by  placing  a  silica  on  the  top  of,  and  on  the  bottom  of, 
either  a  gibbsite  or  brucite  sheet.  Fig.  7  shows  the  mineral  pyrophyllite  made  of  a 
gibbsite  sheet  sandwiched  between  two  silica  sheets.  Fig.  8  shows  the  structure  of  the 
mineral  muscovite,  which  is  similar  to  pyrophyllite  except  that  there  has  been 
isomorphous  substitution  of  aluminum  for  silicon  in  muscovite.  The  net  charge  created 
by  this  substitution  is  balanced  by  potassium  ions  which  serve  to  link  the  three  layer 
sandwiches  together.  The  most  common  three-layer  minerals  in  soil  are 
montmorillonite  and  illite.  Montmorillonite  is  similar  to  pyrophyllite  with  the 
exception  that  there  has  been  isomorphous  substitution  of  magnesium  for  aluminum  in 
the  gibbsite  sheets.  Due  to  the  isomorphous  substitution  the  illite  mineral  has  a 
moderately  high  capacity  for  volume  change  and  montmorillonite  has  a  very  high 
capacity  for  volume  change  (Lambe  and  Whitman  1969). 
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Fig.  7.  (Lambe  and  Whitman  1969)  The  Structure  of  Pyrophyllite 
(a)  Atomic  Structure  (b)  Symbolic  Structure 
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Fig.  8.  (Lambe  and  Whitman  1969)  The  Structure  of  Muscovite 
(a)  Atomic  Structure  (b)  Symbolic  Structure 
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Gilgai 

Clay  soil  left  to  interact  with  surface  climate  conditions  will  develop  an  irregular 
surface  microrelief  pattern  called  gilgai,  which  in  colloquial  terms  is  referred  to  as  "hog 
wallows."  This  phenomenon  is  recognizable  in  areas  where  highly  expansive  soils  are 
known  to  exist  undisturbed  for  extended  periods  of  time.  These  gilgai  may  be 
experienced  while  taking  a  very  undulating  ride:  up  the  fairway  of  an  old  golf  course, 
during  farming  evolutions  in  a  native  hay  field,  or  down  an  old  country  road.  As 
illustrated  in  Fig.  9  from  Spotts  (1974),  the  pattern  develops  in  seven  stages.  First,  the 
sedimentary  clay  is  deposited.  Following  deposition,  drying  occurs,  producing  surface 
shrinkage  cracks.  Further  desiccation  results  in  the  formation  of  major  cracks 
penetrating  far  into  the  clay.  The  extent  of  penetration  and  spacing  are  determined  by  a 
number  of  factors  of  which  the  degree  and  rate  of  drying  and  the  clay  type  are  the  most 
significant.  Subsequent  development  of  microrelief  occurs  as  the  clay  is  wetted  and 
dried  by  cyclic  variation  of  the  climate.  Because  water  penetrates  deeper  at  the  major 
cracks,  greater  swell  occurs.  As  the  next  cycle  of  drying  proceeds,  the  higher  areas  dry 
first,  causing  the  major  cracks  to  reappear  in  approximately  the  same  locations.  Thus, 
over  a  period  of  several  years,  an  undulating  surface  pattern  called  gilgai  develops  on 
the  soil  surface  (McKeen  1981).  This  phenomenon  is  important  to  realize  since  it  is  a 
natural  occurrence  of  differential  soil  movement  that  must  be  understood  by  engineers 
so  that  appropriate  models  of  this  type  of  soil  behavior  can  be  properly  applied  to 
design.  There  are  important  parameters  to  consider  when  addressing  volume  strain  that 
causes  the  gilgai  phenomenon,  such  as:  (1)  The  depth  at  which  soil  interacts  with  the 
climate,  known  as  the  depth  to  constant  suction,  2^.  (2)  The  climate,  or  the  annual 
wetting  and  drying  cycles  at  the  surface.  (3)  The  type  of  surface  conditions,  such  as; 
vegetation  growth,  bare  soil,  slope,  or  barriers  of  some  type.  (4)  The  soil  characteristics, 
namely:  suction  compression  index,  yh(Zj),  unsaturated  permeability,  K0 (zj )lh0 1 ,  and  soil 
diffusivity,  a(z;).  The  roles  that  each  of  these  parameters  play  in  predicting  differential 
soil  movement  are  presented. 
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Fig.  9.  (Spotts  1974)  Gilgai  Stages 
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Suction,  Matrix  and  Osmotic 

Soil  suction  is  the  measure  of  a  soil's  affinity  for  water.  Soil  suction  is  an 
energy:  water  always  moves  down  an  energy  gradient  and  it  takes  energy  to  cause 
volume  strain  in  any  soil.  Soil  scientists  and  arboriculture  scientists  refer  to  soil  suction 
in  terms  of  water  potential.  As  defined  by  the  1960  Moisture  Equilibria  Symposium, 
soil  suction  is  the  negative  gage  pressure  relative  to  the  external  gas  pressure  on  the  soil 
water  to  which  a  pool  of  pure  water  must  be  subjected  in  order  to  be  in  equilibrium 
through  a  semi-permeable  membrane  with  the  soil  water.  Suction  is  the  term  used 
principally  by  engineers  for  the  thermodynamic  quantity,  Gibb's  free  energy  which  is 
inherently  negative,  as  seen  in  (1)  below,  and  generates  tension  in  the  pore  water 
stretching  between  soil  particles  (Lytton  1994).  Soil  suction  is  equal  to  0.0  when  the 

R-T'      (  K\ 

h  = -liJ—  (1) 

m-g       U00/  v  ' 

relative  humidity  is  equal  to  1 00%.  A  relative  humidity  value  less  than  1 00%  in  a  soil 

would  indicate  the  presence  of  suction  in  the  soil. 

Components  of  Soil  Suction 

Soil  suction  is  quantified  in  terms  of  the  relative  humidity  (1)  and  is  commonly 
called  "total  suction."  "Total  suction,"  vy  ,  has  two  components,  namely:  matrix  suction 
(u.),  which  is  due  to  the  attraction  of  water  to  the  soil  particle  surfaces,  and  osmotic 
suction  (tu),  which  is  due  to  dissolved  salts  or  other  solutes  in  the  pore  water  (2). 

\\f  -  U.  +  71.  (2) 

Matrix  suction  will  have  the  most  effect  on  the  volumetric  changes  in  soil  due  to  the  fact 
that  this  component  varies  with  the  changes  in  the  moisture  environment.  Osmotic 
suction  will  not  change  unless  there  is  a  change  in  the  salt  concentration  within  a  given 
soil  layer. 

Typical  Suction  Levels  in  the  Field 

Fig.  10  illustrates  the  suction-vs-volumetric  water  content  curve  for  a  natural  soil  under 
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wetting  and  drying  conditions.  A  common  measure  of  suction  is  the  pF-scale,  in  which 
pF  is  defined  as  (3),  where  p  refers  to  the  logarithmic  value  analagous  to  pH,  and  the  F 
refers  to  free  energy. 

Suction  in  pF  =  log]0|  h  |  (3) 

Where  |  h  |  equals  the  magnitude  of  suction  in  centimeters  of  water,  a  positive  value. 
The  various  suction  levels  corresponding  to  field  cases  are  listed  below: 
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Fig.  10.  (Lytton  1994)  Suction  vs  Volumetric  Water  Content  Curve 


Field  capacity  (hfc  =  2.0  pF) 

Wet  Limit  for  clays  (hwl  =  2.5  pF) 

Plastic  limit  (3.5  pF) 

Wilting  point  of  plants  (hdry  =  4.5  pF  for  vegetation) 

Tensile  strength  of  confined  water  (5.3  pF) 

Air  dry  at  50%  relative  humidity  (hdry  =  6.0  pF  for  bare  soil) 

Oven  dry  (7.0  pF) 
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Field  capacity  is  an  agricultural  engineering  term  used  to  define  the  maximum  amount 
of  water  that  a  soil  can  hold  in  the  field.  Fig.  10  shows  field  capacity  to  be  a  value  of 
2.0  pF,  but  actual  field  values  never  quite  drop  below  a  value  of  2.5  pF  for  clays  (Lytton 
1995  and  Mitchell  1980).  Soil  suction  is  a  term  that  measures  energy  potential  and  is 
directly  related  to  the  soil's  potential  for  volume  change.  This  information  is  referred  to 
repeatedly  throughout  this  document. 
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THE  THEORY  OF  DIFFERENTIAL  SOIL  MOVEMENT 

Some  notable  breakthroughs  have  been  made  pertaining  to  the  topic  of 
predicting  differential  soil  movement  and  are  described  in  this  chapter.  The  real 
challenge  in  geotechnical  engineering,  especially  in  the  area  of  expansive  soils,  is  to 
develop  a  consistent  body  of  theoretical  mechanics  which  adequately  explains  the 
behavior  of  the  soil  mass  and  to  develop  methods  of  measuring  the  relevant  material 
properties  critical  for  making  appropriate  design  decisions. 

In  the  1970's  the  first  attempt  was  made  to  design  slabs-on-ground  using  a 
scientifically  based  theoretical  approach.    In  1973,  Lytton  and  Woodburn  prepared  a 
report  which  chronicles  a  successful  school  building  foundation  design.  The  school 
foundation  design  procedure  considered  differentially  expanding  and  shrinking  clay  due 
to  the  changing  suction  of  the  soil  beneath  the  foundation.  The  suction  profiles  of  the 
soil  were  predicted  using  a  rational  method  which  relied  upon  suction  profiles  measured 
on  samples  taken  on  the  site  at  different  times  in  the  year.  These  suction  values  were 
used  to  predict  the  differential  soil  movement,  and  the  foundations  were  designed  based 
on  the  distortion  modes  caused  by  the  magnitude  and  location  of  the  expansive  soil 
movements.  Lytton  noted  that  this  differentially  expanding  and  shrinking  effect  caused 
"mounds"  to  form  beneath  the  slab.  The  foundation  interacts  with  these  mounds 
pressing  down  on  the  high  spots  and  bridging  the  low  spots.  The  foundation  design 
procedures  considers  the  limit  cases  in  which  the  mounds  produce  two  types  of 
distortion  patterns:  edge  lift  pattern  and  center  lift  pattern,  as  shown  in  Fig.  1 1 .  These 
two  distortion  modes  were  used  to  design  the  slabs  so  that  the  loads  of  the 
superstructure  would  be  carried  across  a  span  generated  by  edge  lift,  or  carried  in  a 
cantilevered  fashion  out  to  the  edge  of  the  foundation  in  the  case  of  center  lift. 
Woodburn  recorded  the  performance  of  the  foundation  systems  after  the  designs  were 
built  and  the  school  was  put  into  service.  Lytton  and  Woodburn  reported  post 
construction  edge  shrink  values  to  be  2.75  inches  which  compared  favorably  with  the 
predicted  value  of  3.1  inches.  Just  as  Lytton  predicted,  the  worst  case  differential 
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movement  occurred  in  the  first  dry  season  subsequent  to  the  construction  of  the  slab 
during  the  wettest  time  of  the  year.  In  this  post  construction  case  the  soil  was  "capped' 
and  the  moisture  content  was  preserved  in  a 


Edge  Heave 


Center  Heave 


Fig.  11.  (Lytton  1973)  Edge  Heave  and  Center  Heave  Distortion  Modes 


state  of  heave  under  the  slab  while  the  perimeter  soil  shrank  as  it  dried.  This  post 
construction  case  produced  the  worst  case  center  lift  mode.  Subsequent  weather  cycles 
did  not  produce  such  extreme  cases  of  differential  soil  movements  as  compared  to  the 
post  construction  case.  However,  subsequent  weather  cycles  did  produce  differential 
soil  movements  as  predicted.  The  model  to  predict  this  type  of  post  construction 
differential  soil  movement  is  presented  in  the  Post  Construction  section  of  this  chapter. 
This  design  and  research  venture  documented  by  Lytton  and  Woodburn's  1973 
paper  marked  the  beginning  of  a  slab-on-ground  design  method  that  has  been  adopted 
by  the  Post  Tensioning  Institute  as  the  standard  for  design  of  post-tensioned  slabs-on- 
ground  on  expansive  soils  (PTI  1980,  1996).  However,  the  objective  of  this  research  is 
to  consider  the  methodology  behind  predicting  the  differential  soil  movement  used  to 
design  such  slabs. 

Soil  Volume  Strain 

Movements  in  expansive  soils  are  generated  by  changes  of  suction  which  are 
brought  about  by  the  entry  or  loss  of  moisture  or  change  in  concentration  of  dissolved 
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salts  in  the  pore  water.  The  volume  change  that  accompanies  the  change  of  suction  (and 
water  content)  depends  upon  the  total  stress  states  that  surround  the  soil.  Within  a  soil 
mass,  a  decrease  in  the  magnitude  of  suction  results  in  an  increase  of  water  content. 
The  volume  of  the  soil  also  increases  unless  the  surrounding  pressure  is  sufficient  to 
restrain  the  swelling.  A  conceptual  graph  of  suction-vs-volume  can  be  drawn  using  the 
relations  of  each  to  water  content.  This  is  illustrated  in  Fig.  12  on  the  plane 
corresponding  to  zero  pressure.  A  similar  graph  can  be  drawn  relating  pressure  (total 
stress)  -  versus  -  volume  on  the  plane  corresponding  to  zero  suction.  The  simultaneous 
change  of  the  magnitude  of  suction  (decrease)  and  pressure  (increase)  results  in  a  small 
change  of  volume,  following  the  path  from  Point  A  to  Point  C  on  the  pressure  -  suction 
-  volume  surface.  The  magnitude  of  suction  decreases  from  Point  A  to  Point  B  while 
the  pressure  increases  from  Point  B  to  Point  C.  The  volume  change  process  can  be 
viewed  as  the  net  result  of  two  processes; 

a.  Increase  of  volume  from  A  to  B  at  constant  mechanical  pressure  or  total 
stress. 

b.  Decrease  of  volume  from  B  to  C  at  constant  suction. 


Volume 


Mean  Principal  Stress 
(Log  Scale) 


Matrix  Suction 
(Log  Scale) 


Fig.  12.  (Lytton  1994)  Mean  Principal  Stress-Volume-Matrix  Suction  Diagram 
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For  small  increments  of  volume  change  on  this  surface,  the  volume  strain,  is  linearly 
related  to  the  logarithms  of  both  pressure  and  |suction| .  The  soil  volume  strain  equation 
which  shows  the  relation  between  matrix  suction,  mean  principal  stress,  and  osmotic 
suction  is  (4),  (Lytton  1977,  1994). 


-y)   =-Yh(zi)logio 


hF(Zi) 
h,(Zi). 


V"i 


-Ya(Zi)l08l 


°F(zi) 


-Yrt(Zi)log 


10 


(4) 


In  the  procedures  described  herein  the  assumption  will  be  made  that  the  osmotic 
component  has  no  change.  If  the  case  exists  that  the  osmotic  component  has  no  real 
change,  which  is  the  case  most  of  the  time,  then  that  term  within  the  volume  strain 
equation  can  be  neglected.  The  suction  compression  index  ,  yh(Zj),  and  the  mean 
principal  stress  compression  index,  y„(Zj),  will  be  discussed  in  more  detail  in  the  section 
titled  "The  Suction/Mean  Principal  Stress  Compression  Indexes".  The  initial  and  final 
suction  profiles,  hI(zi)  and  hF(Zj),  will  be  discussed  in  more  detail  in  the  section  titled 
"Suction  Profiles".  The  initial  mean  principal  stress,  a„  is  the  stress  level  below  which 
there  is  no  suppression  of  volume  strain  by  overburden  pressure.  This  is  the  stress  at 
which  the  curve  in  the  pressure-volume  plane  in  Fig  1 2  begins  to  depart  from  horizontal 
as  pressure  increases  from  zero  and  has  been  observed  to  be  approximately  40 
centimeters  of  soil  (Lytton  1994).  The  final  stress  at  depth  z;  >cxF(zi),  is  the  mean 
principal  stress  at  a  depth  z;  at  the  time  the  volume  strain  is  desired  and  is  given  by  (5). 


°F(zi)  = 


1  +  2-ko 


Y,-Zi 


(5) 


where  k0  is: 


k0  =  0.00  when  the  soil  is  badly  cracked. 

k0  =  0.33  when  the  soil  is  drying. 

k0  =  0.67  when  the  soil  is  wetting. 

k0  =  1 .00  when  the  cracks  are  closed  and  the  soil  is  swelling. 
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Once  the  volume  strain  for  an  incremental  layer  of  soil  is  known  the  vertical  strain  for 
that  layer  is  estimated  using  a  cracked  fabric  factor,  f; ,  as  shown  in  (6).  The  cracked 
fabric  factor  is  the  percent  soil  volume  strain  that  is  directed  vertically  in  a  soil  column. 


m 


(fi-ni^i  (»i 


McKeen  (1981)  back  calculated  vales  of  fj  and  found  that  fj  =  0.5  when  soils  are  drying 
and  fj=  0.8  when  soils  are  wetting.  The  total  vertical  differential  movement  in  a  soil 
mass,  Az ,  is  then  found  by  summing  the  product  of  vertical  strain  and  the  incremental 
depth  to  which  they  apply  as  shown  in  (7). 

Az  =  /J  — —     •  increment  (7) 

i=i  v  HA 

The  volume  strain  theory  previously  described  was  developed  by  Lytton  (1973,  1977, 
1994,  and  1995).  All  symbols  are  defined  in  the  List  of  Symbols. 

A  Simplified  Method  for  Identifying  Predominant  Clays 

There  are  three  areas  of  soil  characteristics  that  must  be  known  before  the 
volume  strain  theory  can  be  applied,  namely:  the  soil  suction  compression  index,  the 
soil  suction  profiles,  and  the  mean  principal  stresses  within  the  soil.  The  mean  principal 
stress  was  described  earlier  and  suction  profiles  will  be  discussed  in  a  later  section. 
However,  the  soil  suction  compression  index,  as  described  in  the  next  section,  depends 
on  proper  identification  of  the  predominant  clays  in  soils.  A  simplified  method  for 
identifying  these  predominant  clays  is  discussed  in  this  section. 

Before  Lytton  developed  the  methodology  to  predict  vertical  differential  soil 
movement,  Pearring  (1968)  and  Holt  (1969)  had  completed  a  correlation  chart  to  aid  in 
the  identification  of  predominant  clay  mineral  in  a  given  soil.  According  to  Pearring 
and  Holt's  research,  the  two  parameters  used  to  classify  the  clay  minerals  are  the  cation 
exchange  activity,  CEAc(Zj),  and  activity  ratio,  Ac(Zj).  Obtaining  the  measure  of  these 
two  parameters  requires  the  plasticity  index,  PI(Zj),  the  cation  exchange  capacity, 
CEC(Zj),  and  the  percent  of  clay  in  the  soil  passing  the  number  200  sieve,  Clay(Zj),  (8), 
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(9),  (10),  (11). 

PI(zi)  =  LL(z,)-PL(zi)  (8) 

O*,,-^  (9) 

^--^  (10) 

CEC(Zi)  nn 

cEAc(zi)=      ;■'  (11) 

Clay^) 
All  symbols  are  defined  in  the  List  of  Symbols.  Once  the  cation  exchange  activity  and 
the  activity  ratio  are  known  for  a  given  clay  mineral  the  type  of  clay  is  determined  from 
Fig  13. 

The  difficulty  comes  in  arriving  at  the  measure  of  cation  exchange  capacity. 
Until  Pearring  and  Holt's  method  of  identifying  predominant  clays  came  to  light  the 
geotechnical  engineer  had  to  rely  on  experience  or  the  more  reliable,  expensive,  and 
arcane  techniques  used  by  clay  mineralogist  and  soil  physicists.    However,  Pearring 
and  Holt's  method  still  required  engineers  to  find  the  cation  exchange  capacity  (CEC) 
of  the  soil,  which  is  a  soil  property  not  typically  evaluated  in  normal  soils  test 
laboratories.  Mojekwu  (1979)  realized  this  problem  and  developed  a  simplified  way  to 
arrive  at  the  CEC.  His  work  involved  data  collected  for  a  wide  range  and  number  of 
soil  samples  taken  in  the  state  of  Texas.  Among  other  data,  Mojekwu  collected  and 
analyzed  the  Atterberg  limits  of  these  soils  and  found  that  very  good  correlations  could 
be  made  that  linked  the  plastic  limit  and  the  liquid  limit  to  the  CEC  of  a  specific  soil. 
The  method  he  developed  requires  the  simple  Atterberg  limits  tests  to  find  the  plastic 
limit,  PL,  and  the  liquid  limit,  LL.  Mojekwu  ran  a  simple  regression  analysis,  the 
dependent  variable  being  CEC,  the  independent  variables  being  PL  and  LL  respectively, 
and  produced  very  high  R  squared  values  of  .9941  and  .9942  for  (12)  and  (13) 
respectively. 

CEC(z1)  =  PL(zi)"7  (12) 

CEC(zi)-LL(z1)09'2  (13) 
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The  larger  of  these  two  approximated  values  of  CEC  is  used  to  find  the  CEAc,  which  is 
then  used  along  with  Ac  to  determine  the  predominant  clay  minerals  in  Fig.  13.  These 
findings  by  Pearring,  Holt,  and  Mojekwu  paved  the  way  for  research  by  McKeen  in 
developing  a  simplified  method  to  identify  the  suction  compression  index,  which  is 
critical  to  predicting  differential  soil  movements. 
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Fig.  13.  (Pearring  1968  and  Holt  1969)  Chart  to  Identify  Predominant  Clays 


The  Suction/Mean  Principal  Stress  Compression  Indexes 

Suction  forces  on  a  soil  element  are  body  forces  rather  than  externally  applied 
forces  at  the  surface  of  the  element.  Because  the  suction  stresses  are  dominating  the 
effective  stresses  in  expansive  soils,  volume  change  studies  are  made  using  total  suction 
values.  It  is  the  slope  of  the  volume  change  versus  suction  curve  that  quantifies  soil 
response  to  moisture  changes.  Thus,  the  suction  compression  index,  yh(Zj),  is  defined  as 
the  slope  of  the  volume-total  suction  curve,  see  Fig  14  (Lytton  1977).  McKeen  (1981) 
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analyzed  soils  data  in  which  he  measured  volume  strain  corresponding  to  soil  suction 
values  to  arrive  at  a  corresponding  suction  compression  index  for  particular  types  of 
expansive  soils.  In  McKeen's  study,  he  measured  the  soil  volume  change  using  the 
conventional  oedometer  for  cases  when  the  soil  was  un-cracked  and  using  the 
coefficient  of  linear  extensibility  (COLE)  method  when  the  clay  was  cracked.  The 
measurements  of  soil  suction  were  taken  utilizing  a  variety  of  methods  as  well.  The 
primary  methods  were  the  thermocouple  psychrometer  and  filter  paper  methods.  Using 
the  basic  concept 


Mean 
Principal 
Stress 

Compression 
Index  (yc) 


Mean  Principal  Stress 
(Log  Scale) 


Suction 
Compression 
Index  (yh) 


Matrix  Suction  (Log  Scale) 


Fig.  14.  (Lytton  1977)  Pressure-Volume-Suction  Surface 


described  above,  McKeen  developed  the  chart  method  to  determine  yh  by  analyzing  a 
large  number  of  Soil  Survey  samples,  adopting  the  methodology  for  determining 
predominant  clays  by  Pearring  (1968),  Holt  (1969)  and  Mojekwu  (1979),  and  produced 
Fig  15  which  gives  values  of  yhl00  without  requiring  suction  tests.  The  value  of  yhl00 
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typically  ranges  from  .033  to  .220.  This  method  requires  only  the  activity  ratio  and  the 
cation  exchange  activity  of  the  soil  to  find  the  corresponding  yh  for  particular  soils.  The 
chart  method  is  used  in  the  procedures  to  calculate  yh  in  the  software  design  format 
procedures. 

Another  factor  to  consider  in  calculating  soil  volume  strain  is  the  effect  of 
applied  loads  in  suppressing  the  suction  volume  strain  by  overburden  pressure.  The 
mean  principal  stress  compression  index  is  required  to  complete  the  term  to  calculate 
the  amount  of  suppression  to  the  suction  volume  strain  by  overburden  pressure.  The 
mean  principal  stress  compression  index  is  calculated  by  (14)  (Lytton  1977). 
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Fig.  15.  (McKeen  1981)  Chart  to  Determine  the  Suction  Compression  Index 
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The  mean  principal  compression  index  is  also  calculated  by  the  commonly  used 
compression  index,  Cc,  by  (15)  (Lytton  1994). 

r-"iS:  (15) 

Typical  values  of  ya  were  calculated  from  results  of  swell  tests  and  swell  pressure  tests 
reported  by  Vigayvergiya  and  Ghazzaly  (1973)  according  to  (14).  The  range  of  these 
calculated  values  of  y0  fell  between  0.02  and  0.09  indicating  that  the  ratio  of  yh/yCT  is 
typically  at  or  slightly  below  1 .0.  Therefore  it  is  slightly  conservative  but  convenient  to 
assume  yh  is  equal  to  yCT  Until  there  is  research  that  develops  a  more  practical  and  cost 
effective  method  for  obtaining  more  accurate  values  of  yOJ  that  coefficient  will  assume 
the  value  of  yh  for  calculations  of  volume  strain  in  the  software  design  format 
procedures  described  in  this  document  (Vigayvergiya  and  Ghazzaly  1 973). 

Estimates  of  Unsaturated  Soil  Properties 

Before  a  clear  explanation  of  the  remaining  sections  is  possible,  estimates  of 
unsaturated  soil  properties  must  be  defined.  Specifically,  definitions  are  required  for  the 
following  soil  properties,  namely:  the  slope  of  the  suction  versus  volumetric  water 
content  wetting  line  shown  in  Fig.  10,  S(Zj),  Mitchell's  diffusion  coefficient,  a(Zj),  and 
Mitchell's  unsaturated  permeability,  ko(z;  )|h0|  orP(Zj).  The  slope  of  the  suction  versus 
volumetric  water  content  wetting  line  is  given  by  (16).  Mitchell's  diffusion  coefficient 
is  given  by  (17).  Finally,  Mitchell's  unsaturated  permeability  is  given  by  (18).  The 
process  flow  sheet  that  describes  procedures  for  calculating  all  soil  properties  required 
for  the  software  design  can  be  found  in  Appendix  A-88. 

S(Zi )  =  -20.29  +  0.1 555  •  LL^ )  -  0.1 1 7  ■  PI(zi )  +  0.0684  •  v(z, )  (16) 

a(Zi )  =  0.0029  -  0.000 1 62  •  S(z( )  -  0.0 1 22  •  y  „  (Zj )  (17) 

i      i  Otfz  )  •  Y    (z  ) 

K0(Zi)-h0  =  0.4343-    ,    ,J     d     '  (18) 

|S(Zj)|-yw 

The  above  equations  were  developed  by  Jayatilaka  et  al.  (1992)  and  Lytton  (1994).  All 
symbols  are  defined  in  the  list  of  symbols. 
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Edge  Moisture  Variation  for  Center  Lift  and  Edge  Lift  Cases 

As  discussed  earlier,  the  effects  of  vertical  differential  soil  movement  on 
commercial  and  residential  types  of  slabs-on-ground  is  dictated  by  the  perimeter 
moisture  effect  cases  and  the  total  amount  of  vertical  movement  between  two  distinct 
locations  separated  by  a  characteristic  distance.  The  characteristic  distance  is  described 
as  the  edge  moisture  variation  distance.  At  some  time  after  construction,  the  soil  under 
the  slab  will  reach  an  equilibrium  moisture  content.  However,  the  soil  at  the  edge  will 
vary  in  moisture  content  dependent  on  the  local  annual  weather  cycle.  Additionally,  the 
soil  underneath  the  slab  will  vary  in  moisture  content  up  to  a  distinct  distance,  the  edge 
moisture  variation  distance,  that  is  controlled  by  the  permeability  and  diffusivity  of  the 
soil.  This  differential  in  moisture  content  of  soil  underneath  the  slab  and  the  edge  of  the 
slab  is  what  causes  the  center  lift  and  edge  lift  distortion  modes.  As  part  of  research 
conducted  by  Jayatilaka  et  al.  (1992)  a  calibrated  finite  element  program  with  coupled 
transient  moisture  flow  and  elasticity  that  had  been  used  in  the  study  of  vertical 
moisture  barriers  provided  an  ideal  means  to  study  the  edge  moisture  variation  distance. 
The  unsaturated  soil  properties  discussed  in  the  previous  section  in  addition  to  the 
Thornthwaite  Moisture  Index  were  used  to  determine  the  relation  of  the  edge  moisture 
variation  distance. 

The  Thornthwaite  Moisture  Index  (TMI),  which  measures  the  location's 
potential  for  evapotranspiration,  is  defined  as  a  climatic  index  that  measures  the  water 
balance  of  an  area  to  determine  if  there  is  a  water  deficiency  or  water  excess. 

Both  edge  lift  and  center  lift  conditions  were  explored  using  several  hundred 
runs  with  the  program  (Jayatilaka  et  al.  1992).  Center  lift  conditions  were  simulated  by 
a  one  year  dry  spell  following  a  wet  suction  profile  condition.  Edge  lift  conditions 
were  simulated  by  a  one  year  wet  spell  following  a  dry  suction  profile  condition.  The 
edge  moisture  variation  distance  was  considered  to  be  that  distance  between  the  edge  of 
the  foundation  and  the  point  beneath  the  covered  area  where  the  suction  changed  no 
more  that  .2  pF  during  the  entire  period  of  simulation  (Lytton  1994). 
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The  dry  and  wet  conditions  used  annual  suction  variation  patterns  that  were 
appropriate  for  each  of  nine  different  climatic  zones  ranging  from  a  Thornthwaite 
Moisture  Index  of -46.5  to  +26.8,  spanning  the  range  found  in  Texas.  The  resulting 
edge  moisture  variation  distances  are  shown  in  Figs.  16  and  17.  No  distance  less  that 
2.0  feet  should  be  considered  for  design  purposes  (Lytton  1994).  In  the  procedures  for 
the  software  design,  the  edge  moisture  variation  distances,  eme(Zj)  and  e^Z;),  are 
calculated  using  Figs.  16  and  17. 

In  1 997  Lytton  developed  a  method  to  determine  the  edge  moisture  variation 
distance  for  a  particular  design  return  period  em  .  In  designing  slabs-on-ground 

recognition  must  be  taken  of  the  length  of  the  time  these  structures  must  be  in  service, 
and  of  the  severity  of  the  weather  patterns  that  may  occur  during  the  expected  life  of  the 
structure.  The  return  period  of  hydrologic  events  is  appropriate  to  use  in  estimating  the 
design  criteria  for  slabs-on-ground  on  unsaturated  soils.  Determining  the  edge  moisture 
variation  distance  for  a  selected  return  period  involves  applying  an  equation  that  relates 
the  edge  moisture  variation  values  for  the  10-year  and  50-year  return  periods  to  the 
design  edge  moisture  variation  value  for  the  selected  return  period,  (19). 


=  e      + 1  e      —  e      I 


rz   -z     ^ 


(19) 


Vz50  —  z10y 

The  edge  moisture  variation  distance  for  a  10  year  return  period  can  be  obtained  from 
Fig.  18.  These  edge  moisture  variation  distances  were  derived  by  back  calculation  from 
slabs  which  were  performing  successfully  in  San  Antonio,  Dallas,  and  Houston.  None 
of  the  slabs  were  more  than  10  years  old  at  the  time.  It  can  be  argued  that  the  design 
values  of  edge  moisture  variation  distance  represent  a  10  year  return  period,  em   .  The 
edge  moisture  variation  distance  charts  shown  in  Fig.  16  and  Fig.  17  represent  a  50  year 
return  period,  em   .  The  use  of  the  Gumbel  probability  density  function,  which  is 

commonly  used  to  represent  the  probability  of  weather  events,  may  be  used  to  establish 
the  risk  level  that  is  desired  for  design  in  accordance  with  the  expected  service  life. 
The  z,  scores  used  in  (19)  are  computed  from  the  Gumbel  cumulative  probability 
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distribution  curve,  (20).  The  return  period  is  represented  by  the  variable  r.  p  and  p  are 
P 


z,  = 


(20) 


In  1 


xJJ 


both  shape  factors  and  can  be  assumed  as  one.  Solving  for  an  edge  moisture  variation 
distance  for  a  particular  return  period  between  10  and  50  years  is  a  simple  two  step 
process:  (1)  simply  substitute  the  return  periods  of  10,  50  and  the  selected  return 
period  r  into  (20)  and  solve  for  z10,  z50  and  z,.  respectively ,  and  (2)  Substitute  values  of 
em    from  Fig.  18,  e      from  Fig.  16  or  Fig.  17,  and  values  of  z10,  z50,  and  z,.  previously 
calculated  into  equation  (19)  and  solve  for  the  em  .  A  common  design  period  for 
residential  and  light  commercial  structures  is  r  =  20  years  (5  percent  risk). 
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Fig.  16.  (Lytton  1994)  Edge  Moisture  Variation  Distance  for  Center  Lift  Conditions  (50  yrs) 
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Fig.  17.  (Lytton  1994)  Edge  Moisture  Variation  Distance  for  Edge  Lift  Conditions  (50  yrs) 
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Fig.  18.  (Wray  1978)  Edge  Moisture  Variation  Distances  for  Edge/Center  Lift  Conditions  (lOyrs) 
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Equilibrium  Suction 

Before  the  equilibrium  suction  profiles,  h^Zj),  wet  and  dry  limit  suction  profiles, 
h^^Zj)  and  h^Zj),  transient  state  suction  profiles,  h(Zi,t),  and  depth  to  constant  suction, 
Zn,  can  be  generated,  the  equilibrium  suction  ,  h,,,,  must  be  known  for  particular  soil 
profiles  and  surface  conditions.  Gay  (1994)  did  extensive  work  in  the  area  of 
calculating  the  mean  volumetric  moisture  content  for  a  given  soil  mass  dependent  upon 
the  soil's  depth  of  available  moisture,  d^,  and  the  location's  potential 
evapotranspiration.  The  depth  of  available  moisture,  d^,  is  defined  as  the  maximum 
depth  of  moisture  available  for  use  by  transpiring  vegetation,  which  is  stored  within  the 
soil  zone  down  to  the  depth  to  constant  suction.  The  Thornthwaite  Moisture  Index  is 
used  to  ascertain  if  a  climate  is  arid  or  humid  (Thornthwaite  1948).  Gay's  theory 
considers  dam  and  TMI.  With  these,  he  developed  a  set  of  functional  relationships  that 
are  used  to  calculate  the  mean  volumetric  moisture  content,  which  are  then  used  to 
calculate  the  equilibrium  suction  value  for  a  particular  soil  profile  and  location.    The 
approach  he  used  to  solve  for  these  equations  was  developed  by  Juarez-Badillo  (1975) 
in  which  functional  domains  for  the  problem  are  first  established  as  shown  in  Fig.  1 9. 
Function  (21)  was  developed  to  satisfy  the  boundary  conditions  in  Fig.  19.  Then 
through  an  assumption  of  linear  proportionality  of  the  rates  of  change  of  the  two 
complete  functions,  Gay  obtained,  (22)  and  (23). 
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Fig.  19.  (Gay  1994)  Domains  of  Functions  f(d)  and  g(T) 


Solving  (23)  for  d  gives  (24). 

d.„ 


d  = 


d     -d, 

am  1 


1  + 
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(24) 


All  symbols  are  defined  in  the  list  of  symbols.  Gay  then  determined  how  the 
coefficients  y,  d„  and  T,  depend  upon  d^  using  the  pattern  search  technique  and  arrived 
at  the  expressions  (25),  (26),  and  (27). 


Y  =  0.039337dam+ 1.357033 


d1  =  0.449079dam  + 0.304560 
T,  =  0.06265  ldam  +  59.53593 


(25) 
(26) 
(27) 


These  relationships  facilitate  the  calculation  of  mean  moisture  depths,  dm,  for  all  values 
of  TMI  for  any  value  of  available  moisture  depth  between  10cm  and  50cm,  as  illustrated 
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in  Fig.  20  (Gay  1994).  The  development  of  (24)  is  important  because  it  is  used  to  solve 
for  the  mean  volumetric  water  content,  0m   The  value  0m  is  applied  as  a  target  to  solve 
for  \\n  in  an  iterative  process  using  Nieber's  (1981)  equation  (34).  Nieber's  equation 
relates  suction  to  volumetric  water  content.  The  method  to  apply  this  iterative  process 
to  solve  for  the  equilibrium  suction  value  can  be  found  in  Appendix  A-90. 
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Fig.  20.  (Gay  1994)  Theoretical  Relationships  for  Mean  Annual  Moisture  Depth  vs  TMI 

The  equation  used  to  solve  for  the  mean  annual  moisture  depth,  (24),  developed 
by  Gay  is  rearranged  to  arrive  at  the  mean  volumetric  water  content,  (32).  Below  is  the 
step  by  step  process  to  obtain  the  equation  to  solve  for  9m.    First  note  the  following 
equations  from  Gay  (1994),  (28),  (29),  (30),  (31).  Fig  21  provides  a  graphical 
explanation  for  some  of  the  symbols  in  these  equations.  All  symbols  are  defined  in  the 
list  of  symbols. 
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Dividing  through  equation  (24)  by  zc  gives  (32). 
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Solving  for  d^/z,.  in  (29)  and  substituting  in  (32)  gives  (33). 


33 


(29) 
(30) 
(31) 


(32) 


9w.-Qdry 
Um  "  f  \   +0dry 


1  + 


d     -d, 

am  1 


(33) 


W  J 

The  steps  to  arrive  at  equations  which  solve  for  0wl,  0dry,  and  0m  using  Nieber's 
equation  (34)  are  given  below. 


h  J  = 


0-0 


0-0. 


L-i 


(34) 


Introduce  dummy  variable  r,  (35). 


r  = 


es-e 

0-0, 


(35) 
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Solving  (34)  for  G  gives  (36). 

e„+r-e. 

G  = 


1  +  r 


(36) 


Substituting 


I         |B        I  |B 


,  or 


dry 


back  into  (36)  for  r  gives  (37),  (38),  or  (39) 


A    '     A  A 

depending  on  which  volumetric  moisture  content  is  to  be  found. 


MOISTURE  CHANGES  IN    UNSATURATED   SOILS 

pF  5.7      pF  4.2  pF  1 .6-3 .2       pF  1 .0- 1 .5 


Matrix 

Potential 

pF 


#air     #root  #fc     #s 

.  Moisture  -  Loss  Roots     i 


Depth 


Q  Volumetric 
Moisture 
Content 


dfc    =  depth  of  moisture   at  field   capacity 
dam  =  depth  of  available   moisture 
Zr    =  rooting   depth  in  soil 


Fig.  21.  (Lytton  1994)  Moisture  Changes  in  Unsaturated  Soils  Between  Typical  Wet  and  Dry  States 
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e  = 


I         |B 

|hm| 

(37) 


1  + 


B 


A 

I      lB 

e.+^M-e. 


e..,  = 


i  + 


K 


B 


(38) 


A 

|B 


hdry 


'+     A 

The  equilibrium  suction  value,  h,,,,  can  then  be  solved  through  an  iterative  process  using 
equations  (33),  (37),  (38)  and  (39)  as  described  in  Appendix  A-90.  The  equilibrium 
suction  profile,  hjz{)  for  any  pervious  soil  profile  is  defined  by  decreasing  the 
equilibrium  suction  value  at  the  surface,  h^  a  centimeter  in  suction  for  every  centimeter 
increase  in  depth  (McQueen  and  Miller  1968).  The  process  for  determining  the 
equilibrium  suction  profile  is  detailed  in  Appendix  A-96. 

Suction  Profiles 

It  is  important  to  note  that  moisture  changes  within  a  natural  soil  are  determined 
by  the  ratio  of  the  period  of  evaporation  to  the  period  of  precipitation.  A  measure  of 
this  period  of  evaporation  to  the  period  of  precipitation  is  determined  through  a  method 
developed  by  Thornthwaite  (1948).  In  the  previous  section  the  equilibrium  suction  is 
determined  based  on  relations  to  the  climatic  water  balance  expressed  in  terms  of  the 
Thornthwaite  Moisture  Index  (1948).  When  a  marked  separation  occurs  between  wet 
and  dry  seasons,  a  large  seasonal  variation  in  soil  moisture  content  occurs,  whereas  in 
areas  which  are  either  predominantly  dry  or  predominantly  wet  for  most  part  of  the  year, 
the  changes  in  soil  moisture  content  is  not  so  marked.  A  location  where  predominant 
clays  are  abundant  in  which  the  climate  is  marked  by  this  large  separation  in  the  periods 
of  evaporation  and  periods  of  precipitation  will  experience  problems  with  differential 
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movement  to  the  greatest  degree.  The  state  of  moisture  is  dependent  on  the  state  of 
suction  and  tends  towards  an  equilibrium  suction  profile,  as  defined  in  the  previous 
section,  with  a  moisture  source  either  at  the  boundaries  or  within  the  soil  mass. 
When  the  state  of  suction  is  measured  at  intervals  of  depth  in  a  soil  profile,  the  resulting 
relation  between  suction  and  depth  is  termed  the  suction  profile.  Mitchell  measured 
actual  suction  profiles  in  the  area  of  Adelaide,  South  Australia  and  found  a  very 
consistent  and  typical  trend  that  could  be  modeled.  In  his  study  he  developed  a  method 
based  on  constitutive  soil  properties  to  model  the  suction  profile  within  a  soil  profile 
base  on  the  local  annual  weather  cycle,  n,  soil  diffusion,  a(Zj),  equilibrium  suction 
value,  hm,  and  the  amplitude  suction  value,  U0  The  maximum  limit  suction  profile, 
ri^Zj),  and  minimum  limit  suction  profile,  h^fe)  begin  at  the  limiting  field  condition 
suction  value  at  the  surface,  and  decrease  at  an  exponential  rate  toward  the  equilibrium 
profile,  hmCZj).  These  limit  suction  profiles  will  approach  the  equilibrium  suction  profile 
at  the  depth  to  constant  suction,  z^,  which  will  be  defined  in  the  next  section.  The  limit 
suction  value  depends  on  the  type  of  surface  conditions  and  is  defined  in  the  section 
titled  "Typical  Suction  Levels  in  the  Field".  The  dry  limit  suction  value  is  usually  4.5 
pF  for  soil  with  vegetation  at  the  surface  and  6.0  pF  for  bare  soil  at  the  surface.  Refer  to 
Fig.  22  for  a  typical  maximum  limit  suction  profile.  The  wet  limit  suction  value,  which 
is  somewhat  drier  than  the  field  capacity  suction  value,  is  usually  2.5  pF.  Refer  to  Fig. 
23  for  a  typical  minimum  limit  suction  profile.  Mitchell  uses  an  idealized  model  of  the 
annual  weather  surface  cycle  which  he  represents  by  (40).  This  is  a  good  representation 
for  design  because  it  assumes  that  soil  at  the  surface  approaches  the  maximum  and 
minimum  moisture  conditions  at  two  distinct  and  separate  times  of  the  year.  As  stated 
earlier,  this  separation  between  wet  and  dry  periods  produces  worst  case  differential  soil 
movement  scenarios.  Mitchell  combines  the  soil  suction  relationship  with  the 
relationship  that  describes  the  local  surface  weather  cycle  and  arrives  at  equation  (42). 

Mitchell  found  that  the  suction  change  due  to  the  effects  of  climate,  drainage  and 
site  cover  is  a  periodic  function  of  time.  The  suction  at  any  time  and  at  any  depth  in  the 
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soil  profile  is  determined  by  solving  the  diffusion  equation,  (4 1 ),  for  a  particular 
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Fig.  22.  Typical  Dry  Limit  Suction  Profile  for  Bare  Soil  at  the  Surface 
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Fig.  23.  Typical  Wet  Limit  Suction  Profile  for  Bare  Soil  at  the  Surface 
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set  of  boundary  conditions.  The  solution  to  arrive  at  the  soil  suction  profile  equation 
(42)  can  be  found  in  Mitchell's  1981  report  (values  of  h(Zj,t),  h^z^t),  and  U0  are  all  in 
units  of  pF). 

h(Zj  =  0,t)  =  cos 


27t-n-t-(Zj  =0)-  I     "'"  (40) 


^  V1       '\a(Zi=0). 

ah         d2h 

-^r  =  a^-y  (41) 

5t         dz2  v     ' 


,    n-7i 


h(zi,t)  =  hm(zi)  +  U0e      ,ow.cos 


a(Zi)J 


(42) 


Based  on  the  surface  and  depth  boundary  conditions,  referred  to  as  moisture  effect 
cases,  (42)  is  applied  to  produce  suction  profiles  for  many  surface  conditions  and  soil 
profiles.  Appendix  A  contains  detailed  algorithms  that  describe  precisely  how  to  apply 
(42)  to  arrive  at  suction  profiles  resulting  from  combinations  of  different  moisture 
effect  cases  and  soil  profiles.  Later  sections  in  "The  Procedures,  in  a  Software  Design 
Format"  chapter  will  provide  a  more  detailed  discussion  with  respect  to  application  of 
(42)  in  producing  suction  profiles. 

Depth  to  Constant  Suction 

The  depth  to  constant  suction  is  the  depth  in  a  soil  profile  to  which  there  is  no 
longer  a  significant  seasonal  change  in  soil  suction.  This  depth  has  been  reported  to  be 
anywhere  from  6  feet  to  as  much  as  26  feet  dependent  upon  the  type  of  soil,  seasonal 
weather  trends,  and  the  type  of  surface  conditions  such  as  vegetation  and  slope.  Mohan 
and  Rao  reported  where  the  active  zone  extended  to  a  depth  of  1 3  to  16  feet,  in  the 
black  cotton  soils  of  India.  Kassiff  reported  movement  of  an  exposed  clay  at  depths  up 
to  20  feet  in  Israel,  with  the  largest  differential  movements  being  measured  within  an 
active  zone  of  6  Vz  feet,  where  moisture  changes  are  maximal  (Mitchell  1980).  As  one 
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uses  the  procedures  contained  herein  to  predict  the  differential  movement  in  soils,  it 

becomes  clear  that  most  of  the  movement  does  in  fact  occur  nearer  the  surface.  This 

reduction  in  the  amount  of  differential  soil  movement  with  an  increase  in  depth  is 

attributable  to  the  combination  of  small  seasonal  changes  in  soil  moisture  content  with 

an  increase  in  depth,  and  the  increase  in  restriction  to  volume  strain  due  to  overburden 

stresses  with  an  increase  in  depth.  The  algorithm  to  solve  for  the  depth  to  constant 

suction  is  based  primarily  on  Mitchell's  equation  (42)  that  is  described  in  the  previous 

section. 

Lytton  and  the  author  are  responsible  for  developing  the  following  steps  to 

arrive  at  the  equation  to  calculate  the  depth  to  constant  suction.  Only  the  limit  profile 

aspect  of  Mitchell's  equation,  (42),  will  be  used.  Setting  the  cosine  term  equal  to  one 

forces  this  equation  to  generate  only  the  maximum  and  or  minimum  suction  envelopes. 

Refer  to  (43). 

( 
*— -~— : — —        \  n'n 

(43) 


h(zi,t)  =  hm(zi)  +  U0e 


-(z,) 


1.0 


Mitchell's  equation  in  the  form  of  (43)  is  used  in  most  of  the  algorithms  to  calculate  the 
volume  strain  for  post  equilibrium  moisture  effect  cases.  Equation  (43)  is  rearranged  to 
isolate  the  exponent  term,  (44). 


h(zl,t)-hm(zi)  =  U0e 


"(Zi) 


a(Z|) 


(44) 


(44)  is  then  divided  through  by  U0,  (45). 


h(zi,t)-hm(zi)  -ta>>a 

U„ 


(45) 
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The  last  step  is  to  take  the  natural  log  of  both  sides  and  solve  for  Zj,  (46). 

1 


zm  =  -In 

m 


^h(Zi,t)-hm(Zi)^ 

u0 


n-7i 


(46) 


Va(Zj) 

The  depth  to  constant  suction,  z^,  has  been  obtained  when  the  amplitude  at  depth  z{  is 
equal  to  0.2  pF.  Later  sections  in  "The  Procedures,  in  a  Software  Design  Format" 
chapter  will  provide  a  more  detailed  discussion  with  respect  to  application  of  (45)  and 
(46)  in  calculating  the  depth  to  constant  suction. 

Horizontal  Velocity  in  Expansive  Soils 

The  ability  to  calculate  horizontal  velocity  of  water  flow  in  soils  is  critical  for 
the  application  of  the  soil  volume  strain  theory.  Predicting  vertical  differential  soil 
movement  when  a  horizontal  barrier  is  employed,  or  when  a  moisture  effect  case  is 
located  a  distinct  distance  from  the  slab's  edge  depends  upon  the  soil's  horizontal  water 
flow  characteristics.  The  horizontal  velocity  is  calculated  and  then  used  to  calculate  a 
suction  profile  a  distinct  distance  from  another  known  suction  profile.  Lytton  and  the 
author  developed  the  horizontal  velocity  equation  by  applying  Darcy's  law,  (47),  to  the 
unsaturated  permeability  relationship  for  clay  soils  developed  by  Mitchell,  (18),  to 
arrive  at  the  horizontal  velocity  equation  for  unsaturated  soils,  (52).  The  steps  to  derive 
(52)  are  described  below.  First  Darcy's  law,  (47),  is  solved  for  dh,  (48). 


,'drA 

•ij  (47) 


Vx(z.) 
dh  =  -      ,    '     dx  (48) 


Darcy's  permeability  constant  is  set  equal  to  Mitchell's  unsaturated  permeability,  (49). 
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k  =  ^^  (49) 

h 

Substitute  (49)  in  for  Darcy's  permeability  in  (48)  and  separate  variables  dependent  on 

distance  and  variables  dependent  on  suction  to  arrive  at  (50). 

V.(Zi)       .       dh 

Integrate  the  left  side  of  (50)  from  zero  to  a  distance  x,  and  the  right  side  of  (50)  from  an 
entry  suction  value,  hentry(Zj),  to  an  exit  suction  value,  hexit(Zj).  Refer  to  equation  (51). 

Vx(z=)     xf         'wrZi)dh 


K0(Zi) 


)      f  r     an 

ir/dx=.  {  t  (51) 


0    0  h0,lly(Zi) 

Solve  equation  (51)  for  Vx(Zj)  to  arrive  at  equation  (52). 


V   /      .        Ko(Zi)ho     , 

vx(zi)  =  — : — ln 


'W^ 


(52) 


KM)  J 

The  entry  suction  value  is  the  suction  value  at  x  =  0.  The  exit  suction  value  is  the 
suction  value  at  x.  Once  the  horizontal  velocity  is  known  between  two  distinct  locations 
in  an  incremental  soil  layer,  the  suction  value  for  any  point  between  those  two  distinct 
locations  is  solved  for  using  equation  (53).  Equation  (53)  is  derived  by  solving  for  the 
exit  suction  in  (52). 


VX(Z;)-X 


hexit(zi)  =  hentry(zi)-e^(z')h°  (53) 


Post  Construction  Theory 

The  worst  case  differential  soil  movement  may  occur  in  the  first  dry  season 
subsequent  to  the  construction  of  the  slab  during  the  wettest  time  of  the  year,  or 
during  the  first  wet  season  subsequent  to  construction  of  the  slab  during  the  driest  time 
of  the  year.  This  type  of  moisture  effect  case  is  referred  to  as  the  post  construction 
case.  The  post  construction  case  is  a  transient  case  for  calculating  volume  strain  of 
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expansive  soils  immediately  after  construction  and  before  the  soil  under  the  center  of 
the  slab  has  reached  an  equilibrium  moisture  content.  For  example,  when  a  slab  is 
constructed  the  soil  is  "capped"  and  the  moisture  content  changes  slowly  toward  its 
equilibrium  value  while  the  perimeter  soil  shrinks  and  swells  with  the  change  in 
seasonal  moisture.  This  type  of  scenario  provides  conditions  for  the  post  construction 
moisture  effect  case.  The  moisture  content  of  the  "capped"  soil  looses  or  gains 
moisture  for  some  period  of  time  until  the  soil  under  the  center  of  the  slab  reaches  an 
equilibrium  moisture  content.  The  rate  at  which  the  soil  approaches  equilibrium  is 
dictated  by  the  soil's  diffusivity  and  unsaturated  permeability.  Once  the  soil  under  the 
slab  reaches  the  equilibrium  moisture  content,  subsequent  suction  distributions  and 
differential  movements  must  be  determined  using  the  post  equilibrium  methodology. 

The  "Damped"  Slab  Suction  Profile 

Lytton  and  the  author  developed  the  following  set  of  equations  to  model  the  soil 
under  the  slab  subsequent  to  construction  in  a  transient  case  in  which  a  soil  suction 
profile  approaches  an  equilibrium  suction  profile.  The  idea  is  to  calculate  the  soil 
suction  profile  under  the  slab  at  the  time  of  construction,  hsiab(Zj ,tconst/dry ) ,  and  then 
dampen  this  profile  until  it  approaches  the  equilibrium  profile.  The  soil  suction  profiles 
under  the  slab  and  at  the  slab  edge  are  solved,  with  respect  to  time,  based  on  some 
knowledge  of  the  number  of  annual  weather  cycles.  The  annual  weather  cycle  will  be 
assumed  to  vary  according  to  the  function  described  by  (40).  To  demonstrate  this 
method,  the  annual  weather  cycle  is  assumed  to  be  one,  n  =  1 .  This  means  that,  in  one 
year  there  will  be  one  wet  season  and  one  dry  season.  The  driest,  or  wettest,  month  of 
the  year  must  be  known.     In  this  case,  June  is  selected  as  the  driest  month  of  the  year. 
Fig.  24  is  an  example  of  a  weather  cycle  modeled  using  (40)  for  bare  soil  at  the  surface 
(n  =  1,  and  June  is  the  driest  month  of  the  year). 

The  approach  used  to  derive  the  equations  to  damp  the  suction  profile  at  time  of 
construction  back  to  the  equilibrium  suction  profile  was  developed  by  Juarez-Badillo 
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(1975)  in  which  functional  domains  for  the  problem  are  first  established  as  shown  in 
Fig.  25.  Function  (54)  was  developed  to  satisfy  the  boundary  conditions  in  Fig.  25. 

For  Bare  Soil  At  The  Surface 
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Fig.  24.  A  Weather  Cycle  Modeled  by  a  Cosine  Function,  n  =  1,  June  is  the  Driest  Month 


f(h)  = 


hm  —  h     hm  —  h c. 

m  m  const 


(Case  1 ,  when  h,^  is  greater  than  hconst.)  (54) 


f(h)  = 


hm  -  h         hm  -  hconst 


(Case  2,  when  ^  is  less  than  hconst.) 


(54) 


The  suction  value  with  respect  to  time  and  depth  for  a  profile  that  is  modeled  to  tend 
toward  an  equilibrium  suction  was  solved  for  cases  1  and  2.  Both  cases  arrived  at  the 
same  expression  (59).  The  steps  to  arrive  at  (59)  are  described  below.  Through  an 
assumption  of  linear  proportionality  of  the  rates  of  change  of  the  function  equations 
(55)  is  established,  and  the  constant  of  proportionality,  (3,  is  equal  to  1. 


df (h)  _    dt 
f(h)  "  P  t 


(55) 
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Fig.  25.  Domains  of  Functions  h,  t,  and  1(h) 


The  left  side  of  equation  (55)  is  integrated  from  an  arbitrarily  chosen  suction  value,  h,  to 

a  suction  value  halfway  between  h,,,  and  hconst,  h.  The  right  side  of  equation  (55)  is 

integrated  from  a  time  corresponding  to  h,,  t,,  to  the  time  required  to  approach  h,  t.  This 

integration  is  given  by  equation  (56). 

1  1 

t 


hm-h 
hm-h, 


hm-hc 

1 


=  (Mn 


(56) 


nm        n  const 

Equation  (56)  is  rearranged  to  collect  all  constant  variables  on  one  side  of  the  equation 
and  set  equal  to  a  variable,  r.  Refer  to  equation  (57). 


/ 


h-h 


hi-hm 
h-h, 


XiJ 


^p 


=  r 


hm-h 
Now  equation  (57)  is  rearranged  to  solve  for  the  suction  value,  h,  given  by  (58). 

h=  hconst  +r'hm 


(57) 


1  +  r 


(58) 
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Substitute  (57)  back  in  to  (58)  to  arrive  at  the  equation  used  to  solve  for  the  suction 
profile  under  the  slab  some  time,  t,  after  the  slab  was  "capped"  as  a  result  of  slab 
construction,  (59). 

h,  -h    ^ 


hcons,(Z>  toons,  )  +  hm(Z)L 


t,-P-tP 


= K     I  s   '" — ' 

h,  -h. 


1  + 


_°l  I  .t  -P  .tP 


^hm  -h. 


The  Time  Equation 

Before  equation  (59)  can  be  applied,  a  relationship  must  be  derived  to  determine 
the  amount  of  time  required  for  water  to  move  upward  or  downward  to  change  the 
suction  at  any  given  depth  halfway  between  its  initial  value  and  the  equilibrium  value 
of  suction.  This  time  equation  for  t,  is  given  by  (61).  The  derivation  of  (61)  is  as 
follows.  Darcy's  law  of  permeability  is  the  relationship  used  to  introduce  the  time 
variable.  Velocity  is  defined  as  a  distance  through  which  water  flows  in  soil,  x,  traveled 
in  some  amount  of  time,  t,.  Velocity  described  in  terms  of  permeability  of  water  in 
expansive  soils  in  which  the  time,  t„  required  to  obtain  half  of  the  change  in  suction  (the 
time  to  obtain  half  of  the  change  in  suction  was  arbitrarily  chosen),  0.5  •  Ah ,  is  given  by 
(60) 

db^         T  0.5  Ah 


jf  x  an--~  u.d  •  An 

^=HFP(z>>-hH  (60) 

Equation  (60)  is  solved  for  t„  (61)  In  equation  (60),  the  variable  x  is  replaced  with 
equation  (44),  where  the  ratio  of  suction  values  is  set  equal  to  one  half. 

-J-^Mn(0.5) 

I    V  n-7r  J 

t,  = K  (z)  h (61) 

°-5'     °04343V(h(zi'tcons,)-hm(z)) 

The  time  variable  represented  by  (61)  is  used  in  (59)  to  solve  for  the  suction  profile 
some  time,  t,  after  soil  is  "capped"  as  a  result  of  slab  construction. 
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THE  PROCEDURES,  IN  A  SOFTWARE  DESIGN  FORMAT 

Introduction  to  the  Software  Design  Approach 

The  tools  described  in  "The  Theory  of  Differential  Soil  Movement"  chapter  are 
all  that  is  needed  to  develop  the  equations  and  terms  used  in  the  procedures  described  in 
this  chapter.  Appendix  B  is  a  graphic  representation  suggested  as  a  flow  of  vertical 
differential  soil  movement  algorithms  in  a  windows  environment  type  software  design 
layout.  The  suggested  software  design  layout,  Appendix  B,  is  modeled  after  the 
PTISLAB  WIN  1.0  software  package  owned  by  Geostructural  Tool  Kit,  Inc.  (1996). 
The  focus  of  this  document  is  not  software  design  or  marketing,  those  skills  are  left  to 
experts  in  that  area  of  work.  The  focus  is  to  present  the  equations  and  terms  in  easy  to 
follow  algorithms,  as  described  in  Appendix  A,  so  that  engineers  or  programmers  can 
apply  this  soil  volume  strain  technology  in  everyday  practice. 

A  question  that  a  practicing  engineer  may  ask  is,  "What  practical  engineering 
cases  can  these  procedures  address?"  The  answer  is,  cases  involving  a  surface  drained 
soil  profile  in  which  the  slope  is  stable  and  surface  conditions  are  those  described 
below: 

1 .  Bare  soil  at  the  surface. 

2.  Grass  at  the  surface  with  shallow  roots. 

3.  A  flower  bed  at  the  surface  with  a  known  depth  of  flower  bed  zone,  dzone. 

4.  A  tree  with  a  known  depth  of  root  zone,  dzone. 

The  procedures  consider  design  effects  that  are  introduced  to  change  the  magnitude  of 
the  vertical  differential  soil  movement.  These  design  effects  include:  vertical  barriers, 
horizontal  barriers,  and  the  effects  of  locating  the  surface  moisture  effect  a  distinct 
distance  from  the  foundation,  dsource.  All  of  these  design  effects  can  be  rigorously 
analyzed  to  enable  an  engineer  the  ability  to  choose  the  effect  which  will  provide  the 
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desired  design  results. 

The  algorithms  contained  herein  address  all  cases  and  design  effects  listed  above 
by  way  of  two  design  approaches:  (1)  the  post  construction  approach,  and  (2)  the  post 
equilibrium  approach.  An  engineer  would  apply  the  algorithms,  used  in  both 
approaches,  to  address  all  possible  combinations  of  cases  and  design  effects  presented. 
The  engineer  compares  the  results  of  these  cases  to  find  scenarios  that  produce  vertical 
differential  soil  movement  that  generates  maximum  bending,  shear  and  deflection.  This 
information  is  used  to  generate  adequate  foundation  designs.  The  most  promising 
scenario  is  to  use  the  power  of  today's  computing  technology  to  perform  the  iterative 
calculations  and  tests  required  to  arrive  at  these  design  quantities.  A  programmer  could 
apply  software  technology,  such  as  PTISLAB  WIN  1.0,  to  model  a  slab's  interaction 
with  the  differential  soil  movements  produced  by  the  algorithms  contained  herein  to 
arrive  at  these  maximum  design  cases. 

The  algorithms  are  presented  in  a  format  that  can  be  easily  programmed  using 
simple  spreadsheet  type  software.  First,  input  required  for  all  calculations  are  entered. 
Soil  properties  and  other  related  soil  characteristics  are  calculated  and  stored  for  use  in 
algorithms  presented  later.  Equilibrium  suction  values,  h,,,,  and  depth  to  constant 
suction  values,  z^,  are  calculated  for  a  specific  location  and  specific  sets  of  surface 
conditions.  With  the  above  complete,  all  moisture  effect  cases  and  design  effects  listed 
previously  are  addressed  using  the  post  construction  solutions  and  the  post  equilibrium 
solutions. 

Input 

When  considering  all  cases  and  all  calculations  required  to  complete  algorithms, 
an  accounting  of  the  variables  required  for  input  becomes  important.  Great  effort  is 
made  to  keep  the  complex  nature  of  data  required  for  input  down  to  a  minimum.  The 
goal,  relative  to  input  of  data,  is  to  ensure  that  all  variables  required  for  input  are  items 
that  are  easy  and  economical  for  practicing  engineers  to  obtain.  This  is  the  same  goal 
authors  had  in  mind  while  developing  the  simplified  method  for  identifying 
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predominant  clay  minerals  (Pearling  1968,  Holt  1969  and  Mojekwu  1979),  the  methods 
for  obtaining  the  suction  compression  index  and  the  mean  principal  stress  compression 
index  (McKeen  1981  and  Viayvergiya  and  Ghazzaly  1973),  and  the  method  for 
obtaining  the  edge  moisture  variation  distance  (Jayatilaka  et  al.  1992).    Below  is  a  full 
list  of  items  required  for  input  to  perform  every  calculation  in  the  algorithms  described 
in  Appendix  A. 

A,B  Gardner's  coefficients  (Appendix  A-90  and  Table  2) 

B  Number  of  boring  samples  (Appendix  A-84) 

Dc  Depth  to  the  cA  suction  measurement  in  a  sample  (Appendix  A-87) 

Ls  Depth  to  the  bottom  of  layer  s  (Appendix  A-84) 

LLS  Liquid  limit  for  layer  s,  in  percent  (Appendix  A-84) 

Mb  Number  of  layers  in  the  b*  boring  sample  (Appendix  A-84) 

N  Number  of  suction  measurements  in  a  sample  (Appendix  A-87) 

PLS  Plastic  limit  for  layer  s,  in  percent  (Appendix  A-84) 

TMI  Local  annual  Thornthwaite  Moisture  Index,  (Fig.  26  and  Appendix  A-83) 

Xb  X  coordinate  for  the  0th  boring  sample  (Appendix  A-84) 

Yb  Y  coordinate  for  the  bth  boring  sample  (Appendix  A-84) 

d^  Available  moisture  depth  of  the  soil  mass  (Appendix  A-90) 

dhbar  Size  of  horizontal  barrier  (Appendix  A-86) 

^source  Distance  effect  is  away  from  the  edge  of  barrier/slab  (Appendix  A-86) 

cL^  Depth  of  vertical  barrier  (Appendix  A-86) 

dzone  Depth  of  the  root  or  flower  bed  zone  (Appendix  A-86) 

hc  cth  suction  measurement  in  a  sample  (Appendix  A-87) 

h^  Dry  suction  limit  for  soils  in  the  field,  bare  soil  at  the  surface  usually 
has  a  value  of  6.0pF  and  soil  with  vegetation  at  the  surface  usually 
has  a  value  of  4.5pF  (It  may  be  requested  in  various  algorithms. 
Appendix  A-90  is  the  first  time  h^  is  requested  as  input). 
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h^  Wet  limit  suction  for  soils  in  the  field,  usually  2.5pF  (It  may  be  requested 

in  various  algorithms.  Appendix  A-90  is  the  first  time  hw,  is  used), 
n  Local  annual  weather  cycle  (Appendix  A-83). 

tconst  The  month  the  slab  was  constructed  (Appendix  A-l  07) 

tdry  Driest  month  of  the  year  (Appendix  A- 1 07) 

t^,  The  month  the  soil  volume  change  is  sought  (Appendix  A- 1 07) 

fj.s  Percent  passing  2  micron  for  layer  s  (Appendix  A-84) 

vs  Percent  passing  number  200  sieve  for  layer  s  (Appendix  A-84) 

More  detailed  explanations  of  how  these  variables  are  used  is  described  in  the  following 
sections  of  this  document.  Appendix  A  gives  a  complete  process  flow  detailing  the  use 
of  the  variables  listed  above  and  the  variables  listed  in  the  list  of  symbols. 

Input,  as  requested  by  the  process  flow  in  Appendix  A  and  as  graphically 
represented  in  the  software  design  layout  sheets  in  Appendix  B,  is  broken  down  into  six 
major  categories: 

1.  Project  Information  input.  This  includes  project  name,  date  the  project  started, 
name  of  the  engineer  assigned  to  the  project,  and  the  site  location.  This  information 
is  basic  information  an  engineering  firm  may  want  to  establish  in  a  database  so  that 
information  can  be  readily  accessed.  The  only  information  critical  to  calculating 
differential  soil  movement  is  the  site  location.  The  site  location  can  be  used  as  the 
means  to  locate  the  TMI  for  a  site  without  requiring  the  user  to  manually  research 
this  bit  of  information.  Refer  to  Appendix  A-82  for  a  process  flow  and  Appendix  B- 
1 17  for  a  program  design  layout.  The  TMI  for  a  particular  site  is  found  in  Fig.26. 

2.  Foundation  Layout.  This  is  basic  surveying  information  to  establish  a  three 
dimensional  grid  for  logging  information  related  to  appropriate  physical  locations. 
Refer  to  Appendix  A-83  for  process  flow  and  Appendix  B-l  18  for  program  design 
layout. 

3.  Soil  Layer  Property  Input.  These  values  include  boring  sample  data  and  soil  layer 
property  data.  The  data  entered  here  will  provide  enough  information  to  generate  a 
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three  dimensional  picture  of  the  soil  profile.  Refer  to  Appendix  A-84  for  a  process 
flow  and  Appendix  B- 1 1 9,  1 20,  1 26  for  program  design  layouts.  The  soil  profile 
used  in  all  example  calculations  is  shown  in  Table  1 . 


Houston       L^  30 
o 


Fig.  26.  (Wray  1978)  Thornthwaite  Index  for  Texas  (20  year  average  1955-1974) 


Table  1.  Soil  Profile  Considered  in  all  Sample  Calculations 


Soil 

Ls 

LLS 

PLS 

vs 

Ms 

Yd 

Layer  No. 

ft 

% 

% 

% 

% 

lbs/ft3 

s  =  1 

5 

56.0% 

22.0% 

65.0% 

25.0% 

120 

s  =  2 

9 

67.0% 

24.0% 

76.0% 

46.0% 

120 

s  =  3 

13 

44.0% 

18.0% 

47.0% 

18.0% 

120 

s  =  4 

30 

53.0% 

20.0% 

64.0% 

31.0% 

120 
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4.  Moisture  Effect  Input.  This  gives  the  user  the  option  to  input  the  type  and  location 
of  the  moisture  effect  cases  considered  in  vertical  differential  soil  movement 
calculations.  Refer  to  the  plan  view,  Appendix  C-132,  for  a  view  of  the  foundation 
used  as  the  example  throughout  this  document.  The  input  for  location  and  moisture 
effect  case  type  can  be  found  listed  on  the  program  design  layout  Appendix  B-121 . 
For  example,  case  "1-5"  describes  the  sidewalk  with  grass  covered  lawn  on  the 
north  and  west  portions  of  the  building  shown  in  Appendix  C-132.  The  number  "1", 
in  "1-5",  describes  this  case  as  the  first  sequential  case  of  this  type.  The  number 
"5",  in  "1-5",  describes  this  case  to  be  grass  at  the  surface  with  a  horizontal  barrier. 
Additional  information  is  shown  in  Appendix  B-121  and  is  described  in  the  process 
flow  sheets  Appendix  A-85,  86. 

5.  Measured  Suction  Profile  Input.  Geotechnical  engineers  may  identify  some 
unique  moisture  effect  cases  that  cannot  be  modeled  in  the  procedures  to  address  the 
typical  moisture  effect  cases  contained  herein.  Those  special  cases  are  addressed 
using  initial  and  final  measured  suction  profiles.  The  input  of  these  measured 
suction  profiles  is  described  in  the  process  flow  sheet  Appendix  A-87,  and  shown  in 
the  program  design  layout  sheet,  Appendix  B-122,  123. 

6.  Time  Input  for  Post  Construction  Cases.  Algorithms  for  post  construction 
solutions  are  time  dependent.  Calculating  the  post  construction  differential  soil 
movements  can  only  be  accomplished  if  the  driest  (or  wettest)  time  of  year  is 
known,  tdry,  the  month  of  year  the  foundation  was  built  is  known,  tconst,  and  the  month 
the  soil  volume  change  is  sought  is  known,  t^.  Input  of  these  quantities  is  described 
in  the  process  flow  sheet  Appendix  A- 107. 

Calculated  Soil  Properties  and  Soil  Characteristics 

Once  the  liquid  limit,  the  plastic  limit,  the  percent  passing  the  number  200  sieve, 
the  percent  smaller  than  two  microns,  and  the  dry  unit  weight  is  entered  for  the  soil 
profile  considered,  the  remainder  of  the  required  soil  properties  and  characteristics  are 
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solved.  The  process  flow  sheet,  Appendix  A-88,  describes  the  algorithm  by  which  these 
remaining  soil  characteristics  are  solved.  Definitions  for  all  symbols  are  found  in  the 
list  of  symbols. 

The  plasticity  index,  PI(Zj),  is  the  difference  between  the  liquid  limit,  LL(Zj),  and 
the  plastic  limit,  PL(Zj).  The  percent  fine  clay,  Clay(Zj),  is  calculated  by  dividing  the 
percent  smaller  than  two  microns,  n(Zj),  by  the  percent  passing  the  number  200  sieve, 
v(Zj).  The  activity  ratio,  Ac(Zj),  is  calculated  by  dividing  the  Plfe)  by  Clay(Zj).  The 
cation  exchange  capacity,  CEC(Zj),  is  the  larger  of  the  values  PL(Zj)1 17  or  LL(Zj)0912. 
The  cation  exchange  activity,  CEAc(Zj),  is  calculated  by  dividing  CEC(Zj)  by  Clay(Zj). 
Using  Fig.  15,  the  matrix  suction  compression  index,  yhl00,  is  the  value  found  at  the 
intersection  of  CEAc(Zj)  and  Ac(Zj).  To  obtain  the  yh(Zj)  simply  multiply  Clay(zi)  by 
yhl00.  The  slope  of  the  suction  versus  the  gravimetric  water  content  line,  Sfe),  is 
approximated  by  (16).  Mitchell's  diffusion  coefficient,  a(Zj),  is  approximated  by  (17). 

Mitchell's  unsaturated  permeability,  K0  (z; )  ■  |h0 1 ,  is  calculated  by  (18).  Using  Fig.  16 
for  center  lift  conditions  and  Fig.  1 7  for  edge  lift  conditions,  the  edge  moisture  variation 
is  found  by  finding  the  curve  that  corresponds  with  the  soil's  diffusion  coefficient,  ct(Zj), 
then  finding  the  intersection  made  by  drawing  a  vertical  line  corresponding  to  the  local 
TMI  (for  locations  in  Texas  the  local  TMI  can  be  found  in  Fig.  26),  and  finally  drawing 
a  horizontal  line  from  that  intersection  across  to  the  y  axis  identifies  the  appropriate 
edge  moisture  variation  distance.  For  each  case  of  center  lift  and  edge  lift,  the  largest 
edge  moisture  variation  value,  considering  all  incremental  soil  layers,  is  the  soil 
profile's  edge  moisture  variation  distance,  emc  and  eme. 

Procedures  to  Calculate  the  Equilibrium  Suction 

The  next  step  in  the  process  of  calculating  the  vertical  differential  soil  movement 
is  to  establish  a  baseline  equilibrium  suction,  h,,,  (units  of  pF).  All  suction  profiles  are 
calculated  relative  to  this  baseline.  The  equilibrium  suction,  h,,,,  is  calculated  for  the  soil 
layer  located  at  the  depth  to  constant  suction.  The  depth  to  constant  suction,  z^,  is 
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dependent  upon  the  h^  Therefore,  if  the  depth  to  constant  suction,  z^,  does  not  fall 
within  the  soil  layer  in  which  the  h,,,  was  calculated,  then  the  h,,,  must  be  re-calculated 
for  the  soil  layer  located  at  the  depth  to  constant  suction.  This  is  an  iterative  process. 
The  process  flow  sheet  which  describes  the  procedure  for  calculating  h,,,  is  found  on 
Appendix  A-90.  The  procedures  for  calculating  z^  are  found  in  the  next  section.  The 
equilibrium  suction  value  is  dependent  upon  the  following  inputs:  TMI,  d^,  h^,  h^,  0S, 
0r,  A,  and  B.  For  the  state  of  Texas  the  TMI  is  obtained  from  Fig.  26.  Values  for  Gs,  0r, 
A,  and  B  for  a  particular  classification  of  soil  are  listed  in  Table  2.  As  reported  by  Gay 
(1994),  a  conservative  estimate  for  the  d^  can  be  taken  as  30  cm  for  all  expansive  soils 
in  which  the  d^  is  not  known.  All  variables  have  default  values.  Refer  to  Appendix  A- 
90.  Once  all  input  is  entered  y,  d„  Tls  0wl,  6dry  are  calculated  using  (25),  (26),  (27),  (38), 
and  (39)  respectively.  These  values  are  used  in  equations  (33)  and  (37)  to  solve  for  h,,,. 
A  value  of  h,,,  is  picked  and  an  iterative  process  for  solving  for  9m  is  performed  until 
equal  values  of  0m  satisfies  both  equations  (33)  and  (37). 

Depth  to  Constant  Suction 

The  depth  to  constant  suction  is  the  depth  in  a  soil  profile  to  which  there  is  no 
longer  a  significant  seasonal  suction  change.  The  definition  of  no  significant  seasonal 
suction  change  is  0.2  pF  from  Lytton's  1994  paper.  Depth  to  constant  suction  is 
dependent  upon  the  diffusivity  of  the  soil,  the  number  of  local  annual  weather  cycles, 
and  the  suction  amplitude  at  the  surface.  As  discussed  in  the  "Soil  Volume  Strain" 
section,  the  soil's  ability  to  change  volume  is  dependent  upon  the  seasonal  change  of 
moisture  and  the  magnitude  of  the  mean  principal  stress  present  to  suppress  the  volume 
change.  Using  Mitchell's  equation  (42)  there  will  always  be  a  reduction  in  the  suction 
amplitude  with  an  increase  of  depth.  Additionally,  an  increase  in  depth  introduces 
added  restraint  to  soil  volume  strain  due  to  the  effect  of  overburden.  Considering  these 
two  aspects  of  the  soil  volume  change  equation  it  becomes  apparent  that  for  depths 
greater  than  30  feet  there  is  usually  no  longer  any  significant  changes  in  suction.  The 
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depth  to  constant  suction,  z^,  will  be  used  as  the  depth  beyond  which  there  is  no 
significant  change  in  suction. 

Table  2.  (Lytton  et  al.  1989)  Gardner's  Coefficients 


Unified 

Void 

Porosity 

Soil 

ratio 

No 

Class. 

e 

AWL 

XWL 

n  =  es 

9=0.1  es 

A 

B 

GM-GC1 

GM-GC 

0.42 

0.004 

0.637 

0.296 

0.0296 

1083.777 

0.637 

GM1 

GM 

0.259 

0.152 

0.269 

0.206 

0.0206 

12.2224 

0.269 

GM2 

GM 

0.382 

0.04 

0.648 

0.276 

0.0276 

111.1578 

0.648 

GM3 

GM 

0.451 

0.066 

0.251 

0.311 

0.0311 

27.00574 

0.251 

GM4 

GM 

0.608 

0.043 

0.478 

0.378 

0.0378 

69.90875 

0.478 

GP1 

GP 

0.255 

0.065 

0.55 

0.203 

0.0203 

54.58668 

0.55 

GW1 

GW 

0.443 

0.039 

0.302 

0.307 

0.0307 

51.39672 

0.302 

GW2 

GW 

0.506 

0.596 

0.318 

0.416 

0.0416 

3.489424 

0.318 

GW3 

GW 

0.252 

0.309 

0.319 

0.355 

0.0355 

6.745925 

0.319 

SM-SC1 

SM-SC 

0.656 

0.013 

0.77 

0.396 

0.0396 

452.9567 

0.77 

SM1 

SM 

0.396 

0.016 

0.562 

0.284 

0.0284 

227.9712 

0.562 

SM2 

SM 

0.619 

0.001 

1.023 

0.382 

0.0382 

10543.87 

1.023 

SM3 

SM 

1.192 

0.11 

0.339 

0.544 

0.0544 

19.843 

0.339 

SM4 

SM 

0.655 

0.039 

0.468 

0.396 

0.0396 

75.32435 

0.468 

SM5 

SM 

0.623 

0.015 

0.835 

0.384 

0.0384 

455.9411 

0.835 

SM6 

SM 

1.024 

0.011 

0.671 

0.506 

0.0506 

426.194 

0.671 

SM7 

SM 

0.72 

0.023 

0.549 

0.419 

0.0419 

153.9119 

0.549 

SM8 

SM 

0.831 

0.208 

0.436 

0.454 

0.0454 

13.12009 

0.436 

SM9 

SM 

0.721 

0.01 

0.835 

0.419 

0.0419 

683.9116 

0.835 

SM10 

SM 

0.818 

0.018 

0.806 

0.45 

0.045 

355.4082 

0.806 

SM11 

SM 

1.132 

0.029 

0.745 

0.531 

0.0531 

191.6911 

0.745 

SM12 

SM 

2.16 

0.042 

0.501 

0.684 

0.0684 

75.46589 

0.501 

SP-SM1 

SP-SM 

0.574 

0.048 

0.769 

0.365 

0.0365 

122.3936 

0.769 

SP-SM2 

SP-SM 

0.583 

0.095 

0.613 

0.368 

0.0368 

43.17938 

0.613 

SP1 

SP 

0.712 

0.0554 

0.79 

0.416 

0.0416 

111.2987 

0.79 

SP2 

SP 

0.736 

0.024 

0.951 

0.424 

0.0424 

372.2106 

0.951 

SP3 

SP 

0.818 

0.042 

0.9 

0.45 

0.045 

189.1258 

0.9 

SP4 

SP 

0.593 

0.076 

0.665 

0.372 

0.0372 

60.83961 

0.665 

SP5 

SP 

0.583 

0.053 

0.809 

0.368 

0.0368 

121.5414 

0.809 

SW-SP1 

SW-SP 

0.773 

0.012 

1.082 

0.296 

0.0296 

1006.512 

1.082 

SW1 

SW 

0.551 

0.162 

0.5852 

0.355 

0.0355 

23.75117 

0.5852 

SW3 

SW 

0.385 

0.111 

0.616 

0.278 

0.0278 

37.21149 

0.616 

CL1 

CL 

1.183 

0.065 

0.417 

0.542 

0.0542 

40.1871 

0.417 

CL2 

CL 

0.621 

0 

0.976 

0.383 

0.0383 

0 

0.976 

CL3 

CL 

0.785 

0.004 

0.61 

0.44 

0.044 

1018.451 

0.61 

CL4 

CL 

0.61 

0 

0.957 

0.379 

0.0379 

0 

0.957 

CL5 

CL 

0.825 

0.018 

0.523 

0.452 

0.0452 

185.2369 

0.523 

CL6 

CL 

0.964 

0 

1.91 

0.491 

0.0491 

0 

1.91 

CL7 

CL 

0.808 

0 

1.593 

0.447 

0.0447 

0 

1.593 

CL8 

CL 

0.972 

0.028 

0.344 

0.493 

0.0493 

78.85731 

0.344 
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Table  2.  Continued 


Unified 

Void 

Porosity 

Soil 

ratio 

No 

Class. 

e 

AWL 

XWL 

n  =  es 

e=o.ies 

A 

B 

CL9 

CL 

0.562 

6.9E-06 

1.551 

0.36 

0.036 

5154077 

1.551 

CL10 

CL 

0.517 

0.024 

0.359 

0.341 

0.0341 

95.23328 

0.359 

CL11 

CL 

0.415 

0.064 

0.247 

0.293 

0.0293 

27.59434 

0.247 

ML8 

ML 

0.858 

0.003 

0.973 

0.462 

0.0462 

3132.411 

0.973 

ML9 

ML 

0.606 

0.022 

0.707 

0.377 

0.0377 

231.514 

0.707 

ML10 

ML 

2.163 

0.038 

0.61 

0.684 

0.0684 

107.2053 

0.61 

ML11 

ML 

0.603 

0.031 

0.722 

0.376 

0.0376 

170.0741 

0.722 

ML12 

ML 

0.821 

0 

2.059 

0.451 

0.0451 

0 

2.059 

CL12 

CL 

0.972 

0.004 

0.474 

0.493 

0.0493 

744.6291 

0.474 

ML-OL1 

ML-OL 

1.721 

0.064 

0.535 

0.632 

0.0632 

53.55747 

0.535 

ML-OL2 

ML-OL 

0.885 

0 

1.362 

0.469 

0.0469 

0 

1.362 

ML-CL1 

ML-CL 

0.692 

0.066 

0.365 

0.409 

0.0409 

35.11204 

0.365 

ML-CL2 

ML-CL 

0.642 

0.013 

0.634 

0.391 

0.0391 

331.1743 

0.634 

ML-CL3 

ML-CL 

0.645 

0.001 

1.052 

0.392 

0.0392 

11271.97 

1.052 

ML1 

ML 

1.015 

0.003 

1.054 

0.504 

0.0504 

3774.668 

1.054 

ML2 

ML 

1.067 

0 

1.257 

0.516 

0.0516 

0 

1.257 

ML3 

ML 

1.119 

0.032 

0.681 

0.528 

0.0528 

149.9167 

0.681 

ML4 

ML 

1.558 

0.065 

0.411 

0.609 

0.0609 

39.63571 

0.411 

ML5 

ML 

0.88 

0 

1.019 

0.468 

0.0468 

0 

1.019 

ML6 

ML 

0.825 

0 

1.474 

0.452 

0.0452 

0 

1.474 

ML7 

ML 

0.783 

0.012 

0.575 

0.439 

0.0439 

313.1978 

0.575 

The  depth  to  constant  suction  varies  dependent  upon  the  soil  strata,  the  local 
weather  conditions,  and  the  type  of  soil  surface  conditions.  The  depth  to  constant 
suction,  z^,  will  be  calculated  for  the  following  cases  of  soil  surface  conditions, 
Appendix  A-9 1 : 

1.  Bare  soil  at  the  surface  (Appendix  A-92). 

2.  Grass  at  the  surface  with  shallow  roots  (Appendix  A-92). 

3.  Flower  bed  at  the  surface  with  a  known  flower  bed  zone  depth,  dzone  (Appendix  A- 
93). 

4.  A  tree  at  the  surface  with  a  known  depth  of  root  zone,  dzone  (Appendix  A-93). 

5.  And  measured  suction  profiles  (Appendix  A-94). 

The  simplest  case  would  be  that  of  bare  soil  or  grass  at  the  surface  with  only  one 
soil  layer  the  entire  depth  of  moisture  active  zone.  For  this  case  equation  (46)  is  all  that 
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is  needed  to  solve  for  z,n.  Simply  replace  the  numerator  in  the  natural  log  function  with 
0.2  pF,  substitute  entered  or  previously  calculated  values  for  all  other  variables  and 
solve  for  z„,.    The  suction  amplitude,  U0,  is  the  difference  between  maximum  and  the 
minimum  surface  soil  suction  values.  For  bare  soils  the  suction  amplitude  is  U0  =  hdry  - 
r^,  where  hdiy  =  6.0  pF.  For  soils  with  grass  at  the  surface  the  suction  amplitude  is  U0  = 
h^  -  h^,  where  hdry  =  4.5  pF.  Refer  to  the  "Typical  Suction  Levels  in  the  Field"  for 
typical  values  of  suction. 

In  most  soil  profiles  there  is  usually  at  least  two  soil  layers  to  consider.  In  cases 
considering  bare  soil  or  grass  at  the  surface,  with  multiple  soil  layers,  the  algorithms  in 
Appendix  A-92  are  used  to  solve  for  z^.  The  procedures  are  simple. 

1 .  First  solve  for  the  ratio  of  the  suction  difference  between  the  top  increment  and 
bottom  increment  of  a  soil  layer,  Ps,  using  equation  (45). 

2.  Multiply  Ps  by  the  suction  amplitude  at  the  top  of  the  layer  which  gives  the  suction 
amplitude  at  the  bottom  of  the  layer,  ULs. 

3.  TestULs: 

a.  If  ULs  is  less  than  0.2  pF,  then  the  depth  to  constant  suction  is  in  this 
layer,  in  which  case  z^  is  solved  using  (46).  The  initial  U0  is  the  same  as 
the  one  described  in  the  previous  paragraph.  All  subsequent  U0's  are 
solved  using  the  procedure  in  paragraph  3.b.  described  below. 

b.  If  ULs  is  grater  than  0.2  pF,  then  set  U0  equal  to  ULs  and  repeat  these  three 
steps. 

The  depth  to  constant  suction  is  the  z^  calculated  for  the  last  layer  added  to  the  value  of 
the  depth  to  the  top  of  this  layer. 

The  procedures  for  calculating  the  z^  for  flower  beds  and  trees  is  basically  the 
same  with  two  distinct  differences:  (1)  The  suction  amplitude  is  the  difference  between 
the  limit  suction  value  and  the  equilibrium  suction  value.  (2)  The  suction  amplitude  at 
the  surface  is  constant  to  the  depth  of  moisture  effect  zone,  dzone.    For  a  flower  bed  at 
the  surface  the  suction  amplitude  is  the  difference  between  the  equilibrium  suction  value 
and  the  wet  limit  suction  value,  U0  =  h,,,  -  h^.  For  a  tree  at  the  surface  the  suction 
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amplitude  is  the  difference  between  the  dry  limit  suction  value  for  vegetation  and 
equilibrium  suction  value,  U0  =  hdry  -  h,n.  Refer  to  Appendix  A-93  for  a  detailed  process 
flow  of  these  procedures.  Refer  to  "Typical  Suction  Levels  in  the  Field"  and 
"Procedures  to  Calculate  the  Equilibrium  Suction"  for  values  of  (r^,  and  hdry)  and  h,^ 
respectively. 

When  making  calculations  of  the  depth  to  constant  suction,  z^,  for  the  case  of 
measured  suction  profiles,  the  depth  to  constant  suction  can  be  taken  as  the  value  of  z^ 
calculated  for  one  of  the  four  previously  explained  cases  if  the  surface  conditions  are  the 
same.    However,  a  check  must  be  made  to  ensure  that  the  z^  in  a  particular  measured 
suction  profile  case  does  not  extend  beyond  those  values  of  z^  in  the  previously 
described  cases.  If  the  z^  extends  beyond  those  values  calculated  using  the  procedures 
described  in  the  previous  four  cases  then  the  z^  must  be  calculated  using  the  procedures 
described  in  Appendix  A-94.  The  procedures  are  the  same  as  those  described  above  for 
bare  soil  and  grass  at  the  surface  with  one  major  difference.  The  difference  is  that  the 
process  is  initiated  at  the  bottom  of  the  measured  suction  profile.  This  z^  value 
calculated  for  a  measured  suction  profile  is  then  compared  to  the  z^  for  the  surface 
conditions  most  similar  to  the  measured  profile  surface  conditions.  The  depth  to 
constant  suction  is  taken  as  the  larger  of  these  values. 
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POST  EQUILIBRIUM  AND  POST  CONSTRUCTION  SOLUTIONS 

At  this  point  in  the  procedures  all  inputs  and  calculated  values  required  to  generate 
suction  profiles  for  soil  volume  strain  calculations  have  been  explained.  The  procedures 
to  generate  soil  suction  profiles  and  the  corresponding  soil  movement  for  moisture  effect 
cases  and  design  effects  are  explained  using  the  following  two  procedural  approaches:  (1) 
the  post  equilibrium  approach,  and  (2)  the  post  construction  approach.  The  post 
equilibrium  approach  refers  to  procedures  to  solve  for  differential  movement  of  expansive 
soils  after  the  soil  under  the  slab  has  reached  an  equilibrium  moisture  content.  When  a 
slab  is  placed  on  an  expansive  soil,  there  is  a  time  period  in  which  the  moisture  content 
under  the  slab  changes  from  the  suction  profile  as  it  was  when  the  slab  was  placed  to 
approach  the  equilibrium  suction  value.  The  algorithms  to  solve  for  differential  soil 
movement  during  this  time  period  are  handled  by  the  post  construction  approach. 

Post  Equilibrium  Approach 

Once  the  soil  under  the  slab  reaches  the  equilibrium  suction  value,  the  differential 
soil  movement  is  calculated  using  the  differential  in  soil  suction  profiles  between  the 
equilibrium  suction  profile  under  the  slab  and  the  suction  profile  at  the  edge  of  the  slab. 
The  suction  profile  at  the  edge  of  the  slab  is  dependent  upon  the  seasonal  moisture  change. 
The  restraint  to  soil  volume  increase  at  any  depth  is  dictated  by  the  amount  of  overburden 
and  surcharge  stress.    The  post  equilibrium  algorithms  to  solve  for  the  limit  cases  are:  (1) 
differential  soil  movement  causing  a  maximum  center  lift  distortion  mode,  and  (2) 
differential  soil  movement  causing  a  maximum  edge  lift  distortion  mode. 

Equilibrium  Suction  Profiles 

The  first  calculations  that  are  made  in  the  post  equilibrium  approach  are  the  base 
line  equilibrium  suction  profiles.  The  equilibrium  suction  profile  is  generated  by  reducing 
the  equilibrium  suction  value  calculated  for  a  specific  moisture  effect  case  by  one 
centimeter  in  suction  for  every  centimeter  of  descent  in  soil  depth.  Care  should  be  taken  to 
convert  the  calculated  equilibrium  suction  value  from  units  of  pF  to  units  of  centimeters. 
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Once  the  equilibrium  suction  profile  is  generated  these  values  of  suction  in  units  of 
centimeters  should  be  converted  back  to  units  of  pF  to  make  further  calculations.  Refer  to 
Appendix  A-96  for  algorithms  to  calculate  equilibrium  suction  profiles. 

Suction  Profiles  Causing  The  Center  Lift  Distortion  Mode 

Center  lift  distortion  modes  are  caused  by  a  differential  in  the  equilibrium  suction 
profile  under  the  slab  and  the  dry  limit  suction  profile  at  the  edge  of  the  slab.  Once  the 
equilibrium  suction  profiles  are  calculated  for  each  case  considered,  the  dry  limit  suction 
profiles  for  the  slab  edge  moisture  effect  conditions  are  calculated.  The  dry  limit  suction 
profiles  for  bare  soil  and  grass  at  the  slab  edge  are  calculated  using  Appendix  A-97,  using 
variables  unique  to  the  specific  moisture  effect  case  and  location.  The  dry  limit  suction 
profiles  for  a  flower  bed  and  a  tree  at  the  edge  of  slab  are  calculated  using  the  Appendix  A- 
98.  The  suction  profiles  are  calculated  by  the  application  of  (43).  The  differences  in  these 
processes,  dependent  upon  moisture  effect  case  and  location  are:  the  surface  suction  values, 
the  soil  properties,  the  baseline  equilibrium  suction  profiles,  and  depth  to  constant  suction 
values.  For  all  moisture  effect  cases,  the  suction  amplitude  change  and  the  depth  are 
initialized  at  every  soil  layer  change.    The  cases  of  a  tree  and  a  flower  bed  at  the  surface 
have  one  other  major  difference,  these  cases  include  the  depth  of  moisture  effect  zone.  The 
suction  profile  within  the  tree  root  or  flower  bed  zone  is  constant  and  equal  to  the  limit 
suction  value  for  the  full  depth  of  the  moisture  effect  zone,  dzone.  Refer  to  Fig.  27  for  an 
example  of  the  dry  limit  post  equilibrium  suction  profiles  for  a  flower  bed  at  the  surface 
with  a  4  foot  moisture  effect  zone  depth.  Refer  to  Fig.  28  for  an  example  of  the  dry  limit 
post  equilibrium  suction  profiles  for  a  tree  at  the  surface  with  20  foot  deep  roots. 

Suction  Profiles  Causing  The  Edge  Lift  Distortion  Mode 

Edge  lift  distortion  modes  are  caused  by  a  differential  in  the  equilibrium  suction 
profile  under  the  slab  and  the  wet  limit  suction  profile  at  the  edge  of  the  slab.  Once  the 
equilibrium  suction  profiles  are  calculated  for  each  case  considered,  the  suction  profiles  for 
the  slab  edge  moisture  effect  conditions  are  calculated.  The  wet  limit  suction  profiles  for 
bare  soil,  grass,  and  trees  at  the  slab's  edge  are  calculated  using  Appendix  A-99.  The  wet 
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limit  suction  profiles  for  a  flower  bed  at  the  edge  of  the  slab  is  calculated  using  the 
algorithm  described  in  Appendix  A- 100.  The  differences  in  these  processes  which  depend 
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Fig.  27.  Post  Equilibrium  Suction  Profiles  for  Flower  Bed  at  the  Surface  (Dry  Limit) 
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Fig.  28.  Post  Equilibrium  Suction  Profiles  for  a  Tree  at  The  Surface  (Dry  Limit) 
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upon  the  moisture  effect  case  and  location  are:  the  soil  properties,  the  baseline  equilibrium 
suction  profiles,  and  depth  to  constant  suction  values.  Note  that  the  suction  profile  is 
constant  and  equal  to  the  limit  suction  value  for  the  full  depth  of  the  flower  bed  zone,  dzone. 
The  suction  amplitude  change  and  the  depth  are  initialized  at  every  soil  layer  change. 
Refer  to  Fig.  29  for  an  example  of  the  wet  limit  post  equilibrium  suction  profiles  for  a 
flower  bed  at  the  surface. 
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Fig.  29.  Post  Equilibrium  Suction  Profiles  for  a  Flower  Bed  at  The  Surface  (Wet  Limit) 


Vertical  Differential  Soil  Movement 

Expansive  soil  volume  change  is  predicted  by  applying  the  theory  explained  in  the 
"Soil  Volume  Strain"  section  of  the  previous  chapter.  The  moisture  effect  cases  considered 
consist  of: 

1 .  Bare  soil  at  the  surface, 

2.  Grass  at  the  surface  with  shallow  roots, 
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3.  A  flower  bed  at  the  surface  with  a  known  depth  of  flower  bed  zone,  dzone,  and 

4.  A  tree  at  the  surface  with  a  known  depth  of  root  zone,  dzone 

The  procedures  also  consider  the  application  of  design  effects  to  provide  the  best  possible 
conditions  so  foundation  design  performance  is  at  a  pinnacle.  The  design  effects 
considered  include: 

1 .  Vertical  barriers, 

2.  Horizontal  barriers,  and 

3.  The  effect  of  locating  the  moisture  effect  case  a  distinct  distance  from  the  slab  edge. 
Explanation  of  the  procedures  to  predict  vertical  differential  soil  movement  is  broken  down 
into  three  categories:  (1)  all  moisture  effect  cases  with  or  without  a  vertical  barrier,  (2) 
moisture  effect  cases  with  or  without  a  horizontal  barrier,  and  (3)  moisture  effect  cases 
located  a  distinct  distance  from  the  slab  edge. 

Post  equilibrium  vertical  differential  soil  movement  can  be  determined  for  all 
moisture  effect  cases  with  or  without  a  vertical  barrier  by  simply  applying  the  algorithms 
described  by  Appendix  A- 103.  This  process  can  be  broken  down  into  three  basic  steps: 
(1)  All  pertinent  variables  and  information  collected  or  calculated  in  previous  procedural 
steps  are  retrieved  for  use  in  this  process.  (2)  Tests  are  performed  on  variables  to  ensure 
the  appropriate  calculations  and  quantities  are  used,  based  on  Lytton's  volume  strain 
theory.  (3)  Incremental  strain  and  corresponding  vertical  soil  movement  is  summed  up  for 
every  soil  increment  within  the  moisture  active  zone  of  the  soil  profile.  The  total  vertical 
differential  soil  movement  quantity  is  defined  by  the  variable  Az.    Refer  to  Fig.  30  for  an 
example  of  bare  soil  dry  limit  suction  profiles  (center  lift  distortion  mode)  for  a  slab  edge 
that  has  a  four  foot  deep  vertical  barrier. 

Post  equilibrium  vertical  differential  soil  movement  can  be  computed  for  all 
moisture  effect  cases,  except  the  tree  case,  with  or  without  a  horizontal  barrier  by  simply 
applying  the  algorithms  described  by  Appendix  A- 104.  The  case  of  the  tree  is  unique  in 
that  the  tree  roots  will  be  active  regardless  of  whether  a  horizontal  barrier  exists  or  not.  If 
the  tree  is  killed  because  of  the  placement  of  the  horizontal  barrier  then  the  point  becomes 
mute.  The  process  in  Appendix  A-104  is  broken  down  into  four  basic  steps:  (1)  All 
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Fig.  30  Bare  Soil  Dry  Limit  Suction  Profiles  for  a  Slab  with  a  4ft  Vertical  Barrier 


pertinent  variables  and  information  collected  or  calculated  in  previous  procedural  steps 
are  retrieved  for  use  in  this  process.  (2)  The  horizontal  velocity  for  each  incremental  soil 
layer  is  calculated  based  on  the  edge  moisture  variation  distance  and  the  size  of  horizontal 
barrier.  This  step  is  based  on  equation  (52).  (3)  The  suction  at  the  slab  edge  is  calculated 
using  equation  (53).  (4)  This  slab  edge  suction  value,  hedge(Zj),  is  used  in  the  place  of  the 
final  suction  value  in  Appendix  A- 103  to  make  vertical  differential  soil  movement 
calculations,  Az.  Refer  to  Fig.  3 1  for  an  example  of  wet  limit  suction  profiles  (edge  lift 
distortion  mode)  for  a  four  foot  deep  flower  bed  at  the  edge  of  a  four  foot  horizontal  barrier 
constructed  adjacent  to  the  edge  of  the  slab. 

If  a  flower  bed  or  tree  is  located  a  distinct  distance  from  the  slab  edge,  the  vertical 
differential  soil  movement  is  calculated  by  applying  Appendix  A- 104.  These  are  the  only 
two  cases  that  can  be  located  a  distance  from  the  slab  and  still  have  an  effect  on  vertical 
differential  soil  movement  under  the  slab.  The  differential  movement  caused  by  the  tree  or 
flower  bed  located  a  distance  from  the  slab  edge  should  be  compared  to  the  differential 
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movement  caused  by  the  surface  condition  that  exists  between  the  slab  edge  and  the 
tree/flower  bed.  The  maximum  differential  soil  movement  produced  by  either  the 
tree/flower  bed  case  or  the  conditions  between  the  slab  edge  and  the  tree/flower  bed  case 
controls.  This  process  is  broken  down  into  four  basic  steps:  (1)  All  pertinent  variables  and 
information  collected  or  calculated  in  previous  procedural  steps  are  retrieved  for  use  in  this 
process.  (2)  The  horizontal  velocity  for  each  incremental  soil  layer  is  calculated  based  on 
the  edge  moisture  variation  distance  and  the  distance  the  moisture  effect  is  away  from  the 
slab  edge,  dsource  This  step  is  based  on  equation  (52).  (3)  The  suction  at  the  slab  edge  is 
calculated  using  equation  (53).  (4)  This  slab  edge  suction  value,  hedge(Zj),  is  used  in  the 
place  of  the  final  suction  value  in  Appendix  A- 103  to  make  vertical  differential  soil 
movement  calculations,  Az.  Refer  to  Fig.  32  for  a  dry  limit  suction  profiles  (center  lift 
distortion  mode)  for  tree  roots  that  are  20  ft  deep  and  located  four  feet  away  from  the  edge 
of  a  slab  with  a  four  foot  vertical  barrier. 
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Fig.  32.  Suction  Profiles  for  a  Tree  4ft  From  the  Slab  Edge  with  a  4ft  Vertical  Barrier 

Post  Construction  Approach 

As  described  in  the  "Post  Construction  Theory"  section  earlier,  the  worst  case  of 
vertical  soil  movement  in  expansive  soils  may  occur  in  the  post  construction  time  period  of 
the  life  of  a  structure.  Lytton  and  Woodburn's  1973  paper  provides  a  perfect  illustration  of 
this  worst  case  scenario  occurring  in  the  post  construction  time  period  during  the  life  of  a 
school  building.  The  procedures  to  apply  the  post  construction  theory  is  basically  a  two 
step  process.  First,  all  suction  profiles  considered  are  calculated.  Second,  these  suction 
profiles  are  used  to  calculate  vertical  soil  movement  at  the  center  of  the  slab,  Azslab,  and 
combined  with  the  vertical  soil  movement  at  the  slab  edge,  Azedge,  to  arrive  at  a  total 
vertical  differential  soil  movement  quantity,  Az,otal. 


Post  Construction  Suction  Profiles 

The  suction  profiles  must  be  calculated  for  two  time  periods:  soil  suction  profiles 
at  the  time  of  construction,  tconst,  and  soil  suction  profiles  at  the  time  the  vertical  differential 
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soil  movement  is  sought,  t^.  Soil  suction  profiles  at  the  edge  of  the  slab  are  generated  for 
the  time  at  which  the  vertical  differential  soil  movement  is  sought.  Soil  suction  profiles 
beneath  the  slab  are  generated  at  the  time  of  construction  and  at  the  time  the  vertical 
differential  soil  movement  is  sought. 

Appendix  A- 107  is  the  process  flow  sheet  that:  describes  the  time  variables,  lists  a 
step  by  step  process  for  entering  the  required  times,  and  lists  the  steps  to  calculate  the  time 
periods  needed  to  make  soil  suction  profile  calculations.  Appendix  B-127  is  a  picture  of 
the  program  design  layout  sheet  for  the  time  input  and  display  of  suction  profiles  both 
under  the  slab  and  at  the  edge  of  the  slab  for  time  that  transpires  from  the  driest  month  to 
the  month  the  volume  change  is  sought.  In  Appendix  B-127,  the  suction  profile  under  the 
slab  at  the  time  of  construction  is  the  thin  solid  line.  The  dashed  line  is  the  suction  profile 
under  the  slab  six  months  after  construction.  Notice  that  this  dashed  line  is  tending  toward 
equilibrium.  The  equilibrium  profiles  are  the  vertical  thin  lines  with  tick  marks.  Finally, 
the  thick  solid  line  is  the  suction  at  the  edge  of  the  slab  with  a  vertical  barrier.  To  generate 
suction  profiles  illustrated  in  Appendix  B-127  the  following  time  inputs  are  required:  (1) 
the  driest  month  of  the  year,  tdry,  (2)  the  month  the  slab  is  constructed,  tconst,  and  (3)  the 
month  the  soil  volume  change  is  sought,  t^.  These  entered  values  of  time  are  required  for 
two  primary  purposes:  (1)  to  establish  a  baseline  time  to  enable  application  of  the  surface 
annual  weather  cycle  model  developed  by  Mitchell,  (40),  and  (2)  to  establish  time  periods 
used  in  the  procedures  to  calculate  suction  profiles  for  each  case  at  times  required  to  make 
calculations  for  the  post  construction  algorithms.  Appendix  A- 1 07  is  broken  down  into  two 
basic  steps.  First,  the  times  required  for  input  are  entered.  Second,  three  time  periods  are 
calculated  using  the  driest  month  of  the  year  as  the  baseline.  The  time  period  used  to 
calculate  the  suction  profile  under  the  slab  at  the  time  of  construction  is  represented  by  the 
variable  tconst/dry,  as  calculated  in  Appendix  A- 107.  The  time  period  between  the  time  of 
construction  and  the  time  at  which  the  soil  volume  change  is  sought  is  represented  by  the 
variable  t^^,  as  calculated  in  Appendix  A- 107.  Finally,  the  variable  t^^  describes  the 
time  transpired  from  the  driest  month  of  the  year  to  the  month  the  soil  volume  change  is 
sought. 
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All  soil  suction  profiles  with  respect  to  surface  conditions,  time  and  depth  are 
calculated  using  the  algorithms  in  Appendix  A- 109  and  Appendix  A-l  10.  The  process  by 
which  post  construction  soil  suction  profiles  are  generated  consist  of  the  following  basic 
steps: 

1 .  All  values  required  to  make  calculations  are  retrieved  for  each  particular  moisture 
effect  case  and  time  the  suction  profile  is  sought,  Appendix  A- 109. 

2.  Suction  values,  for  each  particular  moisture  effect  case,  are  calculated  for  5  centimeter 
thick  incremental  soil  layers. 

3.  For  every  incremental  soil  layer  cycle,  variables  are  re-initialized  and  pertinent 
equations  are  activated  based  on:  the  location  within  the  soil  profile,  the  moisture  effect 
case  type,  and  the  type  and  size  of  design  effect  implemented  (If  a  horizontal  barrier  is 
implemented  or  if  a  moisture  effect  case  is  located  a  distance  from  the  edge,  Appendix 
A- 104  will  be  used  to  augment  the  procedures  described  in  Appendix  A-l  10.) 

Once  all  soil  suction  profiles  are  generated,  the  soil  suction  profile  at  the  slab  center 
is  damped  from  time  tconst  to  time  t^  by  applying  the  algorithms  described  by  Appendix  A- 
111.  This  process  takes  the  soil  suction  profile  under  the  slab  at  the  time  of  construction, 
h.siib(zi,tconstJliry),  and  dampens  it  based  on  equations  (59)  and  (61).  These  equations  model 
the  soil  suction  profile  under  the  center  of  the  slab  at  time  tconst  so  that  the  soil  suction  in 
each  incremental  soil  layer  approaches  the  equilibrium  suction  profile.  The  rate  of 
approach  depends  on  unsaturated  soil  profile  properties,  surface  weather  conditions,  and 
time.  The  process  in  Appendix  A-l  1 1  is  broken  down  into  two  basic  steps.  First,  all 
required  variables  are  retrieved.  Second,  the  dampened  suction  value  for  every  incremental 
soil  layer  is  calculated  for  the  time  transpired  from  the  time  of  construction  to  the  time  the 
soil  volume  change  is  sought. 

Vertical  Differential  Soil  Movement 

Once  all  soil  suction  profiles  are  calculated,  the  expansive  volume  strain  theory 
developed  by  Lytton  is  applied.  All  moisture  effect  cases  are  applied  to  the  same  vertical 
differential  soil  movement  algorithms  which  are  broken  down  into  a  four  step  process. 
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First,  all  variables  and  quantities  required  for  calculations  are  retrieved,  as  described  in 
Appendix  A-l  13,  depending  on  which  moisture  effect  case  and  design  effect  scenario  is 
chosen.  Second,  vertical  differential  soil  movement  is  calculated  for  soil  at  the  edge  of  the 
foundation  by  way  of  steps  described  in  Appendix  A-l  14.  Third,  vertical  differential  soil 
movement  is  calculated  for  soil  under  the  center  of  the  foundation  by  way  of  steps 
described  in  Appendix  A-l  15.  Finally,  the  two  quantities  of  vertical  differential  soil 
movement  are  summed  together  to  arrive  at  the  total  differential  soil  movement  for  the  soil 
under  the  slab  at  a  specified  time  after  construction,  as  described  in  Appendix  A-l  15. 

Appendix  B-128  is  a  program  design  layout  sheet  that  illustrates  the  results  of  the 
post  construction  procedures  when  applied  to  calculate  the  worst  case  center  lift  mode  for 
case  "1-6".  The  moisture  effect  case  considered  in  Appendix  B-128  is  grass  at  the  surface 
with  a  4  ft  vertical  barrier,  case  "1-6".  The  vertical  soil  movement  calculated  at  the  edge 
moisture  variation  distance  under  the  slab  is  Azs,ab  =  -0.13  inches  (shrink).  The  vertical  soil 
movement  calculated  at  the  edge  of  the  slab  is  Azedge  =  -0.61  inches.  The  total  vertical 
differential  soil  movement  for  this  moisture  effect  case  using  Table  2  as  the  soil  profile  is 
Az,ota|  =  0.48  inches  (shrink).  This  quantity,  Az,otal ,  is  the  difference  between  the  values  of 
Azslab  and  Azedge. 
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MEASURED  SUCTION  PROFILES 

When  cases  exist  that  can  not  be  modeled  by  one  of  the  moisture  effect  cases 
discussed  in  the  previous  procedures,  vertical  differential  soil  movement  can  be  calculated 
by  the  measured  suction  profile  method.  This  is  a  method  to  solve  for  the  differential  soil 
movement  by  measuring  the  initial  and  final  suction  profiles  and  applying  those  profiles  to 
the  volume  strain  theory  developed  by  Lytton.  The  procedures  are  broken  down  into  three 
broad  steps:  (1)  The  measured  suction  values  for  a  sample  are  entered  and  stored  for  later 
use,  Appendix  A-87.  (2)  The  entered  measured  suction  values  for  a  sample  are  converted 
into  a  measured  suction  profile,  Appendix  A-101.  (3)  The  measured  suction  profiles  are 
used  to  solve  for  vertical  differential  soil  movement  for  a  specific  column  of  soil, 
Appendix  A- 105  and  Appendix  A- 106. 

The  first  step  to  applying  this  methodology  is  to  input  measured  suction  values  for 
a  sample  of  soil.  There  are  a  number  of  methods  to  measure  the  matrix  suction  in  a  soil 
sample.  One  of  the  more  accurate  means  for  measuring  total  soil  suction  is  the 
thermocouple  psychrometer.  The  thermocouple  psychrometer  sensor  and  data  logging 
equipment  is  not  inexpensive  and  requires  some  expertise  and  maintenance  to  operate.  For 
these  reasons  the  filter  paper  is  becoming  a  more  widely  accepted  method  in  practice.  The 
filter  paper  is  inexpensive,  requires  practically  no  maintenance,  and  requires  very  little 
expertise  to  obtain  results  nearly  as  accurate  as  the  thermocouple  psychrometer  (McKeen 
1981). 

The  process  of  inputting  measured  suction  values  occurs  as  part  of  the  process  to 
input  moisture  effect  cases.  Appendix  B-121  is  a  program  design  layout  sheet  that 
illustrates  the  choices  for  the  moisture  effect  case  input.  The  program  user  would  select 
"12.  Measured  Suction  Profile"  and  the  next  screen  that  appears  is  illustrated  by  Appendix 
B-122.  Appendix  B-122  and  123  are  a  program  design  layout  sheets  that  illustrate 
windows  environment  screens  for  inputting  the  measured  suction  profiles.  Appendix  A-87 
is  the  process  flow  describing  measured  suction  profile  input.  It  is  important  to  note  that  a 
surface  suction  value  must  be  entered  or  a  default  value  equal  to  the  equilibrium  suction 
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will  be  assumed. 

Appendix  A- 101  is  the  process  flow  describing  the  steps  to  generate  a  measured 
suction  profile  based  on  the  entered  suction  values.  The  procedures  are  broken  down  into 
two  basic  steps:  (1)  The  slope  of  a  line  connecting  adjacent  suction  values  within  a 
sample  is  calculated.  (2)  The  measured  suction  profile  is  generated,  using  this  previously 
calculated  slope,  by  connecting  a  line  to  adjacent  suction  values  in  the  sample.  Table  3  is 
an  example  of  values  used  to  generate  a  measured  suction  profile.  Fig.  33  is  a  graph  of  the 
initial  and  final  measured  suction  profiles  generated  by  the  values  listed  in  Table  3 

Table  3  Measured  Suction  Profile  Values 


Input  Suction  Values 


Initial 


Final 


(PF) 
Suction 


(cm) 
Suction 


31622.78 
19952.62 
10000 
5011.872 
3715.352 


(cm) 
Depth 


0 

91.44 

365.76 

487.68 

609.6 


(pF)  (cm) 

Suction     Suction 


2.5  316.2278 

3  1000 

3.3  1995.262 

3.5  3162.278 

3.57  3715.352 


Given 

Values 

Zm  = 

20  ft 

hm  = 

3.57  pF 

a  = 

40  cm 

•V>  (dry)  = 

0.33 

■V>  (wet)  = 

0.67 

■  (dry)  = 

0.5 

■  (wet)  ~ 

0.8 

Applying  the  measured  suction  profiles  to  predict  vertical  differential  soil 
movement  is  probably  the  simplest  of  all  procedures  contained  herein.  The  procedures  are 
broken  down  into  two  basic  steps:  (1)  select  an  initial  and  final  suction  profile  to 
accurately  model  the  moisture  change  scenario  being  considered,  Appendix  A- 105,  and  (2) 
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apply  the  algorithms  to  calculate  vertical  differential  soil  movement  based  on  these 
measured  suction  profiles,  Appendix  A- 106. 


Fig.  33.  Initial  and  Final  Measured  Suction  Profiles  Described  in  Table  4 


The  challenge  is  in  deciding  what  suction  profiles  are  used  as  the  initial  and  final 
suction  profile  to  accurately  model  the  soil  conditions  which  produce  reasonable  limit 
conditions  used  for  design.  The  answer  should  come  from  knowledge  about  history  of  the 
site,  to  include:  surface  drainage  conditions,  soil  characteristics  in  the  soil  profile,  surface 
weather  conditions,  and  type  of  structure  placed.  For  example,  consider  the  construction 
site  described  by  the  following:  the  site  is  located  in  College  Station,  the  site  is  surface 
drained,  the  soil  profile  is  described  by  Table  1,  and  the  measured  suction  profile  is  located 
where  a  large  oak  tree  was  removed  as  a  result  of  pre-construction  clearing.  The  measured 
suction  profile  at  this  location  is  a  dry  profile.  This  profile  will  tend  toward  equilibrium 
and  may  approach  the  wet  limit  state  subsequent  to  a  wet  season.  A  conservative  approach 
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for  calculating  vertical  differential  soil  movement  is  to  set  the  initial  suction  profile  equal 
to  the  measured  suction  profile  and  set  the  final  suction  profile  equal  to  the  wet  limit  state 
for  the  surface  conditions  that  exist  after  the  construction  project  is  complete.  This 
scenario  produces  heave  values  of  vertical  differential  soil  movement.  Table  3  and  Fig.  33 
are  values  used  to  obtain  results  for  vertical  differential  soil  movement  in  the  example 
described  above.  The  vertical  differential  soil  movement  in  this  case  is  5.48  inches  of 
heave. 
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DISCUSSION  OF  RESULTS 

The  procedures  to  predict  vertical  differential  soil  movement  were  developed  based 
on  the  breakthroughs  mentioned  earlier.  These  procedures  consider  moisture  effect  cases 
that  are  common  with  light  commercial  and  residential  structures,  such  as:  (1)  bare  soil  at 
the  surface,  (2)  grass  at  the  surface,  (3)  tree  at  the  surface,  and  (4)  flower  bed  at  the  surface. 
Additionally,  these  procedures  include  calculating  the  effects  of  differential  soil  movement 
caused  by  the  introduction  of  design  effects,  such  as,  vertical  barriers  and  horizontal 
barriers.  The  methodology  to  handle  these  moisture  effect  cases  and  design  effects  are 
quite  delineated  in  detail  in  the  process  flow  of  detailed  algorithms,  Appendix  A.  A 
graphical  representation  of  the  program  design  layout  is  illustrated  in  Appendix  B. 

The  algorithms  in  Appendix  A  were  programmed  in  a  spreadsheet  software 
package.  The  author  used  this  package  to  analyze  the  moisture  effect  cases  illustrated  in 
Appendix  C  which  consists  of  the  soil  profile  described  in  Table  1.  Table  4  is  a  summary 
of  vertical  differential  soil  movements  for  the  site  described  by  Appendix  C.  This 
summary  illustrates  the  worst  cases  of  edge  lift  and  center  lift  distortion  modes  for  all 
moisture  effect  cases.  Post  equilibrium  and  post  construction  approaches  were  used  to 
solve  for  all  cases  so  that  the  worst  possible  case  of  vertical  differential  soil  movement  is 
identified.    Additional  analyses  were  performed  using  different  sizes  of  vertical  and 
horizontal  barriers  to  define  the  effects  that  these  design  features  have  on  soil  movement. 
An  engineer  may  use  the  maximum  cases  of  edge  lift  and  center  lift  caused  by  the 
combination  of  moisture  effect  cases  and  the  desired  design  effects  to  design  an  adequate 
foundation. 

Table  4  is  set  up  in  chronological  order  beginning  with  the  post  construction 
differential  soil  movements  and  concluding  with  the  post  equilibrium  soil  movements. 

The  post  construction  numbers  in  Table  4  were  generated  for  June  as  the  driest 
month  of  the  year  and  for  a  slab  built  in  January  (the  wettest  month  of  the  year).  The  slab 
is  placed  on  "wet"  soil  and  this  soil  is  allowed  to  approach  an  equilibrium  moisture  content 
based  on  the  soil  characteristics  and  the  post  construction  theory  described  earlier.  The  soil 
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at  the  edge  of  the  slab  continues  to  fluctuate  with  the  seasonal  weather  changes.  Applying 
this  post  construction  scenario  to  the  "  1  -7  Tree"  case  produces  the  maximum  differential 
soil  movement  compared  to  all  other  cases  which  cause  a  center  lift  distortion  mode.  The 
differential  soil  movement  in  this  case  is  4.32  inches  and  occurs  the  first  August  after 
construction.  Note  that  a  reduction  to  2.44  inches  of  differential  soil  movement  is  made 
possible  by  introducing  a  four  foot  vertical  barrier  at  the  edge  of  the  slab  where  the  tree  is 
located. 

Table  4.  Summary  of  Vertical  Differential  Soil  Movements  for  Appendix  C 

Post  Construction 


Driest  Month  of  the 
year 


Month  of 
Construction. 


June 


January 


Worst  Case  Center  Lift  Distortion  Mode 


Moisture  Effect  Case/Design  Scenario 


1-4  Grass 

1-5  Grass  with  4  ft  Horizontal  Barrier 

1-6  Grass  with  4  ft  Vertical  Barrier 

1-7  Tree 

1-8  Tree  with  4  ft  Vertical  Barrier 

1-9  Flower  Bed 


Differential  Soil  Movement 

Number  of 

(units  in  inches) 

Months  after 

Const.  Worst 

Case  Occurs 

2.34 

6 

1.19 

6 

0.48 

7 

4.32 

7 

2.44 

7 

3.30 

6 

Worst  Case  Edge  Lift  Distortion  Mode 


Moisture  Effect  Case/Design  Scenario 


1-4  Grass 

1-5  Grass  with  4  ft  Horizontal  Barrier 

1-6  Grass  with  4  ft  Vertical  Barrier 

1-7  Tree 

1-8  Tree  with  4  ft  Vertical  Barrier 

1-9  Flower  Bed 

1-1 1  Flower  Bed  with  4  ft  Vertical  Barrier 


Differential  Soil  Movement 

Number  of 

(units  in  inches) 

Months  after 

Const.  Worst 

Case  Occurs 

0.43 

11 

0.31 

11 

0.23 

14 

0.43 

11 

0.23 

14 

0.92 

12 

0.20 

12 

Edqe  Heave  (in) 

Edqe  Shrink 

Ml 

1.65 

1.03 

0.47 

0.56 

0.07 

0.09 

1.64 

3.67 

0.07 

2.36 

2.89 

1.86 

0.68 

0.55 

4.40 

NA 

4.40 

NA 

1.62 

NA 
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Table  4.  Continued 

Post  Equilibrium 

Worst  Case  Differential  Soil  Movement 

Moisture  Effect  Case/Design  Scenario 

1-4  Grass 

1-5  Grass  with  4  ft  Horizontal  Barrier 

1-6  Grass  with  4  ft  Vertical  Barrier 

1-7  Tree 

1-8  Tree  with  4  ft  Vertical  Barrier 

1-9  Flower  Bed 

1-11  Flower  Bed  with  4  ft  Vertical  Barrier 

1-12  Measured  Suction  Profile  (Point  Water  Source) 

2-12  Measured  Suction  Profile  (Point  Water  Source) 

3-12  Measured  Suction  Profile  (Tree  Removal) 

An  inspection  of  Table  4  shows  that  the  differential  soil  movement  is  greatest  when 
there  are  no  vertical  or  horizontal  barriers  in  place.  In  the  "1-4  Grass"  case,  the  edge 
heave  is  1 .65  inches  and  edge  shrink  is  1 .03  inches.  When  a  four  foot  horizontal  barrier  is 
added  the  edge  heave  is  reduced  to  0.47  inches  and  the  edge  shrink  is  reduced  to  0.56 
inches,  as  described  in  the  Table  4  case  "1-5  Grass  with  4  ft  Horizontal  Barrier".  A  further 
reduction  in  soil  movement  is  obtained  by  adding  a  four  foot  vertical  barrier  as  shown  in 
Table  4  case  "1-6  Grass  with  4  ft  Vertical  Barrier".  This  latter  case  reduces  the  edge  heave 
to  0.07  inches  and  the  edge  shrink  to  0.09  inches. 

As  shown  in  Table  4,  the  worst  case  edge  heave  occurs  in  the  post  equilibrium  case 
of  the  measured  suction  profiles  produced  at  the  locations  of  the  leaking  pipes.  The 
leaking  pipes  produced  a  "wet"  spot  that  caused  the  soil  to  heave  4.40  inches  at  those  two 
locations.  This  type  of  heave  at  the  edge  caused  the  worst  case  edge  lift  distortion  mode. 
Refer  to  Appendix  B-130  for  a  conceptual  illustration  of  the  post  construction  case  of 
vertical  differential  soil  movement  that  can  be  predicted  by  applying  the  procedures 
contained  herein.    These  procedures  for  analyzing  vertical  differential  soil  movements 
should  be  applied  to  a  soil  profile  to  enable  an  engineer  to  maximize  resources  ensuring 
that  the  best  foundation  is  designed. 
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SUMMARY  AND  SUGGESTIONS  FOR  RESEARCH  AND  DEVELOPMENT 

The  theory  to  predict  differential  soil  movement  of  expansive  soils  has  evolved 
since  the  1970's.  The  engineering  profession  is  at  a  juncture  today  in  which  engineers 
have  the  technology  to  apply  this  information  to  maximize  resources  and  produce  high 
performing  ground  structures.  The  evolution  of  this  theory  began  in  the  early  1970's  with 
the  advent  of  the  volume  strain  theory  for  expansive  soils  which  was  developed  by  Lytton 
(1973).  In  the  decade  of  the  1970's,  great  strides  were  made  by  Pearring  (1968),  Holt 
(1969),  and  Mojekwu  (1979)  to  enable  engineers  to  identify  predominant  clays.  In  1981, 
McKeen  used  these  newfound  procedures  to  identify  predominant  clays  and  developed  a 
correlation  to  link  groups  of  clays  to  suction  compression  indices.    The  suction 
compression  index  is  used  in  the  soil  volume  strain  equations  developed  by  Lytton  to 
predict  vertical  differential  soil  movement  for  particular  types  of  expansive  soils.  In  the 
1980's,  although  the  volume  strain  theory  was  well  based,  its  application  was  limited  by 
the  lack  of  knowledge  and  methodology  to  predict  the  equilibrium  suction  and  soil  suction 
profiles  for  a  particular  soil  profile  and  location.  This  is  when  work  by  Mitchell  and  later 
Gay  opened  the  door  for  major  advancements  in  this  area  of  expansive  soils.  In  1980, 
Mitchell  developed  and  applied  simple  mathematical  methods  for  predicting  soil  suction 
profiles.  Then  in  1 994,  Gay  developed  a  way  to  predict  the  mean  volumetric  water  content 
for  soils  dependent  upon  the  location  and  climatic  conditions.  Jayatilaka,  Gay,  Lytton,  and 
Wray  (1992)  then  developed  and  calibrated  to  field  data  a  finite  element  method  for 
determining  the  edge  moisture  variation  distance.  They  also  devised  methods  to  predict  a 
number  of  unsaturated  soil  properties  using  simple  index  properties. 

Based  on  the  wealth  of  information  and  breakthroughs  in  the  areas  of  expansive 
soils  listed  above,  Lytton  and  the  author  made  five  critical  advancements  necessary  to 
effectively  and  easily  apply  this  theory: 

1 .  The  means  to  calculate  the  equilibrium  suction  for  a  particular  soil  profile  and  location. 

2.  The  means  to  calculate  depth  to  constant  suction  for  a  soil  profile  with  multiple  layers. 

3.  The  equation  to  calculate  horizontal  velocity  of  water  flow  in  unsaturated  soils. 
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4.  A  method  to  consider  the  transient  case  for  calculating  vertical  differential  soil 
movement  for  expansive  soils  immediately  after  construction  and  before  the  soil  under 
the  center  of  the  slab  has  reached  an  equilibrium  moisture  content. 

5.  Procedures  to  apply  all  of  the  theory  discussed  herein  to  predict  differential  soil 
movement  for  expansive  soils  considering  all  of  the  moisture  effect  cases  and  design 
effects  that  have  been  mentioned. 

The  theory  discussed  herein  represents  physical  behavioral  characteristics  of 
expansive  and  unsaturated  soils  and  is  well  documented.  Those  contributing  to  this 
vertical  differential  soil  movement  theory  have  approached  the  problem  with  a 
theoretically  sound  yet  practical  mindset.  Additionally,  this  expansive  soil  theory  has  been 
developed  in  a  way  to  enable  practicing  engineers  to  run  simple,  common  and  economical 
geotechnical  engineering  lab  tests  to  adequately  define  all  necessary  soils  profile 
information.  This  soil  profile  information  is  applied  in  the  vertical  differential  soil 
movement  theory  and  should  be  used  to  make  sound  geotechnical  engineering  designs. 

Further  development  needs  to  be  made  in  the  areas  of  enhancements  to  the  present 
programs,  materials  characterization  testing,  in  situ  testing,  and  foundation  structural 
analysis  models  for  stiffened  plates  and  drilled  shafts  in  expansive  soils.  The  program 
enhancements  should  include  the  choice  to  use  the  Gardner,  as  well  as  the  Mitchell, 
unsaturated  permeability  options.  Materials  characterization  testing  developments  that  are 
needed  are  more  convenient  methods  of  measuring  the  volume  change,  diffusion,  and  soil- 
water  characteristics  of  expansive  soils  directly.  In  situ  testing  developments  that  are 
needed  include  methods  of  measuring  suction  profiles  on  site  and  preferably  down  hole; 
lateral  earth  pressure  coefficient;  and  of  detecting  the  presence  of  roots  and  voids,  either 
air-filled  or  water-conducting;  in  the  field  and  non-destructively.  Stiffened  plate  models 
which  incorporate  non-rectangular  geometry  need  to  be  coupled  with  this  program  or  its 
enhancements  to  determine  design  quantities  of  moment,  shear,  and  deflection.  Drilled 
shaft  models  with  interface  elements  need  to  be  coupled  with  the  enhanced  program  to 
determine  the  design  quantities  for  both  shrinking  and  expanding  soils. 
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APPENDIX  A 

PROCEDURES  TO  CALCULATE 
VERTICAL  DIFFERENTIAL  SOIL  MOVEMENT 
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Start  of  Program 

and 

Project  Input 

[)  by  Donald  Naiser,  Jr.  5  Aug  97 

-d2 


Note  to  Programmer:  This 
information  should  be  stored  in  a 
database  file  of  "Projects"  in  a 
spreadsheet  type  format.  This 
"Projects"  file  could  then  be  used 
to  organize  all  projects  analyzed 
by  this  software  package. 
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Foundation  Layout  Input 


General  Note 


Note  to  user:  A  three  dimensional  rectangular  (x,y,z) 
coordinate  system  is  used  to  define  foundation  layout.  The  Y 
direction  is  defined  by  an  angle  wrt  true  north.  The  z  direction  is 
depth  of  soil  (increasing  positive  below  surface).  Every  part  of  the 
foundation  structure  must  be  in  the  positive  (x,y)  plane  and  the  slab 
shape  must  be  rectangular.  The  origin  shall  be  defined  by  latitude, 
longitude,  and  elevation  Refer  to  Appendix  B-118, 125,  and  126 
for  graphical  representation. 


Input  the  slab  origin 
(typically  the  west- 
most  and  the  south- 
most  corner  of  the 
slab. 


Specific  Notes 


Inputs  for  latitude  and  longitude 

are  intended  for  data  collection 
purposes  and  to  link  the  site  to  site 
specific  data  such  as;  Thornthwaite 
Moisture  Input  data,  and  edge 
moisture  variation  data. 
Refer  to  the  example  below. 


Input  the  elevation  of 
the  origin. 


The  origin  elevation  is  simply  used  as 
a  reference  for  3D  computations. 
Refer  to  the  example  below. 


This  has  no  bearing  on  calculations. 

This  simply  allows  the  user  to  establish 
building  location/alignment.  The  value 
should  be  entered  clockwise  from  true 
north.  Refer  to  the  example  below. 


Input  the  Local  TMI 
TMI  = 


Input  the  number  of 
Local  weather  cycles/yr 
n  = 


Input  Origin:   Latitude 

|29°1B'N 

Angle  Clok wis e  from 
True  North 

Longitude 
Origin  Elevation 

|  96°28'  W 

|  100.00 

Input  Direction  of  "Y*  Axis 

1  4°  30' 

Input  the  information 
for  the  corners  of  the 
slab. 


Example:  Cut-out  1  from  Appendix  B-118 


Specific  Note 


All  cells  in  the  "corner  information" 
spreadsheet  are  critical.  The  user 
must  not  be  allowed  to  continue  without 
first  entering  this  information. 
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Soil  Layer  Property 
Input 


Input:  Number  of  boring  samples 
B  = 


-Yes- 


S   =S  +  1 


Input  Soil  Layer  Properties: 

Ls= Depth  to  next  layer(ft) 

LLs  = Liquid  Limit  for  layer  s  (%) 

PLS= Plastic  Limit  for  layer  s  (%) 

|is= Percent  passing  2  micron  (%) 

vs= Percent  passing  #200  sieve  (%) 

yds= Dry  density  (pcf) 

Refer  to  Appendix  B-1 20  for  program 
design  layout. 


(Boring  sample 

counter) 

b  =  0 


b  =  b  +  1 
s  =  0  (layer  counter) 


Input  location  of  boring, 

number  of  soil  layers,  and 

depth  of  boring  sample. 

v_ 

Yb=_ 

Mb= (number  of  layers) 

BDb= (depth  of  boring 

sample) 

Refer  to  Appendix  B-1 19  for 

program  design  layout. 


Yes 


Calculate  additional 

soil  property  quantities 

required  to  make 

additional 

computations. 


Give  user  the  option  to  view 

spreadsheet  of  other  calculated 

material  properties  used  in  volume 

strain,  suction  and  water  flow 

velocity  calculations.  Refer  to 

Appendix  B-1 26. 
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Moisture  Effect  Input 


Select  Moisture 

Case:  Refer  to  Cut-out 

2  below  for  possible 

selections,  or 

return  to  start 


Example:  Cut-out  2  from  Appendix  B-121 


Bmmnaaa 


Program  Design  Layout 


■-im>i 


Sutton  PrortU*    I 


Y<w* 


J 


Vertical 

^^  Barrkr 
(lun  Uhng  Oh  Mi 
hrnfth  pf  the  KhI  edge  of 


Example:  Cut-out  3  from  Appendix  B-125 


r 

Start  Location  of  Effect 

End  Location  of  Effect        j 

X  Coord       l       YCoor<i 

X  Coord              Y  Coord. 

I                      HI                       »   | 

0    II                       50 

Blili^HMta^HianBMBMBMtoMMlltfitaWMIItttt^HraBMhK^H^niM 

Example:  Cut-out  4  from  Appendix  B-121 


As  the  Moisture  Effect  Case  types 

are  selected  a  database 
spreadsheet  row  is  generated  anc 
a  counter  is  initiated  for  each  case 


Example 


For  example,  the  first  "flowerbed"  moisture 

effect  case  selected  would  be  given  the 
"case  number"  name  (1-6)  [and  the  second 
would  be  (2-6)  and  so  on].  A  row  would  be 
established  to  contain  all  pertinent  cells  of 
information  required  for  computations  made 
by  this  software  package  for  that  particular 
case.  Refer  to  cut-out  7  on  the  next  page. 


Prompt  user  to  input  a  "start"  and  "endj" 
location  for  all  moisture  effect  cases 


Example 


For  example,  the  "start"  and  "end"  points  for 
the  flower  bed  are  (40,0)  and  (0,50) 
respectively.  If  the  moisture  effect  is  a  point 
source  then  the  two  inputs  will  be  the  same. 
Refer  to  the  cut-out  3  and  cut-out  4  for 
location  of  flowerbed  on  grid  and 
example  of  coordinate  inputs  respectively. 


Moisture  Effect  Input 
Continued 
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Hf  asure  d  Suction  Profile 
Input 


Select*  Input  Moisture  Effect 

Cam 


Select  Type  of  Moisture  Effect  Case 
Start  Location  of  Effect  End  Location  of  EfTect 

X  Coord.      |       Y  Coord.  X  Coord.      Il     Y  Coord. 


40  II  0 

effect  tooted  el  edfe  a*  bendy 
Yes 


Depth  of  Eflrci 


'No 


I  0  II  50 

input  Eitge  Mois  ture  Variation 
Edge  Ltd  I  Center  LU1 


40 


"bar 


"Izl 


ilfefxtivt  nviib*rd»ifl*e.tM  moUmr%  «tttct  *4f%  t«b%  XAdtrtfc*  fcoinick.) 


<     On 


i  nil  it  tuitict 

>  nil  with  hoiizonttlbirrui 

»  toil  with  venial t» m#i 


6  Qnuualh  honzc»tilb«mer 

6  Gn  ii  with  venial  biniM 

f  Tree  mii  kvMnida 

9  Tree  with  verticil  b&mer 

9  Flowerbed  neir  toundittoaedp 

10.  Flower  b«d  uathhoriiontilbuil*! 

11  Flowerbed  with  remcUbtnui 

12  Meisured  sur !  10  il  pro  file  i 


Perform  the  following 

function  for  moisture  effect 

cases  7  thru  1 1 . 


Example:  Cut-out  5  from  Appendix  B-121 


For  example,  the  depth  of  the  flower  bed  in 

this  case  is  4  feet.  Refer  to  cut-out  5  shown 

above. 


Example 


Prompt  user  to  input  the  depth  of  the 
moisture  effect  zone. 
d„„  = 


Refer  to  cut-out  6  for  example  of  edge 
location  input  for  the  flower  bed. 


Is  effect  located  at  edge  t>i  foundation?! 

Yesor 


Depth  of  Effect 


4J0 

If  Ne",  inyut f 

distance  from  edge. 
(Ne gative  number  defines  the 
moisture  effect  edfe  to  be  under  the 
foundation.) 


I  0J  ht.) 


Prompt  user  to  enter 

location  wrt  slab. 

(-)  number  indicates 

the  zone  is  under  the 

slab,  (+)  number 
indicates  the  zone  is 
away  from  the  slab. 


When  yes  is  selected 

this  indicates  the  effec 

is  at  the  edge  of  the 

slab. 


Example:  Cut-out  6  from  Appendix  B-121 


Provide  the  user  opportunity  to 
review  spreadsheet  of  moisture 
effect  input  data.  Refer  to  cut- 
out 7. 


X  Coord 

em 

Coord 

Barrier 

Size 

Orplh 
(B.f 

Ed*      1 

Umlwnl 

I  Com* 

Y Coord 

X  Coord  j  t 

!    1-9 

0 

56 

65     | 

120 

4 

o 

o    i 

!     1-7 

0 

75 

0      | 

no 

20 

-73     | 

j  1J, 

40 

0 

0      | 

50 

4 

o    1 

1  « 

65 

120 

65     | 

0 

4 

0     | 

1    1.12 

65 

100 

1 

|    2-12 

66 

20 

13-12 

45 

60 

1                1 

I 

Example:  Cut-out  7  from  Appendix  B-121 
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Measured  Suction 
Profile  Input 


Input  number  of  Suction 

Measurements  (including 

the  surface  value) 

N  = 


Initialize  Suction 
Measurement  Counter: 

c=1 


Input  Suction  Value 
@  the  surface. 

hc  = (PF) 


Note:  Measurement  one  (h^ 

must  always  be  at  the  surface  and 

D1  is  always  0.    If  no  value  is 

entered  a  default  value  ofhm 
(equilibrium  suction,  calculated 

on  pg  90)   will  be  assumed. 


c  =  c+  1 


-Yes 


Input  the  Depth  to 

Measurement  "C". 

Dc  = (cm) 


Input  the  suction  value 

@  Measurement  "C". 

hc= (pF) 
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Soil  Property 
Calculations 


Input  Quantites 
(Retieve  values  from  pg  3) 

LL(zj)  =  LL, 
PL(z0  =  PL, 

H(Z|)  =  Hs 

v(^)  =  vs 
Vd(^)=f()s 


Initialize  Variables 

Mb  (retrieved  from  pg  84) 

i  =  0  (increment  counter) 

increment  =  5  (cm) 

s  =  0  (layer  counter) 

z.  =  0  (depth  to  the  i  *  increment) 

^(z.)  =  ° 

em.(z,)  =  0 


i  =  i+  1 


s+  1 


Plasticity  Index  (in  percent) 

Pl(z,)  =  LL(zi)  -  PL(20 


Percent  Fine  Clay  (in  percent) 

H(z.) 


K 


Activity  Ratio 


Ac(z,)  = 


PI(z,) 
Clay(Zi) 


Cation  Exchange  Capacity  (in  percent) 

The  Larger  of  the       CEC(z  )  =  PL(z)117 


two  quantities 


CEC(z,)  =  LL(z)09'2 


Cation  Exchange  Activity 
„.  ,  ,  CECfz,) 
CEAC(Z')  =  ^S) 


Slope  of  the  Suction  vs  Gravimetric 
Water  Content  wet  Line 


Yes 


Matric  Suction  Compression  Index 

Use  calculated  values  of  Ac(z  i  )  and 

CEAc(zi)  and  read  values  of  y100  from  the 

chart  in  Figure  15  developed  by  McKeen. 

Yh(z,)  =  Clay(zi)yU)0 


Mitchell's  Diffusion  Coefficient  (cm   2/s) 

a(z,)  =  0.0029-0.000162  S(z,)- 0.0122  yh(zr) 


Mitchell's  Unsaturated  Permeability 

Ko(zi)-|ho|  =  0.4343-°C(Zi)Yd(Zi) 
S(z,)-yw 


Edge  Moisture  Variation  Distance 
Center  Lift 

UseTMl  (input  from  pg  83)  and  a^) 

(calculated  above)  and  go  to  Figure  16 

to  find  the  appropriate  edge  moisture 

variation  distance  for  center  lift  [e  mc(z,)]. 


Edge  Moisture  Variation  Distance 
Edge  Lift 

Use  TMI  (input  from  pg  83)  and  af^) 

(calculated  above)  and  go  to  Figure  1 7 

to  find  the  appropriate  edge  moisture 

variation  distance  for  edge  lift  [e  me(zj)]. 


eme=eJzi) 

em.=em.(Z.) 
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Calculations  for 

Post  Equilibrium 

Moisture  Effect 

Cases 


-Select  1 


Calculate  Equilibrium 

Suction  for  Bare  Soil 

and  Soil  with 

Vegetation  @  the 

surface. 


Calculate  Depth  to 

Constant  Suction(Zm) 

for  all  Moisture  Effect 

Cases. 


Calculate  Suction 
Profiles. 


Calculate  the 

Differential  Soil 

Movement  for  all 

Moisture  Effect  Cases 


Calculations  for 

Post  Construction 

Moisture  Effect 

Cases 


Select : 


Input  Construction 

Time  Constraintsand 

establish  Time 

Baselines 


Calculate  the  Suction 

Profiles  for  all  Moisture 

Effect  Cases. 


Calculate  the 

Differential  Soil 

Movement  for  all 

Moisture  Effect  Cases. 


Equilibrium  Suction 


m 


90 


0 


Input: 

1 .  Input  dam  (units  of  cm), 

2.  Input  hw|  (units  of  pF), 

3.  Input  hdry  (units  of  pF). 


If  input  values  are  not  known,   then 
the  following  default  values  will  be 
used; 

dam  =  30  cm 

hwl  =  2.5PF 
hdry  =  6.0pF  (bare  soil) 
h.    =  4.5pF  (vegetation  @  surface) 

0S  =  0.55 

Gr=0.04 

A  =  160 

B  =  0.68 
If  the  USC  is  known,  then  go  to  Table 
2  to  find  the  corresponding  values  of 
es,  9r,  A,  and  B. 


Retrieve  TMI  from 
pg83 


Calculate    0wl: 


e„  +  ! 


e.,  =  - 


i+! 


Calculate    6 


dry' 


e.+^e. 


6  dry   ~ 


1  + 


iliy 


Calculate  constants: 

T  =  TMI  +  60 

Y  =  0.039337^ +  1.357033 

d,  =  0.4490794m  +  0.304560 

T,  =  0.062651 4m  +  59.53593 


hm  =  hm  +  .01      k- 

m  m 


hm  =  hm  -  .01 

m  m 


Calculate  0m  &9'm: 


e.+Hr-e, 

A 


1  +  1 


1  + 


d„-d, 


Initialize  l\,  &  8: 
6=1 


s^.-e1. 


Depth  to  Constant 
Suction 


m 


0 


The  algorithms  used  to  compute  the 

depth  to  constant  suction  may  vary 

depending  on  the  moisture  effect 

case  being  considered. 
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Depth  to  Constant 

Suction 

For  Bare  Soil 


Depth  to  Constant 

Suction 

For  Grass  @ 

Surface 


Define  Variables 

Uo=  hdry(bare  SOil)  -K\       (Pp) 

s  =        layer  number  considered  given  an 

initial  value  of  zero.  s=0 
zs  =       the  vertical  thickness  of  layer  s  (cm) 
zm  =      depth  to  constant  suction  given  an 

initial  value  of  zero.  zm=0 
Ls=        depth  to  layer  s  (retrieved  from 

pg84) 


Define  Variables 

Uo=       hdry(vegetation)  -  h.,     (pF) 

s  =        layer  number  considered  given  an 

initial  value  of  zero.  s=0 
zs  =       the  vertical  thichness  oflayer  s  (cm) 
zm  =      depth  to  constant  suction  given  an 

initial  value  of  zero.  ^=0 
Ls=        depth  to  layer  s  (retrieved  from 

pg84) 


►      s  =  s  +  1 


Zs=Ls-Zm 


Solve  for  z^, 
.2A 


■al   =  -  H    L 


1 


nrt 
a. 


-Yes 


Solve  for  Ps(the  ratio  of  the  suction 

difference  between  the  top 

increment  and  bottom  increment  foi 

layer  s). 

n  =  Weather  cycles  per  second. 
Mitchell's  diffusivity  for  layer  s. 


P  =e 


Calculate  L^.  The 

difference  between  the 

maximum  and 

minimum  suction 

values  at  depth  l^. 


Depth  to  Moisture  Active 
Zone 
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Depth  to  Constant 

Suction 
For  Tree  (5)  Surface 


Define  Variables 

hdry  (vegetation)  -  hm     (pF) 

layer  number  considered  given  an 

initial  value  of  zero.  s=0 

the  vertical  dimension  of  layer  s  (cm) 

depth  to  constant  suction  given  an 

initial  value  of  zero,  z  =0 

input  value  for  depth  of  root  zone 

(retrieved  from  pg  86) 

depth  to  layer  s  (retrieved  from 

pg84) 


Depth  to  Constant 

Suction 
For  Flower  Bed  @ 

Surface 


Define  Variables 

K-K    (pF) 

layer  number  considered  given  an 
initial  value  of  zero.  s=0 
the  vertical  dimension  of  layer  s  (cm) 
depth  to  constant  suction  given  an 
initial  value  of  zero,  z  =0 

m 

input  value  for  depth  of  flower  bed 
zone  (retrieved  from  pg  86) 
depth  to  layer  s  (retrieved  from 
pg84) 


*•      s  =  s  +  1 


z„  =  z  +  z„    ■* 


U  =  U, 


-Yes- 


No 


Solve  for  z. 


^m    uzone 


-Yes 


Depth  to  Moisture  Active 
Zone 


Solve  for  Ps  (the  ratio  of  the  suction 
difference  between  the  top  and 
bottom  increments  for  layer  s). 

n  =  Weather  cycles  per  second. 
as  =  Mitchell's  diffusivity  for  layer  s. 


P  =e 


difference  between  the 

maximum  and 

minimum  suction 

values  at  depth  Ls. 

UU  =  P.U0 
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Depth  to  Constant 

Suction 

For  a  Measured 

Profile 


s  = 

Zs= 

Zm  = 

K= 

N= 
DN= 

u„= 

Zo  = 


Define  Variables 

layer  number  considered  given  an  initial  value  of  zero.  s=0 

depth  to  layer  s  (retrieved  from  pg  84) 

depth  of  increment  s  (in  cm), 

depth  to  constant  suction  given  an  initial  value  of  zero.  ^=0, 

the  clh  suction  value  for  a  sample  (retrieved  from  pg  87) 

Number  of  suction  measurements  per  sample  (retrieved  from  pg 

depth  to  the  N  lh  suction  measurement  (retrieved  from  pg  87), 

=equilibrium  suction  at  depth  Q  (retrieved  from  pg  96) 

absolute  value  of  (h  N  -  hm(DN)) 

D„ 


Zm=Zm 


-Yes- 


zm=zm  value  for  the  surface 
conditions  most  similar  to 
the  surface  conditions  for 
the  measured  profile  but 
not  less  than  z„  for  bare 
soil. 


is  zm>=zm  value 
of  the  surface  conditions  most- 
similar  to  the  surface  conditions 
vfor  the  measured  profile  but  no> 
sless  than  z„  for  bare  soil?' 


87) 


Yes 


L^-zo 


s  : 

=  s  +  1 

Zm  = 

Ls+Zm 

t 

Zs=Ls-Ls-1 


No 


Depth  to  Moisture  Active 
Zone 

7       —7       +7 

"■m         m       Tinal 


Solve  for  ^,nal 
,0.21 

Zfna^-'nl 


Solve  for  Ps  (the  ratio  of  the  suction 
difference  between  the  top  and  bottom  of 
soil  increment  s) 

n  =  Weather  cycles  per  second. 

as  =  Mitchell's  diffusivity  for  increment  s. 


Calculate  L^.  The 

difference  between  the 

maximum  and 

minimum  suction 

values  at  depth  l^. 

uu=p.u0 
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Calculate  Suction 

Profiles 

Menu  Page 


The  algorithm  used  to  calculate  a 

suction  profile  wrt  depth  varies 

depending  on  the  moisture  effect 

case  and  surface  conditions. 


Go  to  the  Propersuction  profile 

algorithm  depending  on  which  moisturjs 

effect  case  is  being  considered. 
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Equilibrium  Suction 
Profiles 


Bare  Soil 
@  Surface 

Retrieve: 
hm  and  zm  previously 

calculated  by  the 

processes  shown  on  pg 

90  &  92  respectively. 


Grass 

@  Surface 

Retrieve: 

hm  and  zm  previously 

calculated  by  the 

processes  shown  on  pg 

90  &  92  respectively. 


Tree 
@  Surface 

Retrieve: 
hm  and  zm  previously 

calculated  by  the 

processes  shown  on  pg 

90  &  93  respectively. 


Flower  Bed 

@  Surface 

Retrieve. 

hm  and  zm  previously 

calculated  by  the 

processes  shown  on  pg 

90  &  93  respectively. 


The  process  to  calculate 
hm(z,)  is  the  same  for  each 
case.  The  only  difference 

is  the  hm  @  the  surface 
and  the  z^  for  each  case. 


Convert: 

Suction  in  pF  to 

centimeters. 


hm(cm)  =  10 


hm(pF) 


Initialize  Variables: 

i  =  0 

zo  =  0  (cm) 

increment  =  5  (cm) 


-No- 


Yes- 


i  =  i  +  1 

z,  =  zj.1  +  increment 

hm(zi)  =  hm(zn) "  increment 
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Dry  Limit  Suction 

Profile 

For  Bare  Soil  & 

Grass  (a)  Surface 


Bare  Soil 

@  Surface. 

Retrieve  Values  Calculated  in 

previous  process  flow  sheets: 

z,,,  from  pg  92, 

hm  from  pg  90, 

hjz,)  from  pg  96, 

n  from  pg  83, 
a(Zj)  from  pg  88 


Grass 

@  Surface. 

Retrieve  Values  Calculated  in 

previous  process  flow  sheets: 

^  from  pg  92, 

hm  from  pg  90, 
hjzj  from  pg  96, 

n  from  pg  83, 
a(z,)  from  pg  88. 


Input  Dry  Suction  Value  @  the  surface. 

hmax  = (PF) 

h      =  hd    =  6.0pF  (bare  soil  default) 


hma.  =  h. 


dry 


dry 


4.5pF  (vegetation  @  surface  default) 


Initialize  Variables: 

i  =  0  (Increment  counter) 

increment  =  5  (cm) 

z,=  0  (Depth  to  lh  increment) 

KJQ'K*  (Input/default) 

s  =  1  (Layer  counter) 


Retrieve 
The  following  Values 

Mb  from  pg  84 
L  from  pg  84 


i  =  i  +  1 
Zj  =  zl_1  +  increment 


-Yes- 


Uo  ■  hmax(^-l)  "  hm(2|.l) 
S  =  S+  1 


Calculate  Maximum  Suction  at  depth  z,: 


hmax(zi)  =  hm(z,)  +  U0e 


n-7t 
a(z,) 


Yes- 
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Tree  @  Surface. 
Retrieve  Values  Calculated  in 
previous  process  flow  sheets: 

zm  from  pg  93, 
hm  from  pg  90, 

hm(zi) from  P9  96- 
n  from  pg  83, 

a^)  from  pg  88. 


Dry  Limit  Suction 

Profile 

For  Tree  and  Flower 

Bed  @  Surface 
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Flower  Bed  @  Surface. 
Retrieve  Values  Calculated  in 
previous  process  flow  sheets: 

zm  from  pg  93, 

hm  from  pg  90, 
UJz)  from  pg  96, 

n  from  pg  83, 
a(zi)  from  pg  88. 


Input  Dry  Suction 
Value  @  the  surface. 

KM  = (PF) 


dry 


4.5pF 


(Default) 


Retrieve 
The  following  Values 

Mb  from  pg  84 
Ls  from  pg  84 

dzone  fr0m  P9  86 


Initialize  Variables: 

i  =  0  (Increment  counter) 

increment  =  5  (cm) 

Z|  =  0  (Depth  to  ?h  increment) 

hmax(zi)  =  "max  dnput/default) 

U„  =  hmo.  -  h„ 

o         max        m 

s  =  1  (Layer  counter) 
LMb+i  =  °° 


Uo  ■  hmax(ZM)  "  hm(ZHl) 


*-top       T-1 
S  =  S  +  1 


Calculate  Maximum  Suction  at  depth  z,: 


hmax(zi)  =  hm(zi)  +  U0e 


"(Zi    ^'lalz,) 


-Yes 
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Wet  Limit  Suction 

Profile 

For  Bare  Soil, 

Grass,  &  Tree  @ 

Surface 


Initialize  Variables: 

i  =  0  (Increment  counter) 

increment  =  5  (cm) 

2j  =  0  (Depth  to  i  *  increment) 

hmin(2i)  =  ^^Input/default) 

U0  =  hm(z0)  -  h^ 

s  =  1  (Layer  counter) 


i  =  i  +  1 
z{  =  z,.,  +  increment 


Bare  Soil 

@  Surface. 

Retrieve  Values 

Previously 

Calculated: 

z,,,  from  pg  92, 
hm  from  pg  90, 

hm(zi)  from  pg  96, 
n  from  pg  83, 

a(zi)  from  pg  88. 


Grass 

@  Surface. 

Retrieve  Values 

Previously 

Calculated: 

zm  from  pg  92, 
hm  from  pg  90, 

hm(zj)  from  pg  96, 
n  from  pg  83, 

a(Zj)  from  pg  88. 


Input  Wet  Suction 
Value  at  the  surface. 

hmin= (PF) 

"mm  =  \a  =  2-5pF 

(default) 


Retrieve 
The  following  Values 

Mb  from  pg  84 
Ls  from  pg  84 


-Yes- 


Tree  @  Surface. 

Retrieve  Values 
Previously 
Calculated: 

zm  from  pg  93, 
hm  from  pg  90, 

hm(Zi)  from  pg  96, 
n  from  pg  83, 

a(Zj)  from  pg  88. 


U0  =  KfaJ  -  hmin(^-l) 
S  =  S  +  1 


Calculate  Maximum  Suction  at  depth  z  , 

— (  7-  — ' 

hmi„(zi)  =  hm(zi)  +  U0e 


-Yes- 
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Flower  Bed  @  Surface. 
Retrieve  Values  Calculated  in 
previous  process  flow  sheets: 

zm  from  pg  93, 

hm  from  pg  90, 
hj^)  from  pg  96, 

n  from  pg  83, 
a(z,)  from  pg  88. 


Wet  Limit  Suction 

Profile 
For  Flower  Bed  @ 

Surface 


Input  Wet  Suction 
Value  at  the  surface. 

hm,n= (PF) 

hm,n  =  h*  -  2.5pF 
(Default) 


Retrieve 
The  following  Values 

Mb  from  pg  84 
Ls  from  pg  84 

dzone  fr0m  P9  86 


i  =  i  +  1 


Initialize  Variables: 

i  =  0  (Increment  counter) 

increment  =  5  (cm) 

z,  =  0  (Depth  to  ?*  increment) 

IU,W  ■  hmin(lnput/default) 

U0  =  hm(z0)  -  hmin 

s  =  1  (Layer  counter) 

Ztop=0 


U0  =  hm(zH)-hmin(z,1) 


Calculate  Maximum  Suction  at  depth  z. 


hmin(zi)  =  hm(zi)  +  U„e 

mm  V    i  s  m  v     i  /  o 


-U.-z^) 


-Yes- 


Suction  Profile 

For  a  Measured 

Profile 


© 


Retrieve  The  Following  Values 

zm  =  depth  to  constant  suction  (pg  94), 
N  =  number  of  suction  measurements  (pg  87), 

Dc  =  depth  to  measurement  c  (pg  87), 

^cttimN  =  measured  suction  values  from  c  thru  t< 

the  N^  measurement  (pg  87), 

c  =  arbitrary  measurement  counter. 


Initialize  Variables 

i  =  0  (increment  counter) 

increment  =  5  (cm) 

zi  =  0  (depth  to  ith  increment) 

c  =  1  (suction  measurement  counter) 


Initialize  Counter 

c  =  0 


Calculate  Slope 


Slopec+1  = 


*>,+.- h. 


D„ 


D, 


c  =  c+  1 


Yes- 


i  =  i  +  1 
Zj  =  Zj.,  +  increment 


Yes 
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Calculate  Soil 

Volume  change 

Menu  Page 


0 


The  algorithm  used  to  calculate  soil 
volume  change  varies  depending 
on  the  moisture  effect  case  and 
surface  conditions. 


Go  to  the  Proper  soil  volume  change 

algorithm  depending  on  which  moisture 

effect  case  is  considered. 
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Soil  Volume  Change 

for  Bare  Soil,  Grass,  Tree 

and  Flower  bed  at  the 

surface  with 

Vertical  Barrier 


Retrieve  Variable 

dvb«  from  P9  86 


Initialize  Variables 

a,=  40  (cm) 

Az  =    0  (cm) 

i  =  0  (soil  increment  counter) 

z(  =  0  (depth  to  i"1  increment) 

increment  =  5  (cm) 


—j  =-yh(Zi)logJ-— - 


($  -ffl 


Yes 


(AV\      (AV\      (AV}  (AV\      (AV\      ( AV\ 

ItJ.HvJ -M       IvHvHvJ 


3 


lirJ=flvJ. 


Az  =  Az  +  increment 


m 


Initialize  Variables 

ko  =  .33 

f=.5 

h,(z,)  =  hm(z,)  (from  pg  96) 

M2,)  =  h™x<Z,) 
[hmax  (z)  for  bare  soil/grass 

go  to  pg  97] 

[hmax  (z)  for  tree/flower  bed 

go  to  pg  98] 


Initialize  Variables 

ko  =  .67 

f  =8 

"A)  =  nm(z,)  (from  P9  96) 

MZi>  =  hmin(Zi) 

[hmin  (z,)  for  bare  soil,  grass,  and 

tree  go  to  pg  99] 

[hmln  (z,)  for  flower  bed  go  to  pg  1 00] 


i  =  i  +  1 
z,  =  z^  +  increment 


Convert  all  Suction  in  pF  to  centimeters. 
h(cm)  =  10h(pF) 


w = w 

W  =  KJti  wettin9 

W  =rW2i)  dryjn9 

is  z,  <=  40cm?j^:^Yes 


goto 

Vpgioy 


; 


o-F(z,)  =  0 


h(z,)'ogi. 


hF(z,) 
h,(z.) 


Soil  Volume  Change 

for  Bare  Soil,  Grass,  Tree  and 

Flower  Bed  at  Surface  with 

Horizontal  Barrier 
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Retrieve  Values 

Note:  hoizontal  barrier  can  no 
be  used  for  a  tree. 

dhb«  from  P9  86 
d         from  pa  86 

source  Ka 

Default  values  are: 


Initialize  Variables 

i  =  0  (soil  increment  counter) 
Z|  =  0  (depth  to  lh  increment) 

increment  =  5  (cm) 
emc  =  center  lift  (from  pg  88) 
>xit(zi)  =  hmax(zi)  (P9  97  for  grass  &  bare  soi|) 

KJZ)  =  hmax(Zi)  (Pg  98  for  fl0Wer  bed) 

hm(zi)  =  equilibrium  suction  (from  pg  96) 

edge  ~     mc  "     hbar 


Initialize  Variables 

i  =  0  (soil  increment  counter) 
z,  =  0  (depth  to  lh  increment) 

increment  =  5  (cm) 
eme  =  edge  lift  (from  pg  88) 

exit(Zi)  =  hmin(Zi)  (Pg  "  for  graSS  &  bare  Soil) 

he»t(zi)  =  hmin(zi)  (Pg  10° for  flower  bed) 
hm(zi)  =  equilibrium  suction  (from  pg  96) 

edge  ~~     me  "     hbar 


Convert  all 
Suction  in  pF  to  centimeters. 


h(cm)  =  10h<pF) 


Calculate  Horizontal  Velocity 


Vx(Zi)  = 


K0(Zi)-|h0| 

\    mc/e  source/ 


In 


h,(z.) 

KM 


I 


Calculate  Suction  at  Foundation  edge 


hedge(Zi)  =  hm(Zi)-e 


-vx(zi)dedee 

K0{2j)-|h0| 


i  =  i  +  1 
z{  =  zj_1  +  increment 


Yes 


Use  h^   (z,)  in  the  place  of  hF(Z|) 

on  pg  103  to  make  volume  change 

calculations  for  bare  soil  at  the  surface 

with  a  horizontal  barrier. 
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Soil  Volume  Change 
for  Measured  Suction  Profiles 


Goto  pg  101  to  obtain 

values  of 

h,^)  and  h>(zj). 


Mo 

<r      \~  h   - 

hFf>- 

Yes 

1 

(Shrink) 

r 

(Heave) 

" 

Initialize  Variables 

k0  =  .33 

f=.5 

t\JLz)  =  hm(z,)  (from  pg  96) 

hF(z,)  =  hmax(z,) 
[hmax  (z;)  for  bare  soil/grass 

go  to  pg  97] 

[hmax  (z)  for  tree/flower  bed 

go  to  pg  98] 

Initialize  Variables 

k0  =  .67 

f=8 

h,(z,)  =  hm(zj)  (from  pg  96) 

hF<Z,)  =  hm,n(Z,) 

n  (z,)  for  bare  soil,  grass,  and 

tree  go  to  pg  99] 
(z:)  for  flower  bed  go  to  pg  1 00] 

' 

' 

Initialize  Variables 

c,  =  40  (cm) 

Az=    0(cm) 

i  =  0  (soil  increment  counter) 

z,  =  0  (depth  to  i,h  increment) 

increment  =  5  (cm) 

/36 

(  go 

\P9 

(a)\ 
to 
106/ 

Soil  Volume  Change 

for  Measured  Suction  Profiles 

(Continued) 
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fAHA 
Az  =  Az  +  increment  ■  1  — —  1 

t 

i 

(f),=<f), 

ffl  =(f)  At) 


Yes- 


i  =  i 

+ 1 

Zj  =  zi.1  +  increment 

v 

Convert  all  Suction  in  pF  to  centimeters. 
h(cm)  =  10h(pF) 

r-No— <fis  z, <=  40cm?J^^Yes 


o,.(z,)  = 


1  +  2-ko 


aF(z,)  =  0 


AV1  i    m        fhrCZi) 


fAV^  faF(z,; 


Yes 
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Post  Construction  Time 

Constraint  Input  and 

Calculations 


'dry 
const 


Definition  of  Variables 

Driest  month  of  the  year. 
The  month  the  slab  was  constructed. 
The  month  the  soil  volume  change  is  sought 
=Time  elapsed  between  the  driest  month  and 
and  the  month  of  construction. 
'vc/const =  Time  elaPsed  between  the  months  of 

construction  and  the  month  the  soil  volume 

change  is  sought. 

Time  elapsed  between  the  driest  month  and  th£ 

month  the  soil  volume  change  is  sought. 

Dummy  variable 

Dummy  variable 


const/dry 


K/c/dry 

test1  = 
test2= 


Input  the  driest  month  of  the  year 

Input  the  month  of  construction 

Input  the  month  soil  volume  change  is  sought 


'dry  • 
const " 
*vc=- 


Set  text  values  equal  to 
numberical  values  fortdry 


andtvc 


"Jan"=0 
"Feb"=1 
"Mar"=2 
"Apr"=3 
"May"=4 
"Jun"=5 


"Jul"=6 

"Aug"=7 

"Sep"=8 

"Oct"=9 

"Nov"=10 

"Dec"=11 


tconsVdry   =  ^St1  +  test2 


teStlHonst-t, 


onst       dry 


test2  =  12 

^Y 

tvc/cons.  ■  testl  +  test2 

Yes 

r 

test2  =  0 

S 

T(c/dry       Tic/const         const/dry 

' 

' 

Convert  Time  in 

/38 

Sec 

ands 

\pg89 

test1=t^-  t0 
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Calculate  Post 

Construction 

Suction  Profiles 

Menu  Page 


@ 


The  algorithm  used  to  calculate  a 

suction  profile  wrt  depth  and  time 

varies  depending  on  the  moisture 

effect  case  and  surface  conditions. 


Go  to  the  Proper  suction  profile 

algorithm  depending  on  which  moisture 

effect  case  is  being  considered. 
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Retrieve  and 
Intaialize  Variables 
for  Suction  Profile 

Algorithms 


0 


Bare  Soil  at  Surface. 

Retrieve  Values  from  previous 

process  flow  sheets: 

tyc/dry  (SeC)  fr0m  P9  107, 

zm  from  pg  92, 

hm  from  pg  90, 

hm(zi)  from  pg  96, 

n  from  pg  83, 


n  from  pg  83, 
a(Zj)  from  pg  88, 
hdry  from  P9  97, 
hm«(zi)fromP9  9 
hwl  from  pg  99, 
Mb  from  pg  84, 
Ls  from  pg  84, 

*  ~  Sre/dry" 


0 


Grass  at  Surface. 

Retrieve  Values  from  previous 

process  flow  sheets: 

tvC,dry  (sec)  from  pg  107, 
zm  from  pg  92, 
hm  from  pg  90, 

hJ2;) from  P9  96. 

n  from  pg  83, 
a(Zj)  from  pg  88, 
hdiy  from  pg  97, 
hnw»(z,)fromP9  97, 
hwl  from  pg  99, 
Mb  from  pg  84, 
Ls  from  pg  84, 

*  ~  S/c/dry ' 


Tree  at  Surface. 

Retrieve  Values  from  previous 

process  flow  sheets: 

tvc/dor(SeC)frOmP9107, 

2^  from  pg  93, 

hm  from  pg  90, 
hm(zi)  from  pg  96, 

n  from  pg  83, 
a(Z|)  from  pg  88, 

hd    from  pg  98, 

hm«(zi) from  P9  98, 
hw|  from  pg  99, 

hmin(^)from99, 
Mb  from  pg  84, 
Ls  from  pg  84, 

*  ~  \clin- 


0 


Under  Slab  Profile  at  time  of 

construction. 

Retrieve  Values  from  previous 

process  flow  sheets: 

tconst/dry(SeC)frOmP9107, 
—    const/dry 

Use  the  time  variable  listed 

above  (t=tc?n$Wry)  in  the  place  of 

time  variable (t^,.,,,   )  and 

perform  suction  profile 

calculations  depending  on  the 

moisture  effect  conditions  that 

existed  at  the  surface  just  prior  tb 

placing  the  slab. 


0 


Flower  Bed  at  Surface. 

Retrieve  Values  from  previous 

process  flow  sheets: 

t^dry  (sec)  from  pg  107, 

zm  from  pg  93, 

hm  from  pg  90, 

hm(zi)  from  pg  96, 

n  from  pg  83, 
a(z,)  from  pg  88, 
h.    from  pg  98, 

Knty  from  P9  98, 

hwlfrompg  100, 

Mb  from  pg  84, 

Ls  from  pg  84, 

*  "  *Vc/dry  ■ 


Suction  Profile 
Algorithms 


no 


+  1 


Yes 


U.1=hm(ZM>-»U,(ZM) 

Yes-JU^  =  h-(^-,)-hm(^) 

Z.OP  =  ZM 

s  =  s+  1 


Initialize  Variables: 

i  =  0  (Increment  counter) 

increment  =  5  (cm) 

z,  =  0  (Depth  to  ith  increment) 

s  =  1  (Layer  counter) 
z,oP  =  n 


45 


No 


For  tree,  if  h(zl,t)<hm(zl) 

then  U0=U01 

else  all  others 

U=U„, 


U0  =  hmax(z,,)-hm(z,,) 


U„  =  hm(^)-hmin(zM) 


^op       ZM 


Yes 


1            /h(z,.„t)-hm(z,.,) 
t"-  =  2^'aC0l Vm 


No- 


Calculate  Maximum  Suction  atzand  t 

h(zi,t)  =  hm(zi)  +  U0cos^27tnt-(zi-z,op)J^ 


1> 


Calculate  Maximum  Suction  at  zand  t: 


h(zl,t)  =  hm(z,)  +  Ul,e 


cos  2it-nt-(zrz    ) 


a(z,) 


U  =Uo1 

O  01 


No 


3- 


-No- 


hsoil(zj,t)=h(zi,t)  for  bare  soil, 

hgrasS(Z,'t)=h(Zlt)  for  9raSS- 

fVjz.t^h^.t)  for  tree, 
hfb(zj,t)=h(zi,t)  for  flower  bed, 
hs|ab(zj,t)=h(z],t)  for  under  slab, 


"Damped"  Suction 

Profile  Underneath 

the  Slab 

nslab'Zi'S/c/dry/ 


Initialize  Variables: 

i  =  0  (Increment  counter) 

increment  =  5  (cm) 

z.=  0  (Depth  to  lh  increment) 


i  =  i  +  1 
z,  =  zM  +  increment 


,      v        klab(Zi'tco,B,/d^)-hm(Zi)] 

U,(z;)  =  1 l 


Yes 


a  = 


U1(zi)-hm(zi) 


hslab(Zi'tconst/dry)-U|(Zi) 


t,(Zi)' 
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Retrieve  Values  from  previous 
process  flow  sheets: 

hJZj)  from  pg  96, 

hslab(ZAons./dJ  fr°m  P9  109/1 10, 

a(Zj)  from  pg  88, 
n  from  pg  83, 

K„(z,)|ho|from  pg  88 
Ls  from  pg  84, 

tcons«d.y(SeC)fr0mP9107' 

tvc/const^e^frompgl07. 

tvc/dry  (SeC)  fr0m  P9  107, 

z,,,  from  pg  92  or  93, 

p  =1  (an  exponent  in  from 

mitchell's  equation), 

t,(Zj)  time  when  the  value  l^(z-() 

is  reached, 

U^zJ  is  a  suction  value  halfway 

between 

hm(Zj)  and  hslab(Zi,tconst/dry) 


t,(zi)  = 


a(Zj) 


n-7i 


ln(0i) 


K0(z;)-h 


0.8686 


NslabVZi'tconst/dry)       hm(Zj) 


^1slab(Zi'';vc/drv)  ~~ 


hm(Zi)  +  hslab(Zi5t 


slab  V^j '  l  const/dry 

1  +  a- tvc/const 


)-a-t 


-p 


vc/const 
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Post  Construction 

Soil  Volume  Change 

Menu  Page 


0 


The  algorithm  used  to  calculate  soil 

volume  change  varies  depending 

on  the  moisture  effect  case  and 

surface  conditions. 


Go  to  the  Propersoil  volume  change 

algorithm  depending  on  which  moistur4 

effect  case  is  considered. 
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Retrive  Values  for  a 

Particular  Case  of  Post 

Construction  Soil  Volume 

Change  Calculations 


Initialize  And  Retrieve  Values 


o,=  40(cm) 

^edge  =  °  <Cm) 
^slab  "  °  (Cm) 

i  =  0  (soil  increment  counter) 

z,  =  0  (depth  to  ith  increment) 

increment  =  5  (cm) 

retrieve  dvbar  from  pg  86 

retrive  dhbar  from  pg  86 

retrieve  hsoil(zi,tvc,dry)  from  pg  1 10 

retrieve  hslab(zi,tvc/dry)  from  pg  1 1 1 

retrieve  hslab(zj(tconst/    )  from  pg  110 

^Mvc/dry^soillMvc/dry) 


Initialize  And  Retrieve  Values 


a.  =  40  (cm) 


Az. 


Az., 


'edge 


0  (cm) 


-s.ab  -  °  (Cm) 

i  =  0  (soil  increment  counter) 

Zj  =  0  (depth  to  i,h  increment) 

increment  =  5  (cm) 

retrieve  dvbar  from  pg  86 

retrive  dhbar  from  pg  86 

retrieve  h9rass(zj,tvc/dry)  from  pg  1 10 

retrieve  hslab(zi,tconst/d  )  from  pg  1 10 

retrieve  hslab(zj,tvc/d   J  from  pg  1 1 1 


"(Mvc/dry^grass^Ac/dry) 


Initialize  And  Retrieve  Values 

o,  =  40  (cm) 

Azedge  =  °  (Cm) 
A^ab  "  0  (Cm> 

i  =  0  (soil  increment  counter) 

z,  =  0  (depth  to  ilh  increment) 

increment  =  5  (cm) 

retrieve  dvbar  from  pg  86 

retrive  dhbar  from  pg  86 

retrieve  dsource  from  pg  86 

retrieve  htns(^.tv^lry)  from  pg  1 10 

retrieve  hslab(zj)tconsl/d  )  from  pg  1 10 

retrieve  h^z,,^, )  from  pg  1 1 1 

"'^•Src/dry  '""treeV^'V/c/dry  ' 


Initialize  And  Retrieve  Values 

a,  =  40  (cm) 

AZedge  =  °  (Cm) 

AZs,,,,  =  0  (cm) 

i  =  0  (soil  increment  counter) 

Zj  =  0  (depth  to  ilh  increment) 

increment  =  5  (cm) 

retrieve  dvbar  from  pg  86 

retrive  dhbarfrom  pg  86 

retrieve  dsource  from  pg  86 

retrieve  h^z^,.,^)  from  pg  1 10 


retrieve  h^z^t,. 


)  from  pg  110 


slab^T  const/diy'  ■■*»■■■  ra 

retrieve  hslab(zi,tvc/dry)  from  pg  1 1 1 

h(Zl-tvc/dry)=htb(Zi'tvc/dry) 
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Post  construction  Soil 

Volume  Change 

Calculations  With  or 

Without  Barriers 

(Part  I  of  II) 


v 

i  =  i+  1 
zs  =  zM  +  increment 

1 

Yes 

_| 

Convert  all  Suction  in  pF 
to  centimeters. 
h(cm)  =  10h<pF) 
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Post  construction  Soil 

Volume  Change 

Calculations  With  or 

Without  Barriers 

(Part  II  of  II) 


i  =  i  +  1 


Convert  all  Suction  in  pF 
to  centimeters. 
h(cm)  =  10h<pF) 


Yes 


i 

hl(Zi)  ■  "slab^i'tccnst/dry) 

► 

hF(Zi)  =  "slab^i'W 

r~ 

(fl  = 

-yh(Z,)iog,„(h;(z;J 

16 


APPENDIX  B 
SOFTWARE  PROGRAM  DESIGN  LAYOUT 
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Program  Design  Layout 


VOLFLO  WIN  1.0 


Emm 


INPUT 


GRAPHIC       OUTPUT 

j£AlJA— AmA^L— . . .  I     ■ .   ■■        ■  ■■  ■ 


INPUT 


GRAPHIC 
IfflgBUML 


OUTPUT 


Project  Info.       Foundation  Layout 
Site  Location 


Moisture  Effect 


Soil  Layer  Prop.  Input 


COMBO  BOX 
LISTING  PROJECTS 


Project  Name; 
Project  Date: 
Project  Engineer: 

Site  Location:  Address: 
City: 


I  Sample 


|4Apr97  J! 

I  Donald  D.  Naiser,  Jr. 


Today's 


.COMBO  BOX 
OF  ENGINEERS 


d* 


J 400  Marion  Pugh  #1905 
)  College  Station 


State:  frx"    Zip:  1 77840     | 


Program  Design  Layout 


18 


INPUT 


GRAPHIC 


OUTPUT 


Project  Info.    I  Foundation  Layout     ]  Moisture  Effect  Input  |   Soil  Layer  Prop.  Input  | 


Foundation  Input 
Instructions 


.^--'Press  this  button  to  see  this  pull 
\c*~~  down  list  of  instructions.  ^.._ 


A  three  dimensional  rectangular  (X,Y,Z)  coordinate  system  is  used  to  define 
foundation  layout.  The  Y  direction  is  defined  by  some  input  angle  wrt  true 
North.  The  Z  direction  is  depth  of  soil  (positive  is  down).  Every  part  of  the 
foundation  structure  must  be  in  the  positive  X,Y  plane.  The  origin  is  defined 
by  Latitude,  Longitude,  and  Elevation. 


Input  Origin:  Latitude 
Longitude 
Origin  Elevation 
Input  Direction  of  "V  Axis 


29°15'N 


96  28'W 


100.00 


4°30' 


Thornthwaite 
Moisture  Index  I      °  °° 
Weather  Cycles  J      TTZ 
Per  Year1 


(ft.) 

Angle  Clokwise  front 
True  North 


Input  Foundation  Corners  with  respect  to  origin 
using  the  established  coordinate  system. 


Corner  #   1  Corner  Description 

(ft.) 
X  Coord. 

Y  Coord 

(ft.) 
Z  Coord. 

1  Origin 

Southwest  Corner 

0.0                  0.0 

100.0 

2               NW  Corner                                0.0                120.0               100.0 

|        3             |NE  Corner                         |       65.0          J     120.0          |     100.0 

4             |SE  Corner                                65.0                   0.0               100.0 

Program  Design  Layout 
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INPUT 


GRAPHIC 
nifiPIiflY 


OUTPUT 


Project  Into.    |  Foundation  Layout     I  Moist\ire  Effect  Input  |  Soil  Layer  Prop.  Input  1 


III 


Boring  Sample  Location 


Soil  Layer  Properties 


Number  of  Boring  Samples:  | 

Input  "4"  then  a  spreads  he 


Input 
Boring  # 


4        f  .    Maximum  Layers  Per  Sample :  | 
?et  with  4  boring  sample  input  rows  pops  down, 


120 


Program  Design  Layout 


HHE3 


INPUT       GRAPHIC 
i  DISPLAY 


OUTPUT 


Project  Info.   I  Foundation  Layout     ]  Moisture  Effect  Input  j  Soil  Layer  Prop.  Input  I    / 


Goto  3-D 
-Isometric  of 
I     Soil  Layers 


/ 

13 
0 

I 

Mi                       1 

1  -2  li    and  Yd 

akiUi ted  Soil 
roperttes 

IhM  hhhhhi  hem 

Note;  Input  must  be  made  for  LL  PL,  -#20C 

t 

Boring  # 

Layer  # 

Depth  to  (ft) 
next  layer 

LL% 

PL% 

% 
-2|A 

% 

#200 

rd(po£) 

*! 

I        * 

1 

2 

56 

22 

25 

65 

120 

I        1 

2 

6 

67 

24 

46 

76 

120       | 

1        * 

3 

11 

44                IS 

18 

47 

120       | 

j         1             |         4 

30           |      53              20 

31 

1       64        | 

120 

2                     1 

4                   56               22 

25 

65 

120 

2                       2 

8            j      67                24 

46 

I  «   i 

120 

|        2            |         3 

11            |       44               18 

18 

1      47        | 

120      I 

|        2            |         4 

30           |       53               20 

31 

64 

120      j 

|         3             |          1 

5            |      56        |      22 

25     J 

65 

120      1 

MK* 

|         3             |          2 

10           |      67        |       24 

46 

76 

120 

|         3                       3 

11                   44                18 

is       ! 

1       « 

120 

|         3                      4 

30           |       53               20 

31 

|       64 

120 

1         4             |          1 

1            |      56               22 

25 

65 

120 

4                        2                        4                 67                24 

46          1 

1      76 

120 
120      j 

4                        3                      11           |       44               18 

18         ] 

47 

4                       4                      30           |       53               20 

31         1 

64 

120      | 

▼I 

Program  Design  Layout 
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INPUT 


GRAPHIC 
DISPLAY 


OUTPUT 


Project  Into,    j  Foundation  Layout    ]  Moisture  Effect  Input  I  Soil  Layer  Prop.  Input 


Select  &  Input  Moisture  Effect 
Case 


Measure  d  Suction  Profile 
Input 


Select  Type  of  Moisture  Effect  Case 

Start  Location  of  Effect  End  Location  of  Effect 


X  Coord. 


Y  Coord. 


X  Coord. 


Y  Coord. 


Is  effect  tooted  at  edge  01  foundation?!  Input  Edge  Moisture  Variation 

Yes  or, 


Depth  of  Effect 


No 


No 

If"  Ho",  input / 

distance  from  ecige. 
(Ne  pttive  number  defines  tie 
moisture  effect  edge  to  be  under  the 
foundation.) 


Edge  Lift 


Center  Lift 


Bare  soil  at  surface 

Bare  soil  with  horizontal  barrier 
Bare  soil  with  vertical  barrier 
Grass  at  surface 
Grass  ujdti  horizontal  barrier 
Gra  ss  with  vertical  barrier 


Tree  near  foundation  edge 
9.  Tree  with  vertical  barrier 

9.  Flower  bed  near  foundation  edge 

10 .  Flower  bed  with  horizontal  barrier 

11.  Flowerbed  with  vertical  barrier 

12.  Measured  suction  profiles 


I  I 


flL) 


Start 


Case  # 


1-5 


1-7 


1-9 


End 


X  Coord 


40 


Y Coord 


56 


75 


X  Coord 


65 


Y Coord 


120 


110 


50 


„  iMoist 

Barrier  Effect 

Size       Depth 

"■■>  I  «.i 


La 


Display  Moiat 
Effect  Ca32 
/-  Summary  Shi. 


0 


20 


Edge 


-7.5 


1-6 


65 


120 


65 


1-12 


65 


100 


2-12 


65 


20 


1        I 


3-12 


45 


60 


^^Horizontal 
Barrier 
<Trae  wV  roots 
under  slab. 


Leaking " 
Pq>e 


Vertical 
Barrier^" 


Live  Oak 

Tree  Removed 


Leaking 
Pg>e- 


Jk>ysx2z&. 
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Program  Design  Layout 


INPUT 


Projecting.    |  Foundation  Layout     |  Moisture  Elfect  Input  |  Soil  Layer  Prop.  Input  j 


GRAPHIC      niTTPUT 
DISPLAY     jUUlrU1 


HlalB 


Select  &  Input  Moisture  Effect 
Case 

Input  Number  of 
Measurements  per  Profile 


Measured  Suction  Profile 
Input 

Input  Measured 
Suction  Profiles 


Note:  Equilibrium  suction  will  be  assumed  as  default  values  for  the 
initial/final  if  nothing  is  entered. 


Number 

of  Suction 

Me  as  urements 

Per  Sample 


Initial 

Suction 

Profile(pF) 


Final 
Suction 
Profile(pF) 


Surface 
Conditions 
(  Veget  ation/B  are) 


Input 

Default 

Input 


HI 


Input 
-y|  Default 
T-l    Input 


Vegetation 
Vegetation 
Vegetation 
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Program  Design  Layout 


INPUT      GRAPHIC 
[display 


OUTPUT 


Project  Into,   j  Foundation  Layout    ]  Moisture  Effect  Input  I  Soil  Layer  Prop,  Input 


Select  &  Input  Moisture  Effect 

Input  Number  of 
Measurements  per  Profile 


Measured  Suction  Profile 
Input 

Input  Measured 
Suction  Profiles 


Case 
Number 

Suction 

Me  as  urement 

Number 

Depth  to 
Meas  urement 
(ft) 

Suction 
Measurement 

<pF) 

Initial/ 
Final 

1-12-F 

1 

0.0 

3.0 

Final 

1-1 2-F 

2 

3.0 

2.0 

Final 

1-12-F 

3 

12.0 

3.0 

Final 

1-12-F 

4 

16.0 

3.1 

Final 

1-12-F 

5 

20.0 

3.3 

Final 

2-1 2-F 

1 

0.0 

3.1 

Final 

2-12-F 

2 

3.0 

1.8 

Final 

2-12-F 

3 

12.0 

3.0 

Final 

2-12-F 

4 

16.0 

3.1 

Final 

2-12-F 

5 

20.0 

3.2 

Final 

3-12-1 

1 

0.0 

4.5 

Initial 

3-12-1 

2 

3.0 

4.5 

Initial 

3-12-1 

3 

12.0 

4.5 

Initial 

3-12-1 

4 

16.0 

4.0 

Initial 

3-12-1 

5 

20.0 

3.3 

Initial 

HlilES 


124 


Program  Design  Layout 


GRAPHIC  DISPLAY 


MHE3 


INPUT 


GRAPHIC     OUTPUT 
DISPLAY     I 


Site  Location  on 

Thornthwaite  Map 


Foundation  Plan       1  3-D  Isometric  of    j  Differential 

Vehv  of  Moist.  Cases         Soil  Layers  Suction  Profiles       SoflMosement 


Project . 


Name:  |    Sample 
Thornthwaite  — — 

Moisture  Index         1 0.0 

..Location  of  Site 
Latitude     1 29°15'  N 
Longitude     |96°28'W 

Edge  Moisture 
Variation  Based  on 

This  Location   

era(Edge  Lift)  |  6.5  (ft) 


6m(Center  Lift)  10.5  (ft) 
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Program  Design  Layout 


GRAPHIC  DISPLAY 


H@B 


INPUT 


GRAPHIC 


OUTPUT 


Site  Location  on 

Thornt hwaite  Map 


Foundation  Plan 


3-D  Isometric  of 


Veiw  of  Moist.  Cases  Soil  Layers 


Suction  Profiles 


Differential 
Soil  Movement 


Y(ft) 


Tree 


pF4.5 


1)0 


30 


60 


40 


Sidewalk 


20 


40 


Leaking 
pF2.0 


Pipe 


Vertical 
Barrier 
along  the  full 
of  the  East  edge  of 
fou  ndation.) 


Leaking  Pip 
pF2.0 


60 


80 


100. 


X(ft) 


126 


Program  Design  Layout 


GRAPHIC  DISPLAY 


INPUT 


GRAPHIC 
DISPLAY 


OUTPUT 


Site  Location  on      J  Foundation  Plan 

Thornthwaite  Map  |  Veiw  of  Moist.  Cases 


3-D  Isometric  of 
Soil  Layers 


Suction  Profiles 


Differential 
Soil  Movement 


Y(ft) 


Layer  1 


Origin 


_-  X(ft) 


Layer  4 


Goto   Soil  Layer 
Property  Input   _ 


z(ft) 


Go  Back  to 

Soil  Layer 
Input 


Summary  of  Soil  Layer  Properties 

Soil 

%LL        %PL       %-#200        %-2/i          %PI    I  %Fine  Clay  j     CEC 

Layer  # 

% 

1 

56.0%       22.0%       65.0% 

25.0%|     34.0% 

38.46% 

37.21 

2 

67.0%      24.0% 

76.0% 

46.0%|    43.0% 

60.53% 

41.20 

3 

44.0%       18.0% 

47.0% 

18.0%j    26.0% 

38.31 

1%       29.42 

1%J      33.28 

4 

53.0%      20.0%      64.0% 

31.0%     33.0% 

48.44 

Soil 

Ac          CEA       yMoo          7h              7d 

S 

a 

Layer # 

%                          lbs/ft3 

cm2/s 

1 

0.8840     0.9674      0.163    0.062692!      120 

-11.11 

0.003936 

2 

0.71041    0.6806      0.163    0.098658)      120 

-9.70 

0.003268 

3 

0.6789     0.7682      0.163    0.062426        120 

-13.28 

0.004289 

4 

0.6813J    0.6871      0.163    0.078953       120 

1 ; 1 -i 

-11.53 

0.003805 



i ; 

6m                       ©m 

Soil 

p             ko 

THI 

n       |  Center  Liftj     Edge  Lift 

Layer # 

cm2/s         cm/s 

cycles/yr          ft                   ft 

1 

0.000681     2.96E06 

0 

1                   10.2                    6 

2 

0.000648    2.81  E -06         0 

1                      7.6                     5 

3 

0.000621     2. 70  E -06         0 

1                    10.5                  6.5 

4 

0.000635    2.76 

E4J6 

0 

1 

9.5 

5.9 

127 


Program  Design  Layout 


GRAPHIC  DISPLAY 


H1S.1E3 


|l^~|  1e  i 


INPUT 


GRAPHIC 


OUTPUT 


Site  Location  on 
Thornthwaite  Man 


Foundation  Plan        I  3-D  Isometric  of 
Veiw  of  Moist.  Cases  Soil  Layers 


Sue  tkm  Profiles 


Differential 

Soil  Movement 


Post  Construction 
Suction  Profiles 


Post  Equilibrium 
Siiction  Profiles 


January 


Select  the  month  the  slab  is  to  tie  constructed. 

Input  the  driest  month  of  the  year.   I  June 

Select  perimeter  moisture  effect  case  |  15  Grass  urfhoriz.bar. 

you  wish  to  see  suction  profile  for 
Select  the  month  you  wish  to  display 
the  suction  profile  for.  - 


Jan    |   Feb    j  Mar   |  Apr    |  May   |  June 


Underslab  Suction 
Profile 


1 b  Gra  is  urf  vert .  bar. 


1 


1-7  Tree 


1-9  Flower  Bed 


1-12  Measured 


2-12  Measured 


Measured 
Suction  Prof. 

I  |  Oct    |  Nov    j    Dec 


■  Suction  Profile  At  Slab  Edge  With  Barrier. 

Suction  Profile  Under  Slab  in  June. 

Suction  Profile  At  Time  of  Construction. 

Equilibrium  Suction  Profile  3.31  pF. 

Maximum  Suction  Profile. 
'  Minimum  Suction  Profile. 


dvbar 


4ft  Depth  of  Barrier 
:  10.3ft  Depth  to  Constant  Suction 
:  10.5ft  Edge  Moisture  Variation  Distance 


128 


Program  Design  Layout 


GRAPHIC  DISPLAY 


HWQ 


INPUT 


GRAPHIC 
DISPLAY 


OUTPUT 


Site  Location  on 
Thornthwaite  Map 


Foundation  Plan        ]  3-D  Isometric  of 
Veiw  of  Moist.  Cases  Soil  Layers 


Post  Construction  Slab  Deflection 


Suction  Profiles 
Post  Equilibrium  Slab  Deflection 


Differential 
Soil  Movement 


Jan  |    Feb  |  Mar]  Apr]  May)  June]  July  |  Aug]  Sept)  Oct  |  Nov|   Dec|  w«f*tCi«.W€t|w«mCi«ihy] 
Plan  View  Heave  Plan  View  Shrink  Profile  View 


129 


Program  Design  Layout 


GRAPHIC  DISPLAY 


Site  Location  on 
Thorathwaite  Map 


INPUT      jjJJSS?    OUTPUT 

■ ■ IM-AImULpJ|MAi»i I  nimmii  «.i*  i,«.n»iii  m,  ...pi 

Foundation  Plan        I  3-D  Isome trie  of 
Veiw  of  Moist.  Cases  Soil  Layers 


Post  Construction 
Suction  Profiles 


Post  Equilibrium 
Suction  Profiles 


Select  perimeter  moisture  effect  case  for  which  you  wish  to  see  suction  profile. 

W<urM  Cmut  W«t  j  Worit  Cut  Dry] 


Section  A-A 


]    15  Grass  urn/  h'jni. bar 

■|    lb  Grass  nrf  vert  bar.        f 

1    IT  Tree 

)   1-9  Flower  Bed 

|    112  Mea  sured 

j  212  Measured 

|  3  12  Measured 

Slab  Surface 


^ 


r                         ,     i    ,    ■ 

».    ►   i« 

in        b       ig 

Suction  (pF)l 

b 
o 

I 

ro 

o 

b 

o 

70     zone 

105       | 

! 

140 

a 

175  » 
1 

210  $ 

m 

/ 

245  "^ 
280  — 

', 

/ 

315 

f 

350 

Ij 

f 

385 

! 

420 

■ 

Suction  profile  due  to  the  flower  bed. 
Equilibrium  suction  profile. 


em  =  6.5  (ft)  Edge  moisture  variation  distance. 

zm  =13.8  (ft)  Depth  to  constant  suction. 

hm  =  3.57 (ft)  Equilibrium  suction  at  the  surface. 


dzone  =  4.0  (ft)  Depth  of  the  flowerbed. 


130 


Program  Design  Layout 


GRAPHIC  DISPLAY 


liQ 


INPUT 


GRAPHIC 


OUTPUT 


Site  Location  on      1   foundation  Plan        j  3-D  Isometric  of    I  Differential 

Thornthwaite  Map    Veiw  of  Moist.  Cases         SoilLayers  Suction  Profiles        sou  Movement 


Post  Construction  Slab  Deflection  Post  Equilibrium  Slab  Deflection 

Jan  )   Feb  |  Mar}  Apr}  May  |  June  |  July  j  Aug  J  Sept}  Oct}  Nov}  Dec}  w«*c««ww|w«*Cm«'Ot^ 
Plan  View  Heave    j      Plan  View  Shrink  Profile  View 


131 


APPENDIX  C 

SITE  PLAN  OF  MOISTURE  EFFECT  CASES 
USED  IN  SAMPLE  CALCULATIONS 


132 


GRAPHIC  DISPLAY 


Program  Design  Layout 


INPUT      GRAPHIC     OUTPUT 
IdisplailJ— 


ilQ 


Foundation  Example 

Used  to  Illustrate  Program 

Design  Layout 


4"30' 


Depth  of  Tree 
Root  Fiber 
is  20' 


4'     Sidewalk 


Point  of 
7.0  Water  Source 


Vertical    Barrier 
(Length   of  East  Edge) 


Live  Ook 
Tree  Removed 


Depth  of 
Flower  Bed 
Zone  is  4" 


Point   of 
Water  Source 


Project  Name:  Sample 

Project  Date:  4  April  1991 

Project  Engineer:  Donald  D.  Naiser,  Jr. 

Site  Location  Address:  400  Marion  Pug'h  #1905 

College  Station,  TX  77840 
Approximate  Location:  Lat.  29  15'  N 

Long.  96  28'  W 


133 


APPENDIX  D 


REQUEST  FOR  PERMISSION  TO  USE  COPYRIGHTED  MATERIAL 


134 


Donald  D.  Naiser,  Jr. 

400  Marion  Pugh  #1905 

College  Station,  Texas  77840 


August  27,  1997 


Permissions  Department 
John  Wiley  &  Sons,  Inc. 
605  Third  Avenue 
New  York,  NY  10158-0012 

REQUEST  TO  USE  FIGURES  COPIED  FROM  SOIL  MECHANICS  1969 

I  am  a  master  of  science  candidate  at  Texas  A&M  University  and  I  will  be  finalizing  my 
thesis  in  the  next  two  weeks.  Our  library  sends  our  theses  to  University  Microfilms  Inc.,  for 
preparation  of  a  microfilm  copy  of  the  document.  As  I  am  sure  you  are  aware,  UMI  retains  the 
right  to  sell  both  hard  and  microfilm  copies  of  the  document. 

My  reason  for  writing  is  to  request  permission  to  use  the  figures  listed  below  exactly  as 
they  appear  in  SOIL  MECHANICS,  by  Lambe  and  Whitman,  copyrighted  by  John  Wiley  & 
Sons,  Inc.,  1969.  With  your  permission,  I  will  use  the  following  figures  in  my  thesis: 

1 .  Figure  4.7  "The  structure  of  serpentine."  Found  on  page  47. 

2.  Figure  4.8  "The  structure  of  kaolinite."  Found  on  page  47. 

3.  Figure  4.9  "The  structure  of  pyrophyllite."  Found  on  page  47. 

4.  Figure  4.10  "The  structure  of  muscovite."  Found  on  page  47. 

5.  Figure  4. 1 1  "Sheet  silicate  minerals."  Found  on  page  48-49. 

6.  Figure  5.9  "Soil  particles  with  water  and  ions."  Found  on  page  55. 

One  chapter  in  my  thesis  describes  the  magnitude  of  the  problem  of  designing 
foundations  built  on  expansive  soils.  The  figures  described  above  will  be  used  to  illustrate  the 
molecular  structure  of  fine  grained  soil  particles.  The  thesis  has  been  accepted  and  the  author 
proofs  indicate  a  1997  publication  date.  The  thesis  title  is  PROCEDURES  TO  PREDICT 
VERTICAI  DIFFERENTIAL  SOIL  MOVEMENT  FOR  EXPANSIVE  SOILS. 

Please  inform  me  of  any  necessary  action  to  prevent  future  publication  complications.  I 
can  be  contacted  by  phone  at  home  (409)  543-8412  or  at  work  (409)  845-9919.  Thank  you  very 
much  for  your  effort  in  resolving  these  issues. 

Sincerely, 


Donald  D.  Naiser,  Jr. 
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VITA 


Donald  David  Naiser,  Jr.  is  a  citizen  and  resident  of  the  United  States  of 
America,  whose  permanent  mailing  address  is, 

Rt.  3  Box  203-C 
El  Campo,  Texas  77437. 

He  obtained  his  Bachelor  of  Science  degree  in  Civil  Engineering,  from  Texas 
A&M  University  in  1991. 

Since  graduating  from  Texas  A&M  University,  he  has  been  commissioned  as  a 
United  States  Naval  Officer  as  a  member  of  the  Civil  Engineer  Corps.  His  time  served  in 
the  Navy  has  been  exclusively  with  the  SEABEES  performing  contingency  construction 
training  as  well  as  constructing  naval  facilities.  His  tour  highlights  have  been  as  the 
Officer  In  Charge  of  the  Elevated  Causeway  System  (ELCAS)  in  Amphibious 
Construction  Battalion  ONE  and  as  the  Officer  In  Charge  of  Construction  Battalion  Unit 
418. 

In  1 996  he  was  accepted  at  Texas  A&M  University  to  pursue  studies  toward  the 
degree  of  M.S.  within  the  Geotechnical  branch  of  Civil  Engineering. 

In  1 997  he  became  registered  as  a  professional  engineer  in  the  state  of 
Washington. 

"Blessed  are  the  people  whose  leaders  can  look  destiny  in  the  eye  without 
flinching  but  also  without  attempting  to  play  God",  Henry  Kissinger.  His  career 
interests  include  achieving  excellence  through  leadership  in  any  Christian  form. 
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